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5 Atomistic and the Continuum Theories | 
14 of Crystal Elasticity We % = 
By R.8S. Krishnan and E.S. Rajagopal 
Summary 
175 The lattice theory of erystal elastieity is developed with special reference 
to the basic postulates of the lattice models, the role of initial stresses and the 
189 validity of the rotational invariance condition. The elastic behaviour of the 


infinite lattice is investigated theoretically by the method of homogeneous 
196 deformation and by the method of long waves. The exact generalization of the 
classical theory, viz., the existence of surface couples, is deduced from these 
considerations. It has also been proved that the presence of volume couples 
220 does not alter the structure of the classical theory. A complete correspondence 
between the lattice and the continuum approaches is exhibited. Some subtle 
considerations, regarding the extension of the continuum concepts of strains 
and stresses to atomic lattices, and the clear distinction between the initial 
stress and the natural stress problems are also pointed out. 


her 

des 

= § 1. Introduction 

ige The present paper aims at a systematic development of the theory of 
lattice elasticity on order to specify precisely the contradictions between the 
current lattice theory and the classical theory of elasticity and to indicate 
the exact nature of the modifications of the theory of elasticity which are 
required to bring about a correspondence between the two view points. This 
study also brings out some features associated with the extension of the con- 
tinuum concepts of strains and stresses to atomic lattices and enables one to 

me analyse critically the recent papers on the theory of crystal elasticity. 

The continuum theory rests on a few axioms about the elastic bodies, the 

validity of which, in the continuum framework, can be inferred only by a 

sie comparison with experiments. The molecular theories offer the possibility of a 

Fy direct justification of some of the basic axioms of the theory of elasticity. 
However, the atomistic theories themselves are based on the general principles 

n. of mechanics and on several physical concepts. Consequently, it is first neces- 

el- sary to discuss some of the basic postulates of the lattice models ($ 2) before 

en proceeding to the derivation of the elastic behaviour of a lattice (§ 3). One can 


then discuss the modifications of the basic axioms of the continuum theory 
which will lead to the correspondence between the lattice models and the conti- 
je- nuum theories ($ 4). On this basis, one can examine the validity of the modi- 
fications which have been (§ 5). 
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The lattice calculations implicitly assume certain factors regarding the 
applicability of the continuum concepts of strains and stresses to the lattice 
models. Since these factors are of great importance, a brief discussion of their 
limitations is presented in Appendices A, B und C. 

Obviously, it is sufficient to elaborate only such points which are not clearly 
or fully brought out in the recent critical reviews of the subject by Trues- 
dell), Born and Huang?), Huntington?) and others, confining onself to 
the basic concepts ans leaving out the details of applications. This is all the 
more necessary since many of the new suggestions suffer from shortcomings 
not in the development of the basic concepts, but in the basic formulations 
themselves. Moreover, it is sufficient to perform all the theoretical calculations 
on the basis of classical mechanics. ' 

‚ A short report of this work was presented at the first Symposium on Solid 

State Physics, at Bangalore, February 1960 (Part I*). Certain features of the 
calculations regarding the role of initial stresses (Part II®), the validity of 
rotational operations (Part III®), and the presence of interfacial couples 
(Part IV’)), have already appeared in this journal. Some remarks on the in- 
consistencies noted by several workers®) on the calculations of lattice elasticity 
by the current methods will be discussed later. 


Zu $ 2. Conceptual postulates of the lattice models 


A convenient starting point for the lattice calculations is the adiabatic 
approximation?)®) which enables one to treat the problem by considering the 
motion of the nuclei in the potential field of the nuclei and the electrons. Ob- 
viously this potential energy will depend upon the spatial co-ordinates and 
the generalized co-ordinates (like spin, magnetic moment, etc.) of the nuclei., 
In the past, it has been tacitly assumed that the effects of nuclear dipoles 
quadrupoles, etc. are entirely negligible in the calculation of the lattice elasti- 
city so that only nuclear (spatial) co-ordinates need be considered. 

The error in this assumption is easily exhibited by calculating the orders 
of magnitude of the quantities involved. Suppose the displacements 

m /m! 

|) = ") Wap = — Wea (1) 

B I 

are impressed on the nuclei. This is equivalent (see Appendix A and also Part 

III) to a rotation of the whole crystals through e = |w.;| followed by a rotation 

of each dipole through e. If the dipole energy is negligible, equation (1) can be 


1) C. Truesdell, J. Rational Mech. Math. 1, 125 (1952); 2, 593 (1953). 

2) M. Born and K. Huang, “Dynamical Theory of Crystal Lattices’? Oxford 1954. 

3) H. B. Huntington, Solid State Physics 7, 213 (1958). 

4) E.S. Rajagopal and R.S. Krishnan, Proc. Ist. Symposium Solid State Phy- 
sics, Bangalore, 1960, Feb. p 66. 


ae 5) E.S. Rajagopal, Ann. Physik 6, 177 (1960). 
ae *) E.S. Rajagopal, Ann. Physik 6, 182 (1960). 
7) E.S.Rajagopal, Ann. Physik 6, 192 (1960). 


8) E. Forslind, Skrift. Svenska Inst. for Cement och Betong, Stockholm, Handlingar 
No. 21, 1954; K. Miyakawa, Physical Rev. 107, 677 (1957); B. Ya. Moizhes, Fiz. 
tverdego Tela 1, 1770 (1959); L. J. Slutsky and A. Livingston, J. chem. Physics 32, 
1093 (1960). ’ 

®) A.Haug and G. Sauermann, Z. Physik 158, 269 (1958). 
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taken as the rigid rotation of the crystal lattice. One must bear in mind that the 
lattice theory of elasticity is strictly developed for crystals at 0 °K (infact, for the 
“static lattice’), where the nuclear alignments are not disturbed by the ther- 
mal energies *). Also, even if the electronic cloud follows the nuclear motion 
instantaneously the nuclear spins etc. relax only with a finite relaxation 
time. 

Engineering strain gauges capable of measuring strains of the order of 
10-7 — 10% cm per cm are available!) and as early as forty years ago Whid- 
dington!!), developed the methods of measuring displacements of 10-* cm. 
So one takes a practical case of “infinitesimal” strain de ~ 10° cm per cm. 
The strain energy associated with the shear deformation of, say, NaF crystal 
(Na® and F!? have almost equal nuclear magnetic moments of ~ 1.2 - 10-3 
erg/gauss as required below) is ~ 0.5 - 2.4 - 104 - (10-9)? = 1.2 - 107 ergs/cm?. 
The nuclear magnetic moment interacts the other nuclei and also with the 
electronic moments. The potential energy of identical dipoles in simple lattices 
has been calculated by Sauer, Luttinger and Tisza) and others. The 
change in energy between the dipoles oriented along [100] and along [110] is 
~ N? u2 per unit volume. (The number of dipoles per unit volume N is ~ 0.8 - 
108 in the present case.) So the energy change associated with an average 
rotation of 10-9 of each dipole is ~ 10? - (1.2 - 108. 0.8 . 10-3)? = 0.9- 109 
erg/cm?. Nextly, the difference in the dipole energy levels of the interactions 
of the nuclear and electronic moments is ~ 1.2 - 10-8 - 0.9 - 102°/(1.5 - 10-8) = 
3.102 erg. So, for the average rotation of each nuclear dipole through 10-°, 
the energy change is ~ 3 - 10-* erg/cm?. One has only to compare these values 
to realize that the directive properties of the nuclei cannot be neglected. 
Indeed, the effect of applied strain on nuclear quadrupole interactions is too 
wellknown®), to warrant any discussion here. 


The displacements of the nuclei produce energies of the same order of 
magnitude since the elastic constants calculated on that basis are of the same 
order as those observed. One may also consider the bending of valence bonds. 
The bending force constants‘) are of the order of ~ 2.104dyne/cm, so that the 
energy change in bending a valence bond through an angle of 10? is ~ 2 - 104 - 
(2.108 - 10-9)? = 1- 10-2 erg. This is to be compared with the value ~ 3 - 10-9 
erg. obtained above for the energy change in rotating the nuclear dipole against 
the part of the potential energy due to the electronic magnetic moments. So 
the valence bond energy cannot be represented in the lattice theory merely by 


9a) See, for example, M. Born, J. chem. Physics 7, 597, 612 (1939); M. Born and K. 
Huang (loc cit) p. 319. 


10) J. Marin, “Engineering Materials”, Prentice Hall, New York 1952, p. 475; J.C. 
Simmons, Mag. Concrete Research 6, N. 19, 34 (1955). The sensitivity can be increased 
easily by two orders of magnitude by the present materials and techniques, W. P. Mason 
and R.N. Thurston, J. Acoust. Soc. Amer. 29, 1097 (1957). 

11) R. Whiddington, Philos. Mag., 40, 634 (1920). 

12) J. A. Sauer, Physical Rev. 57, 142 (1940); J.M. Luttinger and L. Tisza, 
Physical Rev. 70, 954 (1946). 

13) R. G. Shulman, B. J. Wyluda and P. W. Anderson, Physic. Rev., 107, 953 


14) §.I. Mizushima, Handbuch der Physik, XXVI-I, 187 (1958). 


(1957); E. F. Taylor and N. Bloembergen, Physic. Rev., 118, 431 (1959). a ap es eg 
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specifying the angle between the two directions of the bond"). There is a 
direct experimental evidence for this in the ’’slip effect 1*)* which shows that 
when a molecule is rotating the inner shells slip and point to fixed spatial 
orientations just as in the representation of eq. (1). If the energy changes were 
entirely negligible, the effect would not have been found. 

Thus one admits the case that the potential field experienced by the nuclei 
depend upon the orientational properties of the nuclei, so that operations 
which alter the nuclear orientations are to be discussed separately. One may 
define the second derivatives of this potential energy as the force experienced by 
the nuclei per unit displacement. Since the potential energy of an infinite lattice 
(which one usually considers to avoid the delicate conditions at the boundaries 
of the crystal) is unbounded, one often has difficulties in normalizing the 
energy changes to unit volume. Alternatively one may define the force con- 
stants and then build up the total potential. If one is considering two body 
forces alone, as in the present case, the force constant approach has the merit 
of avoiding the difficulties in normalizing the energy changes in the infinite 
lattice. 


m— 
Let Dy, | B ı) be the x component of the force on the (m,) atom due tothe 
I 


ß displacement of the (p,) atom. This ®,, will be arbitrary except for the condi- 
tions imposed by general physical and mathematical requirements®), for 
example 

(a) the invariance of the potential energy under a body translation __ 


m— p 

ae vo, | = 0; S is over all p,) 2 


me (b) the reciprocity of mutual interactions 


m—p p—m 
= 


(c) the of the lattice 
i? 


From these relations by a simple manipulation, one gets | we rie 

x --.-6 

71 


These are the general relations obeyed by the force constants. Only in some 
restricted calculations can one assume the nature of the force field like Coulomb 
field, first neighbour interactions etc. 

15) This factor is not considered in the usual normal co-ordinate treatment of the bond 
beding force constants used by several workers like Harrison, Moizhes and others‘). 
The discordance between their work and the present work may be in part due to this vio- 
lation. The same is the trouble with the atomistic derivations of the Huangs relations by 
G. Leibfried and W. Ludwig, Z. Physik 160, 80 (1960); L. T. Hedin, Arkiv d. Fysik, 
16, 496 (1960); 18, 369 (1960). 

16) C. H. Townes and A. L. Schawlow, "Microwave Spectroscopy“, McGraw Hill, 
New York (1955), p. 17, 194. 
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Apart from these relations (Actually eq. (4) was not used by Begbie and 
Born!?”) and was first proved by Laval!®), Begbie and Born have used 
another relation viz., the rotational invariance condition 


m— p p m—p p 
b = S9@., ° ) 


This is obtained by subjecting the crystal to the displacements of eq. (1) which 
are supposed to represent the rigid rotation of the body. But as seen above, 
this representation fails when the dipolar energies are considered so that the 
derived result eq. (6) is also not valid. A direct proof of the invalidity of eq. (6) 
is given in part III by framing the equations of lattice dynamics in a coordinate 
system imbedded in the crystal so that the translational and rotational in- 
variance are automatically included. The resulting equations agree with such 
of the usual equations which do not make use of eq. (6) and show that the use 
of eq. (6) introduces additional conditions. This implies that eq. (6) is not 
equivalent the physical restriction of invariance under a rigid body rotation. 

An alternate way of expressing the error in eq. (6) is to calculate the diffe- 
rence {2323}—{2332} (eq. 11) for any lattice. For the simple diamond lattice 19) 20) 
restricting oneself to the nearest neighbours only, 


{2332} = (R — R?/Q)/a; {2323} = (Q — R*/Q)/a. 
The bending force Q—R is —!/, of the bond extensional force constant 
Q+2R so that Q~2R and 


{2332}/{2323} ~ 1/3. | 
According to the classical theory, the ratio must be unity. et ar we 


$ 3. The lattice theory of elasticity fi 


The elastic behaviour of the lattice is usually studied by applying a homo- | 
geneous deformation or by treating the propagation of plane waves. Both the 
methods were used by Born?!) using the simplified central force interactions, 
in giving the first comprehensive atomistic description of the lattice elasticity. 
The method of homogeneous deformation has been discarded by some recent — 
workers on the erroneous assumption that the strain energy cannot be normali- 
zed an that the initial state must satisfy certain conditions. As discussed in 
Appendix C, the question of the magnitude of the initial stresses is completely 
irrelevant for the present discussion. Here both the methods will be developed 
briefly since the mathematical manipulations are well known ?)®)’)?!) and atten- 
tion will be paid to the questions of the validity and the interpretation of the 
results. 

It is well known’) that a homogeneous deformation of a lattice must be 
represented by LEN 


B 


17) G.H. Begbie and M. Born, Proc. Roy. Soc. London 188 A, 179 (1947). 
18) J. Laval, J. Physique Radium 18, 247 (1957). 
1) K.S. Viswanathan, Proc. Ind. Acad. Sci. 89 A, 196; 41 A, 98 (1954). 
*) R.S.Krishnan, V.Chandrasekharan and E.S.Rajagopal, Proc. Nat. 
Inst. Sci. India 26 A, 312 (1960). 
2!) M. Born , ”Atomtheorie des festen Zustandes‘, Encyk. Math. Wiss. V—3, (1923). 
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where the n.(j) term denotes the readjustment of the partial lattices. If there 
are s particles in the unit cell, one can determine 3(s —1) quantities 7,(j) 
(j = 2, 3,...8) after setting n.(1) = 0 (i.e., the relative displacements of 
the partial lattices) by forming the condition that the net force on any atom in 
the deformed state is zero. The force on the (m,) atom due to the (p,) atom is 


> 


and so summing over all (p,) 

Forming the “inverse” @ ( Jar of the singular matrix pring 


one gets easily 
na(h) = (M, in Due (8) 


alk ußö Lk 


_ The stress is calculated (Appendix B) by summing the action of all the 
particles of one half of the infinite lattice on all the particles of the other half 
contained in the infinite parallelepiped erected over unit area, a method used 
by Cauchy, Poincare, Born and others. 


1 a\ = 
Otay = 5-3 fa 
where 6t,, is used to emphasize that one can calculate only the incremental 
stress due to the applied infinitesimal deformation ?1%) Substituting eqs.(7) and 
(8) and using the relations (1—4), one obtains 


where the bracket symbols are defined by 
{x 7, = [a y, B 4] + (xy, B 4) 


1 b b hl 

{ = — M M -1/2 
{x y, B ö} ry (, (M,M,) (; (; k 
(1 


1c) 
In the square bracket symbol [x y, 8 6] one can interchange « with 6 and y 
with 6 so that there are 36 terms. In the round brackets (x y, 6 4) the inter- 


aia) Professor A. E. Green has kindly pointed out that tx)is referred to the co-ordinate 
axes Za of the deformed state and to the unit area of the initial state. In such a case the 
natural stress problem (Appendix C) is greatly simplified. See also A. E. Green and J. L. 
Adkins, Proc. Roy. Soc. London (in press). 
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change of xy with £6 alone is permitted leaving 45terms. So the curly brack- 
ets are 45 in number. When each atom is a center of symmetry the round 
brackets vanish and {x y, ß ö} reduce to 36 independent constants with 


{x y, Bd} = {x 6, By} = (By, « = {8 6, ay}. 


Introducing the strain components (eq. 32), and in the absence of the rotational 


components 


where degs is written to emphasize that it is the applied incremental strain. 
One now notices that 


(12) 


bay 2 bye (13) 
(Note an unfortunate slip in eq. (IV—9a) where an inequality sign must be 


present. ) 
One can easily calculate the strain energy density, 


Sine (" is the force constant, 
m — m — —m 
2 i +&k kj 


is the work done in the relative displacement of the (m,) and (p,) atoms. One 
sums up the work done by all the atoms of a unit cell against the displacements 
of all the atoms of the lattice. This divided by the volume of the unit cell gives 


directly the inremental strain energy density. So rt; 


W=-— ®, u 14 
Substituting eqs. (7, 8), and eq. (11) to rearrange the terms, one has 


W => 7, BO} Uny (15) 


This clearly shows that the iid of W leads exactly to the same result 
as the calculation of dt,, as expected. (Note that eq (IV—10) is strictly 
öW = X {a y, Öug,s the notation there being quite ambiguous.) 

The wave propagation ®) is treated by the usual long wave method. One 
starts from the equation of motion 


M;u,(” 


with the harmonic solutions 


The perturbation procedure of solving of solving this equation has been dis- 
cussed in great detail by Born and Huang?) and the method has been used by 
Laval!®) and Viswanathan”) also. Following exactly the same prodecure, 
the zeroth order and the first order equation are automatically satisfied and 
the metered condition for the second order perturbation equation yields 


= 5 {x7 B 8} Yy Yo Up. (16) 
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N This i is exactly the same equation as that of Begbie and Born”) but the only 
Rx ‘symmetry of {x y, ß ö} is the interchangeability of « y with B 6. So eq. (16) 
x involves 45 constants and if one puts eq. (16) into the more familiar Christ offel 


form 


(17) 
2° a >= the wave equation contains 36 independent combinations of the {a y, ß 6} 
constants. 

_ Before discussing the relations between the atomistic calculations and the 
_ continuum theories, it is necessary to point out two features of the present 
calculations. The first is the direct applicability of the homogeneous deforma- 
tion method to infinite lattice. In particular the strain energy density is dedu- 
ced very easily, completely contraticting the proposition of Huang?) that 
this is impossible. It is difficult to visualise any continuum theory, which the 
lattice calculations must justify, that cannot have a strain energy function. 
The procedure used by the recent workers, viz., starting from the inifinite po- 
tential energy of the infinite lattice, is always plagued by the divergence 
difficulties. At least in this case, the force constant approach is superior. 

The second factor to be noticed is that the constants are 45 in number and 
only in the special case of central forces do they reduce to 21 independent con- 
stants. Huang?) has tried to invoke the initial stress conditions to reduce the 
general 45 constants to 21 constants, but as shown in Appendix C, the arguments 
is not applicable. In the general case {xy,ß 6} remain 45 independent constants, 
and the atomistic Christoffel equation eq. (17) involves 36 independent 
combinations of these constants. 


| $ 4. The correspondence between the atomistic and the continuum 
theories 


The general lattice calculations reveal two important facts, viz., tay + t,. 
and the Christoffel equation contains more than 21 independent constants. 
It is necessary to clarify how these facts can be understood on the continuum 
basis. 

The asymmetry of the stress tensor has been attributed to the presence 
of volume couples G, i. e., 

tea — tap =0 (18) 
by several workers (§ 5). But this cannot be so for two reasons. Firstly, as 
discussed in Appendix B, the lattice method of calculating the stresses as the 
sum of pairwise interactions can give only surface dependent terms and not 
volume dependent terms. So the asymmetry must be due to the presence of 
surface couples 6 (6 dS acting across any hypothetical surface of area dS) 


baa — bap 2 Jr. =% («By in lexical order) (19) 
(If, in a continuum, volume couples are present, as in polarised media, a term 


o Gg must be added to the left side of eq. (19). But this, as proved below, is 
not very significant for the symmetry of the elastic behaviour). Secondly, even 


22) K. Huang, Proc. Roy. Soc. London 208 A, 178 (1950). te 
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if volume couples are present, the elastic behaviour is adequately described 
by 21 constants and the Chrisstoffel equation cannot contain more than 21 
independent constants. This statement is easily verified from the expression 
for the strain energy. 

The strain energy density is obtained from the work done by the surface and 
volume terms 


dv =f ft, - du dS +f 6, -dwodS +[[feT- du dv + ff feG-dadv 
ou, 
öW = (=) — = bw, (tag tga) (20) 


where « ßy are in lexical order in the second term $,. re the svm- 
metric and antisymmetric parts 


Tap = (tap +tpa)/2; Tas = (tag — tea)/2 


and eqs. (18, 32), one has 5 


so that 


= = Tsp OE 


Using Hamilton’s principle, the equation of motion becomes 


and the bowidney 

>> En cos ß + n = x-component of force across unit area of normal n. (24) 


One sees that the fundamental equation (22—24) which adequately des- 
cribe the elastic behaviour of the body contain only Tt, and e,; and so the 21 
constants connecting these symmetric tensors are sufficient to describe the 
entire elastic behaviour of the body. Even in the presence of volume couples 
all the equations of the classical theory are valid if one replaces the t,; in the 


usual equations by T,, (i. e., slightly redefined shears) (Tiffen and Steven- 
2 


son?3), In particular, the wave propagation is described by e = = 


identical with the classical Christoffel equation unlike the atomistic case 
of eq. (16). 

Thus the differences between the exact atomistic theory and the classical 
theory of elasticity are due to the presence of interfacial couples, which are 
ignored in the latter. In the presence of the surface couples 


cw, 


(Note a slip in eq. (IV, 18) oir the second summation term involves 6,,.) 
For the plane wave propagation, the quantities 


Dap,y = Epy,« — Eay,B (26) 
2) R. Tiffen and A,C. Stevenson, Quart. J. Mech. Appl. Math. 9, 306 (1956). 
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do not vanish and so the surface couples interact with the rotational all 
leading to additional energy terms. Thus the increase in the number of con- 
stants in the wave equation is naturally to be expected. 

Before discussing the recent papers on the theory of elasticity it is well 
to emphasize that there is complete harmony between the atomistic and the 
continuum theories of elasticity. The inconsistencies noted by the earlier 
workers are due partly to the errors in the lattice calculations and partly to 
the neglect of the interfacial couples in the theory of infinitesimal elasticity. 


§ 5. Critique of the recent theories of crystal elasticity 


The foundations of the theory of elasticity have recently been examined by 
several workers. Most of the workers have not approached the subject abinitio 
from the continuum and the atomistic viewpoints, so that several subtle argu- 
ments have been overlooked by them. In particular, all their calculations are 
strictly valid only for infinitesimal deformations‘). 

The credit of conceiving the basic framework of the general lattice theory 
of elasticity with general interactions goes to Born and co-workers?)!?), 
However, there were some unsatisfactory features in their calculations which 
led to the subsequent re-investigations. Born and co-workers have overlooked 
that the rotation of a lattice cannot be represented by displacements alone and 
their use of the rotational invariance condition leads to relations which are 
not valid in the general theory. A second error in their work is the ad hoc 
introduction of the considerations of initial stresses in the natural stress problem 
in an attempt to introduce into the lattice calculations some of the symme- 
tries needed in the classical elasticity. 


Laval) and Viswanathan”) re-examined the theory of crystal elasti- 
city from the lattice point of view and showed for the first time that the wave 
equation should contain 45 constants. They used the homogeneous deformation 
method also. But both the authors have neither referred to the considerations 
of initial stresses advanced by Huang”) nor noticed the subtilities in the 
representation of the lattice rotations. Laval!*) has, however conjectured 
that the rotational invariance condition may be unsound, whereas Viswana- 
than has overlooked this point. 


Le Corre*), Raman and Viswanathan”) and Wooster and Joel?) 
have modified the continuum theory in order to reconcile it with the lattice 
theory and suggested the necessity of taking volume couples into considera- 
tion. Their use of t,3 as the 9 stress components and u,,, as the 9 strain com- 
ponents leads to the situation (eq. 21) that the strain energy can be altered 
by rigidly rotating the crystal, in obvious disagreement with the observed 
phenomena. As already pointed out in § 4, the presence of volume couples does 
not lead to any departure from the old classical behaviour. 

The present authors have given in this set of papers, a thorough discussion 
of the basic postulates of the theory of elasticity pointing out the inconsis- 


24) J. Laval, “L’etat solide” Congress Solvay, Stoops, Brussel (1951), p. 273. 

3) Y. Le Corre, ‘‘Elasticite et Piezoelectricite’’ Theses, Paris 1953; J. Physique 
Radium 19, 541 (1958). 

%) C. V. Raman and K. S. Viswanathan, Proc. Ind. Acad. Sci. 42 A, 1, 51 (1955). 
in W. A. Wooster and N. Joel, Nature 180, 430; 182, 1078 (1958). 
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tencies present in the earlier treatments. This covers the basic concepts of 
strains and stresses in the elastic continuum and in the elastic lattice. The 
subsequent development of the theory of elasticity rests on Hooke’s law, 
generalized by Cauchy, that the stress components are linear functions of 
the strain components. This is a separate independent axiom which in 
some cases (eg. nonlinear elasticity) may not be correct. (The statement that 
the lattice theory proves Hooke’s law is not quite correct. The first order 
approximation leads to Hooke’s law, but higher order terms can be incor- 
porated to give the third order elastic constants’). The extension of the Hoo- 
ke-Cauchy law is intimately connected with the equation of wave propaga- 
tion in the continuum mechanics and can be settled only by further rigorous 
studies. 

In order to test the validity of the theory, a detailed discussion of some 
clear cut experiments using such simple crystals like MgO, ZnS, Na(C1O,), 
Ba(NO,),, KDP, MgSO, etc. has been given recently by the authors”). 


Appendix A 
The representation of strain in a continuum and in a lattice 

The essential point’ brought out here are (7) if two states S and S’ differ 
by infinitesimal deformations there is no difference between the description 
of the strains with respect to S and that with respect to S’ and (ii) whereas the 
deformations of a continuum can be described by affine transformations 
(i.e., displacements of the points of the continuum), the deformations of a 
lattice must be represented by contact transformations (i.e., displacements 
and rotations of the lattice points). 

Strain is the geometrical description of the deformation of one state of 
the continuous body (called the deformed state S’) with respect to another 
state (initial state S). One assumes a one-to-one correspondence between the 
points of S’ and S so that in some fixed axes of reference, the position vectors 
of the points in S’ and S are uniquely connected by r’ = r + u. It is also sup- 
posed that uw is continuous and has partial derivatives of as high an order as 
required. 

The length ds’ of a line element in S’ and its length ds in S are compared. 


1 
One has dx, = dr, day so that (ds’? — ds?) can be 
B B 


= dt, dx; where 
1 (ou, dus, dx, Oug 


=F dz, du, 22). 
This ¢,g is the Lagrangean (material) strain tensor. If one had set 

= — ds?) = nap das 


written as 


al 


where 
one has the Eulerian (spatial) strain. p 2 4 
%) F. Birch, Physic. Rev. 71, 809 (1947). 


2) R.S. Krishnan, V.Chandrasekharan and E.S. Rajagopal, Ind. Inst. Sci. 
Golden Jubilee Res. Vol. p. 150 (1959). 
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" ee The material and the spatial descriptions coincide if wert! 


so that their powers > 1 may be neglected. Then one has the infinitesimal 
strains and rotations 


1 1 
Exp + Wap = — Wa,p) - (32) 


The general condition for the infinitesimal deformations is the smallness of the 
displacement gradients. It is not sufficient to take the strains as small. As an 
example, consider a deformation where the strain is smail but the associated 
rotation is large. In that case u,,, are not small, and the descriptions using 
the initial co-ordinates (eq. 29) and the final co-ordinates (eq. 30) differ ap- 
reciably. Thus if the displacement gradients are infinitesimal there 
is no difference between the description of the strains in the ini- 
tial and deformed co-ordinates”), 

Since strain is a geometrical description, the strain components obey the 
geometric condition that the space is locally Euclidean and the vanishing 
of the Riemann-Cristoffel tensor leads to the usual Saint Venant com- 
patibility conditions. The significance of these conditions have been questioned 
sometimes*!)32) but time and again they have been convincingly refuted *3)*), 

Whereas in a continuous body the deformations are specified by the dis- 
placements u or by affine transformations, the extension to the system of 
discrete particles requires some consideration. One passes from the continuum 
to the lattice by marking of the appropriate volume elements and contracting 


# 


c) a) 
Ben Fig. 1. Deformations and rotations in a continuum and in a lattice a 


%) I. S. Sokolnikoff, "Mathematical Theory of Elastieity‘ MaGraw Hill, New York 
(1956), p. 30, 64. See also Z. Karni and M.Reiner, Bull. Res. Council, Israel Se, 


89 (1960). 
31) C. Eckart, Physic. Rev. 78, 373 (1948). 1- on 

3) E.H. Brown, J. Res. Nat. Bur. Standards 59, 421 (1957). mit: a :" 
8) C. Truesdell, J. Rational Mech. Math. 1, 145 (1952). a RR 
3%) J.L. Ericksen, Math. Revs. 20, No. 557 (1959). aw 
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each volume element to the corresponding point of the composite lattice. To 
visualise an anisotropic continuum, one attributes different properties along 
the different directions emanating from each particle (Hellinger®)). A part 
of a body is represented in Fig. 1a and its equivalent lattice representation 
in Fig. 1b. The result of some arbitrary affine transformation on the body 
la be 1c which will be represented in the lattice form in 1d. 

One sees that in the case of the lattice each lattice point suffers the corre- 
sponding displacements and in addition the corresponding rotation. Thus one 
cannot represent the operations by mere displacements alone. One needs an 
infinitesimal contact transformation (Corbon and Stehle**)) to represent 
the rotations in a noncentral force field. This is being realized in recent years?®®). 

This point has an important bearing on some of the recent lattice dynamic 
calculations (§ 2). 


The notion of stresses in a continuum and in a lattice 


The essential points brought out here are (i) if two states S and 8’ differ by 
infinitesimal deformations, there is no difference in the descriptions of stress 
with respect to S and with respect to S’ and (ii) in an atomistic lattice, the con- 
siderations of pair-wise interactions leads only to surface dependent reactions 
and not to the volume terms. 

The Cauchy stress principle is the basis for the development of the ideas 
of stress. The reactions in an elastic body consist of a surface part represented 
by a force tdS on the area dS and a volume part o Tdv on the volume element 
dv. The direction of t is in general different from that of the normal n of dS and 
across the interface, equal and opposite forces balance each other. Then one has 


Strictly the equilibrium must be considered in the strained state S’ and eq. 
(23) must contain x’ instead of x. But 


so that, in view of eq. (30) and the continuity of stress 0t,/Ozx, ~ Ot. /dx,,. 
Thus if the displacement gradients are infinitesimal, there is no 
diffenrence between the descriptions of the stress in the initial and 
the deformed states (Sokolnikoff*®)). 

In the text, the presence of volume couples and surface couples is consi- 
dered. The existence of the interfacial couples constitutes an extension of 
Cauchy’s principle. 


%) E. Hellinger, Encyk. Math. Wiss. IV—4, 609 (1914). 

%) H.C.Corbon, and P. Stehle, "Classical Mechanics“, Wiley, New York (1950), 
p. 324. 

36a) N.Oshima in ”Unifying study of the basic problems of engineering“ Ed. 
K. Kondo, Tokyo Vol I, p 563 (1959); C. Truesdell, Appl. Mech. Revs. 12, 77 
(1959); E. L. Aero and E. V. Ruvshinskii, Soviet Physics—Solid State 2, 1272 
1960). 
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In the case of a lattice, the concept of stress must be introduced as the 
mean value over an area containing a sufficiently large number of unit 
cells but still sufficiently small on a bulk scale. The lines of action of the forces 
between the particles on either side of the area will pass through this area. The 
force resultant of all such lines of forces acting from one part of the lattice 
on the other is the resultant traction over the surface and per unit area this 
is taken as the stress (Love?”?)). In the lattice calculations it is convenient to 
perform the calculations on infinite lattices to avoid the conditions of the 
boundary atoms. One takes a hypothetical interface dividing the crystal into 
two infinite halves and computes the stresses across this interface. Since the 
lines of interaction of any atom of one half with any of the other half passes 
across this interface, all possible interactions between the two halves are in- 
cluded. It is very important to realise that all the interactions give in 
this mode of calculation only the surface dependent terms. 

In the lattice calculations, the volume dependent terms can arise due to the 
interaction with external (gravitational or electromagnetic) fields which are 
such that for some interfaces one cannot take the forces as being between the 
two halves. A possibly different way of getting the volume terms is to calculate 
in a dipolar lattice, the potential of one half of the lattice when one can obtain 
the forces acting on the ions in a given volume. But then one has replaced the 
lattice by a smoothed continuum and one is not considering pair-wise inter- 
actions. As soon as one considers the interactions pair-wise (in the case of two 
body forces), one lands only in the surface dependent terms. ; 4 


> Appendix € 
apes The states of the elastic body 


So far the discussion has been on the relationship between the initial 
state Sand the deformed state S’ which differ through an infinitesimal defor- 
mation. Frequently one considers a hypothetical stress-free state, the natural 
state S°. Since one is interested, in any case, in the final state 8’ two possible 
cases, the natural stress problem connecting S’ and S and the initial 
stress problem connecting S’ and S° are of interest. The initial stress prob- 
lem occurs when one superposes strains on a body having ’’frozen“ stresses as 
in many technical and geophysical applications. The natural stress problem 
is due to the very concept of the solid state. The idea of cohesion presupposes 
that the two halves of the body are held together from flying apart. In the 
absence of the body and the d’Alembertian forces, the equilibrium condi- 
tion is Nz0t,5/0x; = 0. If the body is finite, the condition on the traction-free 
surface Sat.gß-n = 0 leads to t,; = 0. But if the body is infinite, the boun- 
dary conditions do not exist and t,; = constant, is possible. This is the state 
S of the body and all its properties are supposed to be known and one is Te- 
quired to connect S and 8”. 

One can show that while the initial stress problem is sophisticated, the 
natural stress problem is almost trivial. First, to clarify the difference between 
the two cases, consider a one dimensional strain, an element having the lengths 
1,1+.¢, and 1+¢,+ 6; (a <1-+e,) in the states S°, S and 8’ respecti- 


ma. E. H. Love, “Mathematical Theory of Elasticity”, Dover, New York (1944), 
p. 617. 
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vely. In the initial strain case, the strain in the state S is [(1 + &,)? — 1?]/2.1 = 
& + 7/2 while in S’ the strain is [(1 +, +? — 17)/2.1 e +e, + 
&/2 so that the additional strain ¢, + ¢,¢, is coupled to the original strain 
& leading to complications. In the natural stress case, one needs the strain of 
S’ with respect to S which is [(1 +2, + €,)? — (1 + &,)?]/2- (1 +4) = e and 
so no coupling with e, occurs. Again, the strain energy of S’ will be of the form 


W=4e +z08 so that the derived stress is t = OW/de = A + Ce. Here, 


one is interested in the case of uniform initial stress A = const. Then in the 
equations of motion the constant term vanishes on finding öt/0x and the beha- 
viour is indistinguishable from the A = 0 case. N, doubt the boundary condi- 
tions will involve the constant term A, but in the infinite lattice, the boundary 


conditions do not exist. rh 

The state S is characterized by the equations ve 


and in fact, the uniform initial stresses f,(x) (as the functions of the Langran- 
gean co-ordinates x,) satisfy 


Now the state S’ is defined by Sob pat ee 
> = = (36a) 
3 


en t'(x’) refers to the total stress in the Eulerian c-ordinates x;. This 
ti;(a’) is the resultant of the applied stress r(x’) and the initial stress t°(x’) 
now referred w the x; co-ordinates. So 


From (35a) (seo also Green?!*), one 


(878) 


so that 
;(2") = T(x’). (37b) 


The homogeneity of m deformation yields ee 
at? = 0; = 0 (38a) 
and in particular 
The strain energy density is 
= 


The applied strain is linearly connected to the applied stress and does not 
change the magnitude of the original stress so that 


ty= Exı- 
Then, from the second derivative of W, one gets 
Cijkl = Cklij- 
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The relations (37b), (38b), (39) show that the ratio of the applied stress to 
the applied strain has all the symmetries of the ordinary constants. If only 
surface and volume forces are present there will be only 21 constants. Of 
course, the symmetry of the state S will operate on the elastic response, but 
the magnitudes of the initial frozen stresses do not alter the symmetry of the 
stress and strain tensors. As examples, one may cite the recent work of Bha- 
gavantam and Chelam®®) who have shown that a cubic crystal with an arbi- 
trary triaxial stress or shear (symmetry of the initial state S is monoclinic) 
will need 13 constants only. 

This is a conclusion clearly stated only in a few places*), As Truesdell 
has remarked the extensive literature on the initial stress problem is not al- 
ways clear or precise. Most of the workers have discussed the problem of small 
strains superposed on initial strains and have taken the unstrained state as the 
reference, mentioning the natural strain problem as being trivial. As mentioned 
earlier, most of the lattice dynamic calculations are valid only for infinitesimal 
deformations. For example, the wave equation eq. (16) is valid only for infini- 
tesimal deformations and in the case of finite amplitude waves additional 
terms occur even if the initial state is stress free#0)41), 
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Uber die Struktur von Caesium-Wismut-Fotokathoden 


Die — der Bahia des Systems Cs—Bi wird an dünnen 
Schichten durch Elektronenbeugung untersucht. Auswertbare Aufnahmen 
werden nur für die Verbindung Cs,Bi nach längerem Tempern der Schichten 
erhalten. Als Ergebnis wird in Übereinstimmung mit aus Röntgenstruktur- 
untersuchungen bekannten Angaben gefunden, daß Cs,Bi auch in dünnen 
Schichten nach dem Cu,Al(DO,)-Typ kristallisiert, mit a = 9,292 + 0,013 Ä. 


Für die die Struktur von Caesium-Wismut-Verbindungen ergibt sich aus Rönt- 
genstrahlinterferenzanalysen dieser intermetallischen Verbindungen an kom- 
pakten Materialien nach Klemm und Mitarbeitern!), daß CsBi, in einer 
Laves-Phase vom MgCu,-Typ (C 15-Typ, O,—Fd3m, M = 8, a = 9,741 A) 
und Cs,Bi kubisch nach dem Cu, Al- Typ (DO,-Typ, 0h—-Fm3m, M = 4, 
a = 9,291 A) kristallisiert. 

Die Atome dieser hexakisoktaedrisch kristallisierenden A!BY-Verbindung 
sind in folgender Weise auf die Gitterplätze verteilt. An den Stellen (0,0,0) 
(0,1/2,1/2) (1/2,0,1/2) (1/2,1/2,0) befinden sich die Bi-Atome und an den Gitter- 
stellen ((0,0,0) (0,1/2,1/2) (1/2,0,1/2) (1/2,1/2,0)) + (1/4,1/4,1/4) (1/2,1/2,1/2) 
(3/4,3/4,3/4) die Cs-Atome. Eine Veröffentlichung über eine Bestimmung dieser 
Struktur an dünnen Schichten ist uns nicht bekannt. 

In einer Arbeit über das Verhalten der Alkalimetalle zu den Elementen 
der 5. Hauptgruppe des periodischen Systems gibt Gnutzmann?) auf Grund 
thermischer Analysen das binäre Zustandsdiagramm für das System Caesium- 
Wismut an, dieses stimmt recht gut mit dem von Zhuravlyov?°) aufgestellten 
überein. Der Widerstand dünner Cs—Bi-Schichten als Funktion des Verhält- 
nisses der beteiligten Metalle zeigt Maxima, die sehr wahrscheinlich den Stellen 
stöchiometrischer Zusammensetzung entsprechen und die in Analogie zu dem 
binären Zustandsdiagramm für CsBi, und Cs,Bi besonders ausgeprägt sind. 

Für die technische Anwendung der Fotokathoden kommt nur das stöchio- 
metrische Verhältnis A!: BY = 3:14)5) in Frage. Gelegentlich wird sogar 


1) W. Klemm, Angew. Chemie 66, 148 (1954). N 
2) G.Gnutzmann, Diss. Münster 1953. 
3) N.N. Zhuravlyov, Z. exper. teor. Fiz. 34, 827 (1958). a he ae 
4) A.Sommer, Proc. Roy. Soc. London 55, 145 (1943). rey Me = 

5) H. Miyazawa, J. Phys. Soc. Japan 8, 169 (1953). at i 
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noch zusätzlich durch Hinzufügen von Sauerstoff die Empfindlichkeit der 


Fotokathoden erhöht, jedoch soll darauf in diesem Zusammenhang nicht ein- = 
gegangen werden. barı 

Nach einer Arbeit von Sommer®) liegt das Maximum des äußeren Foto- ver] 
effektes für AIBY nicht genau beim stöchiometrischen Verhältnis, sondern mu) 
bei einem Zustand thermischen Gleichgewichtes, der bei einem bestimmten Sta 


Caesium-Unterschuß liegt. Im Gegensatz dazu finden Jack u. Wachtel”), 
daß bei Cs,Sb der Homogenitätsbereich bei Cs3 + „Sbı — „ liegt, mit x x + 0,05. 900 

Unter der Voraussetzung, daB im Formierdiagramm das Minimum der 
Leitfähigkeit einer Zusammensetzung im stöchiometrischen Verhältnis ent- Pla 
spricht, sollen die Strukturen von dünnen Cs,Bi-Schichten untersucht werden. 


Die Versuchsanordnung 


Die Atomformamplituden für die Röntgenbeugung liegen in der Größen- buy 
A ordnung von 10-1, die für Elektronenbeugung bei 10-*. In die Intensität der bon 
*F Reflexe geht das Quadrat der Atomformamplitude ein, d. h., die Intensität pn 
, der Elektronenbeugung ist ungefahr um den Faktor 10% größer als die der 2 
af Röntgenbeugung. Wegen dieser geringen Wechselwirkung sind für die Rönt- e 
_: genbeugung an dünnen Schichten Schichtdicken von 10%... 10° A notwendig. a 
® Die untere Grenze dieser Schichtdickenangabe ist gleichzeitig die obere Grenze . , 
. der Schichtdicke einer Fotokathode, von der noch ein merklicher Fotoeffekt r 

> zu erwarten ist. Im Gegensatz zur Röntgenstrukturanalyse eignen sich Elek- es 

N tronenstrahlen zur Strukturbestimmung der vorliegenden dünnen Schichten = 

I besser. Es genügen schon Schichtdicken von 10... 100 A für eine merkliche hen 
Streuung der Elektronen. Me 
Die von der Industrie hergestellten Elektronoskope sowie die Elektronen- int 


mikroskope, mit denen auch Beugungsaufnahmen gemacht werden können, - 
sind für unsere Untersuchungen nicht verwendbar, da, bedingt durch die große 


Anzahl von Dichtungen und Schliffen und durch die Fotoplatten, höchstens po 
ein Druck von 5 - 10-5 Torr erreicht werden kann. Dieses Vakuum reicht zur} Be 
sauberen Herstellung der Caesium-Wismut-Fotokathoden nicht aus. Es wird a 
daher eine Apparatur verwendet, die vollständig aus Glas besteht und im Hoch- 90 
vakuumteil bis zu den Diffusionspumpen keine Schliffe oder Dichtungen be- Er 
sitzt. 

Das verwendete Versuchsrohr hat einen Durchmesser von 80 mm. Die Auf- > 
nahmen werden in der Registrierebene auf einem sehr feinkörnigen Leucht- He 


schirm (75°/, der Körner < 5,4) erzeugt und von außen fotografiert. Die At 
Glaswand zwischen Anode und Leuchtschirm und der Leuchtschirm selbst bi 
sind zur Vermeidung von Aufladungserscheinungen mit einer Aluminium- 


schicht bedeckt. Auf der AuBenseite des Rohres wird, um eine hohe Isolation = 
zu erreichen, zwischen Kathoden- und Anodenzuleitung eine Silikonlackschicht P 
eingebrannt. U 
. Zur Strahlerzeugung wird eine Fernfokuskathode nach Steigerwald Ca 
. verwendet. Die Anodenblende aus Molybdänblech hat einen Durchmesser von nit 
* 50 u. Die Objektträgerhalterung und damit auch der Objekttrager, falls er 
i} sich im Strahlengang befindet, sind durch eine Wolframwendel bis 600 °C 
i ‘) A. H. Sommer, J. Appl. Phys. 29, 1568 (1958). 
7) K.H. Jack u. M.M. Wachtel, Proc. Roy. Soc. London A 289, 46 (1957). W 
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t der 
ein- 


heizbar. Die Objektträgertemperatur wird mit einem Thermoelement gemessen. 
Für die Beschleunigungsspannung wird eine bis 100 kV kontinuierlich regel- 
bare Gleichspannung mit einer Welligkeit von 6 - 10% verwendet. Der Strahl- 


Foto- verlauf wird magnetisch durch eine Justierspule korrigiert. Für die Bestim- 
dern mung der Aufnahmekonstanten L-A werden NaCl- und Bi-Schichten als 
mten Standardpräparate verwendet. 
el"), Der Objektträger besteht aus zwei Teilen, der Trägerfolie, ein ungefähr 
0,05. 200 Ä dickes Kohlenstoffhäutchen, das auf einem feinmaschigen (6000 Maschen 
| der pro cm?) Silberdrahtnetz liegt und einer Meßplatte aus Glas, auf der zwei 
bai Platinkontakte zur Widerstands- und Fotoeffektmessung aufgebracht sind. 
rden. 
Versuchsdurchführung 
Um eine eindeutige Interpretation der Versuchsergebnisse zu ermöglichen, 
Bel wird auf möglichst reine Ausgangsmaterialien großer Wert gelegt. Die spek- 
+ aie tralanalytischen Untersuchungen ergeben, daß in dem verwendeten Wismut 
sität noch Spuren von Silber und Kalzium und im Caesium im wesentlichen noch 
 dü Spuren von Kalzium, Kupfer, Magnesium und Silizium vorhanden sind. Zur 
ont Herstellung der Schichten wird bei einem Druck von < 10-* Torr mit Hilfe 
ndig eines Glühsenders das Wismut auf den Objektträger aufgedampft. Um eine 
a für das Verhältnis von Reflexintensität zu Untergrundintensität optimale 
fekt Schichtdicke des Cs,Bi von ungefähr 500 A zu erreichen, muß der beim For- 
‘ek. mierprozeB auftretende Quellfaktor bei der Erzeugung der Bi-Schicht berück- 
— sichtigt werden. Aus dem reduzierten Molvolumen von Cs,Bi und dem Atom- 
liche volumen des Bi ergibt sich ein Quellfaktor beim Übergang von der kristallinen 
Wismutschicht zur ausformierten Cs,Bi-Schicht von 5,68. Danach muß die 
ui Schichtdicke des Wismuts ungefähr 85 Ä betragen. Die Schichtdicken werden 
m interferometrisch bestimmt und der Widerstand der Bi-Schicht in Abhängig- 
nail keit von der Dicke gemessen. Um eine bestimmte Schichtdicke zu erhalten, 
heal oN zum Erreichen des entsprechenden Schichtwiderstandes Wismut auf- 
| gedampft. 
al Zur Einleitung des Formierprozesses wird ein Triimmerventil zu einer der 
Bi; angesetzten Caesiumampullen durchgeschlagen und anschließend bei 150 bis 
wu 200 °C das Caesium in die Schicht einformiert. Für einige Proben wird beim 
| Erreichen des ersten Formierzustandes der Formierprozeß unterbrochen und 
Aut. die Probe in den Elektronenstrahl gebracht. Jedoch sind keine indizierbaren 
cht- Debye-Scherrer-Ringe von CsBi, festzustellen. Auch der Versuch, durch 
Die Heizen der Probe ein Wachstum der Kristalle zu begiinstigen, ergibt keine 
Ibst Aufnahme mit indizierbaren Ringen. Bei Durchführung des Formierprozesses 
we, bis zum Erreichen des zweiten Widerstandsmaximums ergeben sich zunächst 
vu auch für die Cs,Bi-Schichten keine indizierbaren Aufnahmen. Nach 20stün- 
‘cht digem Heizen bei einer Temperatur von 210 °C wächst bei einigen Proben die 
Größe der Kristallite, und die Debye-Scherrer-Ringe heben sich besser vom 
ald Untergrund ab. Nach dem 20stiindigen Tempern der Probe wird noch einmal 
res Caesium einformiert, wobei der Widerstand der Schicht zunächst wenig zu- 
~~ nimmt und dann wieder abfällt. 
°C MeBergebnisse 
Abb. 1 zeigt die Abhängigkeit des Widerstandes des Systems Caesium- 
Wismut während des Formierprozesses von der Formierdauer. Von einem 
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_ Anfangswiderstand des metallischen Wismuts von 2700 Ohm — das entspricht 


für unsere ungefähr der gewünschten Schichtdicke der ausformier. 
> = o ten Cs,Bi-Schicht von 550A — steigt ents 
der Widerstand für den ersten Formier. 
| of zustand CsBi, bis 5 - 10 Ohm an, 


darauf bis 2 - 106 Ohm abzufallen. An. 


| | schlieBend zeigt die Formierkurve ein 4 
kleineres Maximum von 3,7 - 106 Ohm | 2¥#! 
für den zweiten Formierzustand Cs,Bi. erha 
Diese Form der. Formierkurven ist 


für das System Cs-Bi charakteristisch. 4 
Fiir verschiedene Proben unterscheiden sir 

sich die Höhen der Maxima bei gleichem | “8! 
Anfangswiderstand his zu einem Fak. f tan 
tor zwei. Die Kurve II in Abb. 1 zeigt, +0 


„ daß erstens der Widerstand der Schicht | ™& 
6 während des Temperns um etwa zwei f &°™ 
Tm Zehnerpotenzen zugenommen hat und 


ke 


Br. zweitens zu Beginn des nochmaligen Ref! 
0 6 3045 60 75 BE: 05 120 Zuführens von Caesium der Widerstand | Sch 
der Schicht wenig zunimmt, um dann 
a Abb.1. Formierkurve des Systems Cs-Bi abzufallen. Die Widerstandszunahme star 
während des Temperns ist wahrschein- f 1 

lich durch teilweises Abdampfen der Schicht zu erklären. Dabei ist die f elek 
Kurve II in das Formierdiagramm so eingetragen, daß die Lage ihres Maxi- f Pt 

mums mit der Lage des zweiten Maximums der Kurve I übereinstimmt. 
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Abb. 2 . AuBerer Fotoeffekt und Widerstand einer diinnen Stu 

= Cs, „Bi,_,-Schicht bei Cs-Abgabe wir 


Abb. 2 zeigt die Abhängigkeit des Widerstandes und des äußeren Foto- Un 
Ss effekts einer Css _ ,B,.,-Schicht von der Zeit, während der die Schicht einer 
Temperatur von 200...210 °C ausgesetzt war. Charakteristisch dafür ist, 
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daß für einen Zustand mit Caesiumunterschuß der Fotoeffekt sein Maximum 
erreicht, während für höheren, wahrscheinlich dem stöchiometrischen Verhältnis 
entsprechenden Caesiumgehalt die Schicht ein Widerstandsmaximum aufweist. 
Bei einigen der untersuchten Schichten wurde nach längerem Tempern der 
äußere Fotoeffekt nahezu gleich Null. 


Auswertbare Beugungsaufnahmen werden erstens nach dem Tempern und 
zweitens nach der Caesiumzugabe an der Stelle des Widerstandsmaximums 


erhalten. In Tabelle 1 sind in 
Spalte 1 die Indizes, in Spalte 2 
die gemessenen Werte für 
2sind/2 angegeben. Daraus 
ergibt sich für die Gitterkon- 
stante ein Wert von a = 9,292 
+ 0,013 Ä. Die Werte in Klam- 
mern in Spalte 2 sind aus dem 
gemessenen a berechnet. 

Um die Intensitäten der 
Reflexe aus den gemessenen 
Schwärzungskurven berechnen 
zu können, müssen zwei Um- 
stände beachtet werden : Erstens 
ist bei der Umwandlung des 
elektronischen Signals in ein 
optisches die Abhängigkeit der 
Lichtausbeute beim Fluores- 
zenzmaterial von der Elek- 
tronenstrahlstromdichte und 
zweitens die Abhängigkeit der 
Schwärzung der fotoempfind- 
lichen Schicht von der Beleuch- 
tungsstärke zu berücksichtigen. 
Für die Anregung von Fluo- 
reszenzmaterial durch Elek- 
tronenbestrahlung gilt für den 
Fall nicht zu hoher Strom- 
dichte, d.h. bis kurz vor Be- 
ginn des Sättigungsgebietes, 
daß die Fluoreszenzmaterialien 
ein lineares Verhalten zeigen ®)®). 
Zur Bestimmung der Abhängig- 
keit der Schwärzung von der 


Tabelle 1 
Wertefür 2sin 2 undfür die beobachteten 
Intensitäten 
Beschleunigungsspannung U = 50 kV 
hkl 2sind | J, 

111 0,01003 | 1800 1600 
200 001152 700 | 420 
290 0,01628 | 17800 | 16500 
311 0,01916 | 30 
222 (0,01993) _ _— 
400 0,02310 1320 1320 
331 0,02517 15 
420 0,02575 13 we 
422 0,02830 721 721 
511,333 0,03001 5 4 
440 0,03257 93 98 
531 0,03408 4 3 
600,442 0,03457 2 ? 
620 0,03642 91 
533 (0,03774) _ 
622 (0,03820) _ 
444 0,04000 15 
711,551 (0,04114) _ 
640 | (0,04154 = 
642 0,04276 43 
553 (0,04431) _ 
800 | -0,04602 3 
733 | (0,04719) _ 
820,644 | (0,04753) _ 
822,660 0,04835 13 
751,555 (0,04989) _ 
622 (0,05024) _ 
840 0,05144 4 


Die beobachteten Intensitäten sind für den Reflex 
422 einander gleichgesetzt. 


Intensität werden in n-Stufen aufgeteilte Eichplatten verwendet, wobei die v-te 
Stufe einer Platte mit f - 2” s belichtet wird. Für die Schwärzungsmessungen 
wird ein Schnellphotometer G II von Zeiss, Jena, verwendet. Um die Inten- 
sität einer Linie eines Debye-Scherrer-Ringes unter Berücksichtigung des 
Untergrundes zu erhalten, muß die Intensität des Untergrundes von der 


8) W. Schnabel, Arch. f. Elektrotechn. XXVIII, 789 (1934). 


°) F.Eckart, E lektronenoptise he Bildwandler und las re Leipzig 1956. 
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Gesamtintensität (Linie + Untergrund) abgezogen werden. Die Intensitäts. unt 
differenzen ergeben sich über die Schwärzungs-Intensitätskurve aus den sch 


Schwärzungsdifferenzen bei geeigneter Mittelwertsbildung für die Untergrund. die 
schwärzung. Es wird dabei nach einem vereinfachten Rechenverfahren von kri 
Honerjäger-Sohm und Kaiser”) gearbeitet. Die sich danach ergebenden dal 
Intensitäten J; (nach dem Tempern der Schichten) und J; (nach nochmaliger n-l 
Caesium-Zufuhr an der Stelle des Widerstandsmaximums) sind in Tabelle | da; 
in Spalte 3 und 4 eingetragen. 


Diskussion 


Entsprechend den Strukturangaben für Cs,Bi von Klemm!) ergibt sich Ki 
für die Strukturamplitude wegen 


= Xn fen 

d.h. F, ner = (fei + 4 4 

4 gritty tts) 4 +hy)), 

1. Für h, k,l gemischt (gerade und ungerade) ist F, ,,, = 0 

2. Fir h+k+1=4n ist F, = 4 fai + 12 fos 

3. Fir h+k+l=4n+4+2 und 
für h, k,l alle ungerade ist F, = 4(fsi — fcs)- 

Ein Vergleich mit Tabelle 1 zeigt, daB die auftretenden Reflexe den Be- 


dingungen 1, 2 und 3 genügen. Die Intensität der nicht gefundenen Reflexe, 
die Punkt 3 entsprechen, 


Tabelle 2 
Vergleich des Leitungstypes mit der Kristalli- der 
sationsart der I-V-Verbindungen ndlichkeitsgrenze er 
Meßanordnung. 
BY Leitungstyp | Sristalli- | Operschuß Die dünnen CsBi 
Schichten kristallisieren 
Na,Sb n | hexagonal 
K,Sb n hexagonal ER stimmung mit dem aus 
Rb,Sb n hexagonal A Röntgenstrukturunter- 
: suchungen bekannten Er- 
a,Bi exagona 
K,Bi ? hexagonal ? se dem Cu,Al 
3 ad 

Cs,Sb p kubisch B Abschließend sei noch 
Cs,Bi_ ? kubisch B auf folgenden, uns be- 
Rb,Bi ! kubisch ? . merkenswert erscheinen- 
Li,Bi ? kubisch d Z 

ß-Li,Sb ? kubisch Baden en Zusammen g 
Na,KSb p kubisch ee gewiesen. In Tabelle ? 
03,KSb p ? sind von einigen A!BY. 


Literatur zur Kristallisationsart: 1, 2, 7, 11, 12. Verbindungen Kristalli- 
_ Literatur zu Leitungstyp und Überschußmetall: 6, 13. sationsart, Leitungstyp 


10) M. Honerjäger-Sohm u. H.Kaiser, Spectrochim. Acta 2, 396 (1944). 
4) J.J. Scheer u. P. Zalm, Philips Res. Rep. 14, 143 (1959). 
12) W.Dorn, Diss. Münster 1955. 

13) W,E. Spicer, Physic. Rev. 112, 114 (1958). 
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und Angabe, welches Metall im thermischen Gleichgewichtszustand im Uber- 
schuß auftritt, zusammengestellt. Nach dieser Tabelle erscheint es so, daß 
die hexagonal kristallisierenden I—V-Verbindungen n-leitend, die kubisch 
kristallisierenden p-leitend sind. Verstärkt wird diese Auffassung dadurch, 
daß die beiden hexagonal kristallisierenden Verbindungen K,Sb und Na,Sb 
n-leitend sind und A-Überschuß aufweisen, die Multialkaliverbindung KNa,Sb 
dagegen, die nach Scheer und Zalm!!) kubisch kristallisiert, p-leitend ist. 


Herrn Prof. Dr. OÖ. Hachenberg danke ich recht herzlich für die Anregung 
zu der vorliegenden Arbeit und für sein förderndes Interesse. Herrn Dr. C. 
Kunze danke ich sehr herzlich für zahlreiche fördernde Diskussionen. 


Berlin-Adlershof, Heinrich-Hertz-Institut für Schwingungsforschung 
der Deutschen Akademie der Wissenschaften zu Berlin. 
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Berichtigung und Erganzung der Arbeiten 
„Die Richtungsabhangigkeit der Kristallenergie“') 
und „Die Richtungsabhängigkeit der Magnetostriktion?)“ 


Von W. Döring und G. Simon Re 


In den genannten Arbeiten wurde die allgemeine Gestalt einer skalaren 
und tensoriellen Funktion einer Richtung berechnet, welche mit der Symmetrie 
einer vorgegebenen Kristallklasse verträglich ist. Darin wurden die Richtungs- 
kosinus der Richtung wie die Komponenten eines Vektors behandelt. Bei der 
Anwendung der Ergebnisse auf den Ferromagnetismus ist aber zu beachten, 
daß die Magnetisierung ein Pseudovektor ist. Bei einer Spiegelung ändern 
sich bei Pseudovektoren die Komponenten senkrecht zur Spiegelebene nicht, 
während die Komponenten parallel zur Spiegelebene ihr Vorzeichen ändern. 
Ein Pseudovektor verhält sich also bei einer Spiegelung genauso wie ein Vektor 
bei einer Drehung um 180° um eine Achse senkrecht zur Spiegelebene. Eine 
Funktion der Komponenten eines Pseudovektors, welche gegenüber den 
Symmetrieoperationen einer Kristallklasse X invariant ist, hat daher die gleiche 
Gestalt wie eine Funktion der Komponenten eines Vektors, welche invariant 
ist gegenüber den Symmetrieoperationen einer anderen Kristallklasse, die 
aus K durch Ersetzen aller Spiegelebenen durch zweizählige Drehachsen 
senkrecht dazu und Streichen der Inversion hervorgeht. Daraus ergibt sich 
z. B., daß für die gegenüber O, und 7, invarianten Funktionen eines Pseu- 
dovektors dasselbe gilt wie für die gegenüber O invarianten Funktionen eines 
Vektors. 


Da ein Tensor 2. Stufe sich gegenüber einer Spiegelung und einer 2zähligen 
Drehung senkrecht zur Spiegelebene gleich verhält, gilt für einen invarianten 
Tensor das gleiche wie für einen Skalar. 


Beim Magnetismus kommt allerdings noch eine weitere Symmetrie hinzu. 
Eine Umkehr aller Komponenten der Magnetisierung entspricht einer Umkehr 
aller Geschwindigkeiten. Eine solche läßt bei fester Lage der Atomkerne und 
Abwesenheit eines äußeren Magnetfeldes die Energie des energetisch tiefsten 
Zustandes invariant, ändert also die Kristallenergie und den Deformationstensor 
der Magnetostriktion nicht. Ein Gitter verhält sich also hinsichtlich dieser 
magnetischen Eigenschaften stets so, als ob außer der Symmetrie, welche in 
der Anordnung der Atomkerne zum Ausdruck kommt, auch noch die Inversion 
Symmetrieelement wäre. 


Daraus ergibt sich die folgende Tabelle. Links steht die Kristallklasse, 
die röntgenographisch bestimmt werden kann und der Anoränung der Kerne 


Bu) a W.Döring, Ann. Physik (7) 1, 102 (1958). 
} ane W. Döring, u. G. Simon, Ann. Physik (7) 5, 373 (1950). 
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entspricht. Rechts steht die Kristallklasse, unter welcher die Form der Funk- 
tion für Kristallenergie und Magnetostriktion in den genannten Arbeiten zu 
finden ist. 


eines Vektors die gleiche Gestalt haben 
wie die Funktionen des Pseudovektors der 
Magnetisierung in der links angegebenen 


| Kristallklasse, bei welcher die Funktionen 


Kristallklasse der Atomkernanordnung 


lasse 
Cys Sys Cun 
Covi Poa D,; Dy Du = 
Zürich (Schweiz), IBM Forschungslaboratorium Adliswil und re, 
Californien (USA). 
Bei der Redaktion eingegangen am 31. Oktober 1960. A 
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Von Witlof Brunner und Harry 


Inhaltsübersicht 


Die von den Verfassern in früheren Arbeiten nacheinander entwickelten 
Vorstellungen zur Theorieder Kernspaltung werden systematisch dargestellt. Aus 
der Grundkonzeption einer schalenstrukturabhängigen Deformation der Bruch- 
stücke ,,im Augenblick“ ihrer Trennung folgen magische Effekte für das Kern- 
kraftpotential zwischen den beiden Fragmenten und die mittlere kinetische 
Energie der Bruchstücke, deren Zusammenspiel die beobachtete Spaltungs- 
asymmetrie quantitativ verständlich macht. Zur Illustration wird die spontane 
Spaltung von Cf 2% diskutiert. Weiterhin wird die Tatsache der symmetrischen 
Spaltung von Bi2® (durch Deuteronen) verständlich gemacht. Schließlich 

| wird eine qualitative Erklärung für die bekannte Abschwächung der Spaltungs- 


asymmetrie bei Erhöhung der Anregungsenergie des sich spaltenden Kerns 
auf der Grundlage der entwickelten allgemeinen Vorstellungen gegeben. 


In der vorliegenden Arbeit sollen zunächst die früher in Einzelunter- 
suchungen 1)*)3) schrittweise (am Beispiel der Spaltung von U2%5 durch ther- 
mische Neutronen) entwickelten Vorstellungen noch einmal zusammengefaßt 
und systematisch dargestellt werden. 


Die Grundlage unserer Überlegungen bildet die Vorstellung, daß die 
Deformation der Bruchstücke (,,im Augenblick“ ihrer Trennung) von ihrer 
Schalenstruktur abhängt. (Es tritt eine minimale Deformation auf, wenn die 
Neutronenzahl eines der Fragmente magisch ist.) Aus diesem Bild ergeben sich 
als Konsequenzen zwei magische Effekte: 


1. Das Kernkraftpotential zwischen den beiden Fragmenten ist, bei festem 
Wert des Schwerpunktsabstandes s der Bruchstücke betrachtet, (betragsmäßig) 
vermindert für diejenigen Werte des Massenverhältnisses m = Ag: A, (As und 
A, Atomgewicht des schweren bzw. leichten Bruchstiicks), bei denen eines 
der Fragmente eine magische Neutronenzahl (N = 50 oder N = 82) besitzt. 


1) W. Brunner u. H.Paul, Ann. Physik (7) 6, 267 (1960). 
3% W.Brunner u. H.Paul, Ann. Physik (7) 7, 326 (1961). 
id TS 3) W. Brunner u. H. Paul, Ann. Physik (7) 7, 333 (1961). 
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2. Die Kurve der mittleren kinetischen Energie 7' der beiden Fragmente 
weist im Gegensatz dazu an denselben Stellen m Maxima auf. 

(Für N = 82 sind die Effekte jeweils stärker ausgeprägt als für N = 50, 
was damit zusammenhängt, daß im zweiten Fall die Wirkung der geringen 
Deformation des magischen Fragments durch eine stärkere Deformation des 
anderen Bruchstücks abgeschwächt wird.) 

Das Zusammenspiel der Effekte 1. und 2. gibt dann eine Erklärung der 
Spaltungsasymmetrie, wenn man den Spaltungsvorgang als Tunneleffekt 
(Durchdringung der aus Coulomb-Potential und Kernkraftpotential be- 
stehenden Potentialschwelle V) auffaßt, d.h. die relative Spaltungswahr- 
scheinlichkeit mit der Durchdringungswahrscheinlichkeit von V identifi- 
ziert. 

Die allgemeinen Überlegungen sollen an Hand der Verhältnisse bei der 
spontanen Spaltung von Cf? verdeutlicht werden. Für dieses Beispiel wird 
die Deformation der Fragmente, das zwischen ihnen herrschende Kernkraft- 
potential und ihre mittlere kinetische Energie berechnet (bezüglich U#® 
vgl. Ref.!)2)2). Man erkennt überall deutlich die magischen Effekte in 
guter Übereinstimmung mit den experimentellen Daten, während die Spaltung 
von Bi? (durch Deuteronen), die ebenfalls kurz diskutiert wird, durch das 
Fehlen solcher Effekte (im interessierenden m-Bereich tritt kein magisches 
Bruchstück auf) gekennzeichnet ist. 

Die Verhältnisse bei höheren Anregungsenergien werden kurz diskutiert, 
wobei gezeigt werden soll, daß die experimentell bekannte Abschwächung 
der Spaltungsasymmetrie bei Erhöhung der Anregungsenergie des sich spalten- 
den Kerns als Folge einer generellen Vergrößerung der Deformation der Frag- 
mente zumindest qualitativ verstanden werden kann. 


Diskussion der magischen Effekte 


Wir studieren den Vorgang der Kernspaltung bei niedrigen Anregungs- 
energien. Er erscheint im Rahmen des Trépfchenmodells des Atomkerns als 
ein Prozeß der Einschnürung des sich spaltenden Kerns, der bis zur völligen 
Trennung in zwei Bruchstücke fortschreitet. Die Wahrscheinlichkeit, daß 
die beiden Fragmente in einem bestimmten Massenverhältnis m gebildet 
werden, hängt nur schwach von m ab (s. z. B. Frankel u. Metropolis®)), 
und wir sehen sie im folgenden als konstant an. Wir machen weiterhin die 
idealisierende Annahme, daß bei einer Spaltung im Massenverhältnis m eine 
„ungeänderte Aufteilung‘‘ der gesamten Ladung des sich spaltenden Kerns 
auf die beiden Bruchstücke stattfindet (Zs : Zn, = m, wenn Zg und Z, die 
Kernladungszahl des schweren bzw. leichten Fragments bedeuten), was in 
guter Näherung der Forderung maximaler Energiefreisetzung (s. Ref.!)) ent- 
spricht. 


Jedes der Fragmente besitzt eine Anregungsenergie E, (¢ = 1, 2), die sich 
(im Bild des einfachen Trépfchenmodells) als Deformationsenergie äußert. _ 
Die dieser Deformationsenergie entsprechende Deformation der Bruchstücke 

‚im Augenblick“ ihrer Trennung ist als Grundlage der weiteren Überlegungen 
zu bestimmen. 


an 
4) S.Frankel and N. Metropolis, Physic. Rev. 72, 914 (1947). we 
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Die gesamte mittlere Anregungsenergie E, + E, der beiden Bruchstücke 
läßt sich im Prinzip sehr einfach aus der Energiebilanz für den Spaltungsvor- 
gang ermitteln. Auf Grund des Energieerhaltungssatzes ist die Energie des 
Gesamtsystems für einen Zeitpunkt vor dem Auftreffen des Einschußteilchens 
und für den „Augenblick‘‘ der Trennung der beiden Fragmente dieselbe, d. h. 
es gilt (da bei niederenergetischer Kernspaltung die Neutronenverdampfung 
erst nach der Spaltung einsetzt) die Beziehung 
(1) 
Hier bezeichnen M, M, und M, die Masse des Targetkerns, des Einschuß- 
teilchens und des Bruchstücks 7, E, ist die (eventuell vorhandene) Anregungs- 
energie des Einschußteilchens, x, seine kinetische Energie, V „ die potentielle 
Energie der sich gerade berührenden Fragmente und x ihre eventuelle kine- 
tische Energie „im Augenblick‘ ihrer Trennung. Da sich V , bei der endgül- 
tigen Trennung der beiden Bruchstücke völlig in kinetische Energie umsetzt, 
ist statt (1) zu schreiben 
wobei E, die mittlere Anregungsenergie des Bruchstücks i und 7 die mittlere 
kinetische Energie der beiden Fragmente bezeichnen soll. Zur Berechnung 
der Massen M, und M, werde die Fongsche°) Massenformel benutzt. Mit den 
experimentellen T-Werten folgt dann aus (2) unmittelbar die gesamte mittlere 
Anregungsenergie E = E, + E,. (Wegen der beschränkten Genauigkeit der 
Fongschen Formel und der Meßwerte für T sind die möglichen Fehler dieser 
Rechnung allerdings beträchtlich ®).) 
Unter der Annahme, daß die Aufteilung von E auf die beiden Bruchstücke 
im Verhältnis ihrer Massen stattfindet und daß die Anregungsenergie E, mit 
der Deformationsenergie D, des Trépfchenmodells identisch ist, ist aus E die 
Deformation zu bestimmen. In quadratischer Näherung gilt nach Bohr und 
Wheeler’), wenn man außerdem noch der (von den genannten Autoren nicht 
berücksichtigten) Abhängigkeit der Oberflächenspannung o von der Schalen- 
struktur des Kerns Rechnung trägt, 
2 2 
D, Aj® — aj (3) 


ı 


2 
(e elektrische Elementarladung, e = (=) w 50, ry 1,4- 10-8 em, 0, = 


krit. 
Wert der Oberflächenspannung im einfachen Trépfchenmodell). 
Der Deformationsparameter x, legt die Form des Bruchstücks i durch die 
Gleichung 
R,(#) = RP [1 +, P,(cos 9)] 


5) P. Fong, Physic. Rev. 102, 434 (1956). 
6) Eine Möglichkeit einer direkteren Bestimmung von E liegt in der Messung der Zahl 
der von den einzelnen Bruchstücken emittierten Neutronen in Abhängigkeit von m (8. 
z.B. Apalin’)). Hiermit ist sogar E, und E, einzeln zu ermitteln. Allerdings ist die 
Genauigkeit eines solchen Verfahrens ebenfalls nicht sehr groß. 
”),B.®. Ananun, W. Il. No6psiunn, B. TI. 3axaposa, Y.E.Kyrukos, 
cl. A. Mukaaaau, AD 8, 15 (1960). 
8) N. Bohr and J. Wheeler, Physic. Rev. 56, 426 (1939). 
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für den Kernradius fest (9 Polarwinkel, P, Legendresches Polynom). Die 
aus neueren Untersuchungen im Rahmen des Kollektivmodells (s. z.B. Alder?)) 
ermittelte Variation von o mit der Nukleonenkonfiguration des Kerns zeigt 
einen deutlichen magischen Effekt: o wächst auf etwa das Vierfache, wenn die 
Neutronenzahl N des Kerns magisch wird. (S. Abb. 1, in welcher die Kurve 
c/o, aufgetragen ist, wie man sie — wegen der starken Streuung der Meßwerte 
nicht ohne eine gewisse Willkür — aus den Daten von Ref.?) entnehmen 
kann; vgl. dazu auch Ref.?), Abb. 5, wo diese Daten mit angegeben wurden.) 

Die Deformation der Bruchstücke, in Abhängigkeit von ihrem Massenver- 
hältnis betrachtet, wird also im wesentlichen durch zwei Dinge bestimmt, 
erstens durch die Variation von E und zweitens durch die Variation von o 
mit m. Der zweite Effekt ist der ausschlaggebende; er bewirkt das Auftreten 
ausgeprägter tiefer Minima in der Kurve, in der x über der Neutronenzahl N 
eines Bruchstücks aufgetragen ist, bei den magischen Neutronenzahlen N = 50 
und N = 82, während der Verlauf von E(m) sich im wesentlichen auf die Lage 


24 


vO 50 60 70 80 90 100710 v0 50 60 70 80 90 100 710 120 
N 
Abb.1 Die Abhängigkeit der Ober- Abb.2. Aus der gesamten Deformations- 
flächenspannung o von der Neutronen- energie E der beiden Fragmente berechnete 
zahl N nach Ref.?) Es wurde mit Deformation eines Bruchstücks (in Abhän- 
0 = 3,83 - 10° MeV/cm? als Wert der gigkeit von seiner Neutronenzahl N) für 


Oberflichenspannung im einfachen 


Tröpfchenmodell gerechnet) die spontane Spaltung von Cf®?, Für E 


wurde dabei einmal ein (bezüglich m) kon- 
stanter Wert angenommen (Kurve a), und 
ay Les einmal wurde der aus Gl. (2) folgende E- 


und Form der Maxima der genannten Kurve auswirkt, deren charakteristi- 
schen Verlauf jedoch nicht beeinflußt. (Das ist einfach dadurch bedingt, daß 
die relative Änderung von E(m) wesentlich geringer ist als die von o(m).) 
Zur Illustration s. Abb. 2, in welcher für die spontane Spaltung von Cf2#2 


%) K.Alder, A. Bohr, T.Huus, B.Mottelson and A. Winther, Rev. Mod. 
Physics 28, 432 (1956). 
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die Werte «®(N) dargestellt sind, die einmal unter der Voraussetzung = 


const (in bezug auf m) und einmal aus dem E-Verlauf, der sich aus Gl. (2) 
ergibt, berechnet wurden. 

Wir können somit feststellen: Bezüglich ihrer Neutronenzahl magische 
Bruchstücke besitzen im Vergleich zu den anderen Fragmenten eine minimale 
Deformation. Dieser Tatbestand, der in einem engen Zusammenhang mit 
dem starken Anwachsen der Oberflachenspannung fiir magische Neutronenzahlen 
steht, erscheint uns als das Grundphänomen, aus dem heraus sich die wesent- 
lichen Züge des Kernspaltungsprozesses, nämlich der Verlauf der mittleren 
kinetischen Energie der Fragmente und die Asymmetrie in der Massenvertei- 
lung, erklären lassen. 

Wir studieren daher im folgenden systematisch die Auswirkungen der 
Schalenstrukturabhängigkeit der Deformation der Fragmente auf die einzelnen 
für die Spaltung charakteristischen Größen. 


B. Die Kernkraftwechselwirkung zwischen den Fragmenten 


Die beiden Spaltstücke stellen nach ihrer Bildung ein Zweiteilchensystem 
dar, dessen potentielle Energie aus der Wirkung der abstoßenden Coulomb- 
Kräfte und der anziehenden Kernkräfte resultiert. Während (entsprechend 
der Wignerschen Matrixtheorie der Kernreaktionen) die Wirkung der Kern- 
kräfte (infolge ihrer kurzen Reichweite) in den bisherigen Theorien zur Kern- 
spaltung stets vernachlässigt wurde, spielen nach unserer Meinung diese 
Kräfte eine wichtige Rolle beim Spaltungsprozeß. Das liegt vor allem daran, 
daß die oben diskutierte Schalenstrukturabhängigkeit der Deformation der 
Bruchstücke sofort zu einer Schalenstrukturabhängigkeit des von den Kern- 
kräften herrührenden Potentials V „ zwischen den beiden Fragmenten Anlaß 
gibt. 

Wegen der kurzen Reichweite der Kernkräfte können wir annehmen, daß 
V, für den Fall der Berührung der beiden Bruchstücke genähert gleich dem 


Potential — V x zweier sich berührender Kugeln ist. Wie wir in Ref.!), Anhang, 
theoretisch abschätzten, hängt Vx nur schwach von m ab und kann daher als 


m-unabhängig (V x © 135 MeV fiir die Bruchstücke von U) angesehen wer- 
den. Wenn man die Ortsabhängigkeit des Zweinukleonenpotentials in der 
Form e”*" (r Nukleonenabstand, u = M, c/h mit M, = Masse des n-Mesons) 
annimmt, ist für V, genähert zu schreiben 

V (m;8) = 3,(m)). (5) 
Hier bedeutet s den Schwerpunktsabstand der beiden Bruchstücke und s, 
den Wert von s, bei dem gerade eine Berührung der beiden Bruchstücke statt- 
findet. (Gl. (5) wurde in Ref.!) für den Fall zweier Kugeln näherungsweise 
hergeleitet; wir führen einen weiteren, durch die Kurzreichweitigkeit der Kern- 
kräfte gerechtfertigten Näherungsschritt durch, indem wir sie auf den Fall 
zweier deformierter Kerne übertragen.) 

Von der Deformation der Bruchstücke nehmen wir an, daß diese stets als 
eine Dehnung in Richtung der Verbindungslinie der Schwerpunkte der beiden 
Fragmente erfolgt, d.h. # in Gl. (4) ist von dieser Verbindungslinie aus zu 
zählen. Dann ist offenbar ih. 


(m) = [1 +0, (m)] RO + [1 + &,(m)] RS, 
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und Gl. (5) läßt sich umschreiben in 
Ve (m; s) =— V . {8-[1+ a, (m)] RY-[1+ (m)] RE} (7) 


Man erkennt hieraus, daß das Kernkraftpotential, bei festem Schwerpunkts- 
abstand s betrachtet, eine (hauptsächlich über &, und &,) mit m variierende 
Größe ist. Der magische Effekt im «-Verlauf führt zu einer Verminderung von 
V,(bei festem s) für die Werte von m, bei 2 
denen ein Bruchstück eine magische Neu- I W=50 
tronenzahl besitzt. Diese ,,magischen m- Fi 

Werte“ hängen offenbar von der Masse des  ~ jo 
Targetkerns ab, und zwar gilt für U5 3 | 
m = 1,3 (entspr. N = 82) und m = 1,9 (ent- Sao 
spr. N = 50), für Cf?® m = 1,1 (entspr. 120- 
N = 82) und m = 2,1 (entspr. N = 50), und 100 
für Bi2® m = 1,9 (entspr. N = 82) sowie PR 
m = 1,5 (entspr. N = 50). In Abb. 3 ist der 7 
Wert von V , an der Stelle s = 14,0- 10-13 cm a Aa RIO 
fir den Fall der spontanen Spaltung von 
Cf in Abhängigkeit von m aufgetragen, Abb. 3. Kernkraftpotential V; 
wie er aus den Überlegungen des AbschnittsD (Betrag) für die spontane Spaltung 
folgt (s. u.). Man erkennt in der Tat ausge- von Cf*? in Abhängigkeit vom 
prägte Minima bei den magischen m-Werten. 
(In der entsprechenden Kurve für fehlt 10cm) nach GI. (2) und 
das Minimum fiir m = 1,9, wenn man das den a-Werten der Abb. 5 

an dieser Stelle liegende zweite Maximum 

(entsprechend N = 50) im T-Verlauf nicht berücksichtigt; es tritt dagegen 
auf, wenn man das genannte Maximum in Rechnung stellt. S. Ref.?)) 

Die bei der Deformation der Bruchstücke in Erscheinung tretenden ma- 
gischen Effekte geben somit zu magischen Effekten im Verlauf von V, (bei 
festem s) Anlaß, die sich auf das aus Coulomb-Potential und Kernkraft- 
potential sich zusammensetzende Gesamtpotential 


8 


(8) 


übertragen. Eine direkte Folge dieses Sachverhaltes ist das Auftreten magi- 
scher Effekte im Verlauf der mittleren kinetischen Energie T der Fragmente, 
wie im folgenden näher ausgeführt werden soll. 


C. Die mittlere kinetische Energie der Fragmente 


Der Zerfall des Zweiteilchensystems der beiden Bruchstücke findet über 
ein virtuelles Niveau der Energie e statt. (In e ist die Anregungsenergie der 
beiden Bruchstücke nicht enthalten!) Da sich e nach dem Auseinanderfliegen 
der beiden Fragmente völlig in kinetische Energie umsetzt, ist e mit der mitt- 
leren kinetischen Energie T der Bruchstücke zu identifizieren. Die virtuellen 
Niveaus sind dadurch ausgezeichnet, daß die zugehörige Wellenfunktion y(s) 
des Zweiteilchensystems für s — oo eine reine auslaufende Welle repräsentiert. 
Zur Berechnung von ¢ benötigt man offenbar die Kenntnis der Randbedingung 
für y an der Stelle s,, die sich im Prinzip aus der Lösung des A-Teilchen-Pro- 
blems ergibt, sich aber praktisch kaum berechnen läßt. 
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Obwohl wir die Randbedingung explizit nicht kennen, scheint es uns doch 
möglich, eine wesentliche Aussage über sie zu machen: Da für s = s, die phy- 
sikalische Situation (Kontakt der beiden Fragmente) unabhängig von m jeweils 

P dieselbe ist, halten wir die Annahme, daB auch die Randbedingung von m 
nicht wesentlich abhängt, für plausibel. Mit dieser Annahme läßt sich dann 
€(m) bereits bis auf eine m-unabhängige additive Konstante (deren Wert eben 
durch die genaue Form der Randbedingung bestimmt ist) berechnen, wie 
wir gleich zeigen wollen. 


Für die Festlegung von e sind hauptsächlich die V-Werte im (von s = s, aus) 
ansteigenden Teil der Potentialschwelle maßgebend, während die genaue Form 
des Abfalls von V(s) bei größeren s-Werten nicht sehr wesentlich eingeht. 
Man kann daher in dem interessierenden s-Bereich das erste Glied in (8) nach 
Potenzen von s’ = s — s, entwickeln und die Entwicklung mit dem zweiten 
Glied abbrechen. Am (8) folgt dann genähert 


q a 20 2 


®,) 0 


Im zweiten Term ist die Variation von Z, Z,/32 mit m in guter Näherung zu 
vernachlässigen. Daher ist (genähert) die Form des Potentialverlaufs für alle 
Werte m die gleiche, bei Übergang von einem m, zu einem m, findet lediglich 
eine Verschiebung der V (s’)-Kurve um den Wert Z, (m,) Z,(m,) - 89 (m,) -e2 — 
Z, (my) Z,(m}) - so (m}) e? statt. Aus der obigen Annahme der m-Unabhängig- 
keit der Randbedingung für y folgt dann sofort, daß das virtuelle Niveau 


ebenfalls mit verschoben wird, d.h. es gilt 2 
Z,(m,)Z,(m,)e?  Z,(m,) Z,(m,) e* 
€ (m,) €(m) = 8 so(m, ) ry (10) 
0 


wobei A eine m-unabhängige Konstante bezeichnet, deren genauer Wert nur 
über die Lösung des A-Teilchen-Problems zu beschaffen wäre. (In Ref.) 
haben wir eine Abschätzung von A dadurch versucht, daß wir die durch die 
vermutliche Existenz eines C!2—C!2.,,Molekiils‘‘ [vgl. die Experimente von 
Almgvist!P) et al.] nahegelegte Randbedingung p(s = 8) = 0 verwendet 
haben. Auf diese Weise kann aber nur die Größenordnung von A erhalten 
werden.) 


ö Die Gleichung (11) beschreibt aber bereits vollständig die Variation von 
€ mit m, und man erkennt deutlich die Auswirkung der magischen Deformations- 
effekte auf den Verlauf von ¢ und damit auch (s. o.) von T': Bei den oben defi- 
nierten ,,magischen m-Werten‘‘ ist der Nenner des ersten Gliedes in (11) wegen 
der geringen Deformation des magischen Bruchstücks besonders klein, damit 
tritt für solche m-Werte (die wir oben für U*>, Cf?5? und Bi? explizit ange- 
geben haben) ein Höcker in der T-Kurve auf. Das ist der experimentell 


10) D.A. Bromley, JA. Kuehner and E. Almqvist, Physic. Rev. L. 4, 365 (1960). 
E. Almqvist, D. A. Bromley, and J.A.Kuehner, Physic. Rev. L. 4, 515 (1960). 
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doch f hekannte magische Effekt im T-Verlauf (s. z. B. Protopopov!!) et al.), der 
phy- hiermit eine Erklärung über die Schalenstrukturabhängigkeit der Deformation 
weils | der Fragmente gefunden hat. In Abb. 4 er- 


N M kennt man, wie gut Gl. (11) die experimen- ” 
dann § zellen 7'-Werte für Cf252 (nach Fraser und 190 7 
eben § Milton!2)) wiedergibt (vgl. auch die theo- 780 4 
Wie retischen Werte für in Ref.3)). (Es 
wurde in (11) RP + RE = 11,9- 10-8 cm = 
aus) | gesetzt und 4 aus der Forderung, daß die oa ne 
‘orm theoretische Kurve an der Stelle m = 1,0 Pn 
zeht, | mit dem experimentellen Wert überein- 10 12 14 16 18 20 22 


m 
Abb. 4. Die mittlere kinetische 
Energie der beiden Fragmente für 
die spontane Spaltung von Cf*?. 


nach | stimmt, zu A = 40,71 MeV ermittelt. Für x, 
eiten | und x, wurden die im folgenden Abschnitt D 
berechneten und in Abb. 5 dargestellten 


Werte benutzt.) Kurve: theoretische Werte nach 
9 Damit verbleibt die Diskussion des Ein- Gl.(11)mit R?+ R2=11,9-10-%em. 
(9) ) = 40,71 MeV und den x-Werten 


flusses der magischen Effekte im Verlauf 


der Abb. 5. 
von Vx und T auf die Spaltungswahrschein- 


experimentelle Werte nach 


| Jichkeit Ref.” 
. 
- alle ‘ 
glich D. Die Spaltungswahrscheinlichkeit j 
Gestiitzt auf die Experimente von Bigham) und Almqvist) nehmen 


1818- | wir an, daß beim Spaltungsprozeß nur sehr wenig Ausgangskanäle vorhanden 
veau | sind, so daß die relative Spaltungswahrscheinlichkeit durch das Produkt aus 
der Wahrscheinlichkeit der Herausbildung der beiden Spaltstücke (infolge 
Einschnürung) und der Durchdringungswahrscheinlichkeit der Potential- 
(10) | schwelle (8) gegeben ist. Da erstere als m-unabhängig angesehen werden kann 
(s.0.), läßt sich somit (in WBK-Näherung) für die relative Spaltungswahr- 
scheinlichkeit w der Ausdruck 


(11) 
w= exp s) — eds 12 
- aufschreiben. (vy reduzierte Masse des Zweiteilchensystems, s, und s, klassische 
bx Umkehrpunkte der Bewegung.) Damit ist nun sehr einfach der Einfluß der 
bs magischen Effekte zu studieren, die sich im Potentialverlauf (s. Abschnitt B) 
2 und in der Höhe des virtuellen Niveaus, in der die Durchdringung statt- 
ten 


findet, (s. Abschnitt C) äußern. Wenn man zunächst einmal von dem Kern- 
kraftanteil V, der Potentialschwelle V absieht, führt der N = 82 entspre- 
von | chende Höcker im 7-Verlauf zu einer Verkleinerung der über dem Tunnel 
ions- | liegenden ‚„‚Bergkuppe‘“ (durch Anheben des Tunnels) und damit zu einer Ver- 
defi- | größerung der Durchdringungswahrscheinlichkeit für das N — 82 entspre- 
egen | chende Ma8senverhältnis. Man würde somit (bei Abwesenheit von V ,) asym- 
amit | metrische Spaltung mit einem Maximum der Spaltungswahrscheinlichkeit für 


nge- 
1, A.H. Tporononos.H.A.Bapanos, BIL Diiemonr. 
KITM 36, 1932 (1959). 
2) J.S. Fraser and J.C.D. Milton, Proceedings of the Second United Nations 
960). International Conference on the Peaceful Uses of Atomic Energy, Geneva, P/199 (1958). 
960). 13) C,B. Bigham, Chalk River report CRP—642—A, 111. 


4) E. Almqvist, Chalk River report, CRP—642—A, 142. 
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das kleinste magische m als Folge des magischen Effektes im 7'-Verlauf er- 
warten. (Der zweite, N = 50 entsprechende Hocker der kinetischen Energie 
T fällt in einen m-Bereich, in dem die Spaltungswahrscheinlichkeit schon stark 
abgeklungen ist; außerdem ist er schwächer ausgeprägt als der erste, so daß 
er für die Spaltungswahrscheinlichkeit keine Rolle spielt.) Betrachten wir den 
zusätzlichen Einfluß der Kernkraftwechselwirkung! Sie bewirkt generell eine 
Verkleinerung des Berges über dem Tunnel, also eine Erhöhung der Durch- 
dringungswahrscheinlichkeit. Für ein magisches m ist V „ (bei festem s) aber 
(betragsmäßig) besonders klein, daher ist die Vergrößerung der Durchdringungs- 
wahrscheinlichkeit gerade für magische m-Werte keiner als für die anderen 
Werte des Massenverhältnisses. Auf diese Weise kommt es zu einer Ver- 
schiebung des wahrscheinlichsten Massenverhältnisses zu größeren Werten 
(bei U5 von 1,3 zu 1,45; bei Cf*52 von 1,1 zu 1,3). 

Doch nicht nur die Lage des Maximums der Massenverteilung, sondern 
deren gesamter Verlauf wird durch die entwickelte Theorie richtig wiederge- 
geben. Der wesentliche Erfolg unserer Überlegungen besteht gerade darin, 
daß die experimentellen Werte sowohl der Spaltungswahrscheinlichkeit w als 
‚auch der mittleren kinetischen Energie T der Bruchstücke mit den gleichen 
Werten der Deformationsparameter der Fragmente reproduziert werden, wobei 
der charakteristische Verlauf dieser Werte seinerseits durch die aus dem Experi- 
ment folgende Abhängigkeit der Obe ‘lichenspannung o von der Schalen- 
‚struktur der Kerne physikalisch völlig .. vständlich ist und außerdem sehr gut 
a s der im Grundzustand vorhandenen(ebenfalls experimentell ermittelten) 

E Tendenz zur Deformation übereinstimmt (s. Ref.?)). 


Zusammenfassend können wir damit festhalten, daß die entwickelte Vor- 
‚stellung ein zwangloses physikalisches Verständnis der Spaltungsasymmetrie 
auf der Grundlage einer schalenstrukturabhängigen Deformation der Spalt- 
stiicke „im Augenblick‘ ihrer Entstehung liefert. Im besonderen ist klar, wie 
die Abhängigkeit des wahrscheinlichsten m-Wertes von der Masse des Target- 
_kerns zustande kommt. Sie ist primär durch den Umstand bedingt, daß das 
magische Massenverhältnis von der Masse des Targetkerns abhängt. 
Man kann nun auch die experimentell bekannte Tatsache der symme- 
trischen Spaltung leichterer Targetkerne, etwa von Bi?® durch Deuteronen, 
zwanglos verstehen: Da in dem interessierenden m-Bereich auf Grund der 
_ kleinen Masse des Targetkerns Bi2® kein magischer Wert des Massenverhält- 
-nisses liegt, (das kleinste magische m hat den Wert m = 1,5, s. o.), hat man 
nach unseren Vorstellungen keinerlei magische Effekte in den Spaltungsdaten 
_ (für den interessierenden m-Bereich) zu erwarten: Es wird weder ein magischer 
Hocker im T-Verlauf noch eine Einsenkung im Verlauf von V „ (für festes s) 
_ auftreten, damit entfällt die Ursache für eine mögliche Spaltungsasymmetrie ®). 
| Es sei noch erwähnt, daß Gl. (12) die genaueste Berechnung der Deforma- 
_ tionsparameter x, und x, (aus den experimentellen Daten für w tind T) er- 
laubt, da sich eine Änderung der x-Werte (wegen ihres Auftretens im Exponen- 


15) Man muß allerdings beachten, daß es sich hier nicht um einen Spaltungsvorgang 
bei niedrigen Anregungsenergien handelt, den wir unseren bisherigen Überlegungen 
_ immer zugrunde gelegt haben. (Die kinetische Energie der eingeschossenen Deuteronen 
_ beträgt 22 MeV.) Die Tatsache der symmetrischen Spaltung von Bi®® wird daher auch zu- 
_ mindest zum Teil durch den unten diskutierten Einfluß der höheren Anregung auf die 
 Spaltungswahrscheinlichkeit verursacht sein. 
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ten der e-Funktion) stark auf die Größe von w auswirkt. Aus Gl. (12) kann man 
offenbar die Summe a, RP? + «, RY ermitteln, die in dem interessierenden 


m-Bereich praktisch durch = (R? + RS) (x, + &,) ersetzbar ist. (Ein geeignetes 


Näherungsverfahren zur Auflösung von Gl. (12) nach V, - exp {u (& R + 
x, RZ)} wurde in Ref.!) angegeben.) Durch die Forderung, daß sich die Anre- 
gungsenergien (3) der beiden Bruch- 


stücke wie ihre Massen verhalten, sind 4 9 

dann x, und x, einzeln festgelegt. In „ 2] 

Abb.5 sind die so erhaltenen Werte 02 0 7- 

für ein bei der Spaltung von Cf ent- S "eo x 

stehendes Bruchstiick aufgetragen; sie 8- . 
stimmen in ihrem charakteristischen a 4 
Verlauf mit den aus der gesamten Defor- 57 : 

mationsenergie E berechneten und in 31 

Abb. 2 dargestellten Werten befriedigend : ] 


überein, woraus die innere Konsistenz 

unserer Überlegungen ersichtlich ist. 40 50 60 70 80 90, 100 110 
. . Y . 

Eine quantitative Übereinstimmung Abb. 5. Die Deformation eines Bruch- 

kann wegen der z. T. recht groben stücks in Abhängigkeit von seiner Neu- 

Näherungen, die bei der Berechnung tronenzahl für den Fall der spontanen 

der letztgenannten Werte x? gemacht Spaltung von Cf®?, berechnet aus der 


, 1 Spaltungswahrscheinlichkeit und der 
wurden (Fongsche Massenformel™), mittleren kinetischen Energie nach 


Willkür in der Wahl des o-Verlaufs, Be- Gl. (12) 
schränkung auf quadratische Glieder im 
Ausdruck (3) für die Deformationsenergie), nicht erwartet werden. — Mit den 


x-Werten der Abb. 5 wurde oben nach Gl. (7) V„(m;s = 14- 10-8 em) be- 
rechnet und in Abb. 3 aufgezeichnet. 

Endlich seinoch kurz auf die Frage eingegangen, inwieweit die experimentell 
bekannte Abschwächung der Massenasymmetrie bei Erhöhung der Anregungs- 
energie des sich spaltenden Kerns auf der Grundlage unserer Vorstellungen 
über die Deformation der Spaltungsfragmente (zumindest qualitativ) zu ver- 
stehen ist. 


Spaltung bei höheren Anregungsenergien 


Die Problematik der Spaltung bei höheren Anregungsenergien (wir denken 
an kinetische Energien des eingeschossenen Teilchens größer als etwa 10 MeV) 
besteht zunächst darin, daß neben der ,,first chance‘‘-Spaltung (Spaltung ohne 
vorherige Neutronenemission) ein beträchtlicher Anteil an ,,second-chance*- 
Spaltung (Emission eines Neutrons vor der Spaltung), ,,third-chance‘‘-Spaltung 
usw. auftritt (vgl. z.B. Leachmann!”)). Die folgenden Überlegungen be- 
ziehen sich auf die ,,first-chance‘‘-Spaltung. (Die anderen Anteile repräsen- 
tieren offenbar Spaltungsvorgänge bei niedrigen Anregungsenergien.) Aus der 
experimentell bekannten Tatsache, daß sich die kinetische Energie der Frag- 


16) Die Fongsche Massenformel liefert die absolute Größe der Deformationsenergie 
und damit auch der Deformationsparameter nur relativ ungenau, woraus sich der große 
Unterschied der Absolutwerte der Deformationsparameter in den Abbn. 2 und 5 erklärt. 

7) R. B. Leachman, Proceedings of the Second United Nations International 
Conference on the Peaceful Uses of Atomic Energy, Geneva, P/2467 (1958). 
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mente bei Erhéhung der kinetischen Energie x, der EinschuBteilchen fast 
nicht ändert (s. z. B. Frisch!8)), ziehen wir den Schluß, daß die zusätzlich 
zugeführte Energie völlig in Anregungsenergie der Fragmente umgesetzt wird. 
Wenn wir wieder Anregungsenergie mit Deformationsenergie identifizieren 
(oder doch zumindest annehmen, daß eine Vergrößerung der Anregungsenergie 
zu einer Vergrößerung der Deformation führt), findet somit bei einer Steigerung 
von x, ein Anwachsen der Deformationsparameter x, und x, statt, und damit 
kommt es zu einer starken Vergrößerung des Kernkraftpotentials V, an der 
Stelle des Tunnels. Das bedeutet aber, daß die Potentialschwelle V stark 
erniedrigt wird. Da sich die Lage der virtuellen Niveaus e (die mit 7 identisch 
sind) praktisch kaum ändert, findet de facto für höhere Anregungsenergien 
keine Durchdringung mehr statt, sondern 


4 | _ Anregungsener gie ken die Potentialschwelle wird „überflogen‘“ 
| J (In Abb. 6 sind die Verhaltnisse schema- 
4 V 4 tisch dargestellt.) 
|= 1 Es ist klar, daß die Unterschiede in 
der Spaltungswahrscheinlichkeit für ver- 
schiedene Werte des Massenverhältnisses 
#3 Pr um so kleiner werden, je stärker die 
4 egungsenergie gro 
Schwelle erniedrigt wird, d. h. letzten 
Endes, je stärker der sich spaltende Kern 
gs 5 


angeregt ist. Diese Tendenz in der Ände- 
tialschwelle bei Erhöhung der kine- Anregungsenergie 
tischen Energie der Einschußteilchen stimmt genau mit der experimentell fest- 
über eine Vergrößerung der Defor- gestellten überein (s. z. B. Leachman!"’)). 
mation der Fragmente (schematisch) Wir sind somit auch zu einem qualitativen 

Verständnis der experimentell bekannten 
Abschwächung der Spaltungsasymmetrie (das Verhältnis der Spaltungswahr- 
scheinlichkeiten für symmetrische und wahrscheinlichste Spaltung vergrößert 
sich mit zunehmender Anregungsenergie) gelangt. 

Eine quantitative Diskussion erscheint schon aus mathematischen Gründen 
sehr kompliziert. Zur Berechnung der Wahrscheinlichkeit des ‚‚Überfliegens“ 
der Potentialschwelle ist die WBK-Näherung nicht geeignet und müßte durch 
eine wesentlich genauere und kompliziertere Methode ersetzt werden. Weiter- 
hin sind für große Deformationen die Gleichungen (3) für die Deformations- 
energie und (5) für das Kernkraftpotential unzulässige Näherungen; an ihre 
Stelle müßten wesentlich kompliziertere Formeln treten, die auch noch weitere 
Deformationsparameter (die Koeffizienten von P,, P, usw. in der Entwick- 
lung (4), die nun auch nicht mehr mit dem zweiten Glied abgebrochen werden 
darf) enthalten müßten. Abgesehen davon erscheint nun auch die Identifika- 
tion von Anregungsenergie und Deformationsenergie des Trépfchenmodells 
recht problematisch; eine genauere Untersuchung wäre im Rahmen des 
Kollektivmodells durchzuführen. Wir begnügen uns mit der oben gegebenen 
qualitativen Diskussion. 

Abschließend sei noch auf ein Experiment hingewiesen. das eventuell 
einen unabhängigen Schluß auf die Deformation der Fragmente zu ziehen ge- 
stattet. Wir denken hierbei an die Messung der Winkelverteilung der von 


Abb. 6. Die Erniedrigung der Poten- 


8) O. R. Frisch, Progress in Nuclear Physics Vol. II, London 1952, S. 167. 


es 


ei 
A 
| er 
| pe 
In 
m 
+ 
% 
= 
= 
| 
| 
| 


Dr W. Brunner u. H. Paul: Zur Theorie der Kernspaltung. IV 


einem Bruchstück emittierten y- Quanten in Abhänigkeit vom Massenverhältnis. 
Allerdings wird die theoretische Analyse derartiger Meßergebnisse dadurch 
erschwert, daß Bruchstücke möglicherweise zusätzlich eine große Spinkom- 
ponente senkrecht zu ihrer Flugrichtung besitzen. Vergleiche hierzu etwa die 
Interpretation des Experiments von Hoffmann) (Messung der über m ge- 
mittelten y-Winkelanisotropie) durch Strutinski?®). 


Herrn Prof. Dr. G. Richter sind wir für sein ständiges förderndes Inter- 
esse an unseren Überlegungen sehr zu Dank verpflichtet. 


19) M.M. Hoffmann. Bull. Am. Phys. Soc. 3, 6 (1958). y. x 
2) V.M.Strutinski, JETP 37 (10), 613 (1960). Ld 


Zeuthen bei Berlin, Kernphysikalisches Institut der Deutschen Akade- 
mie der Wissenschaften. 


_ Bei der Redaktion eingegangen am ‚11. November 1960. 
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Sondenmessungen bei hé jheren Stromdichten 


Mit 17 Abbildungen 


4 
[ 
u | 
Inhaltsübersicht 


Es werden in einer Argonentladung bei 8 em Rohrdurchmesser und 0,1 Torr 
Fülldruck die Plasmaparameter höherer Stromdichten mittels Sonden unter- 
sucht. Die Entladung wird stationär betrieben. Mittels Sonden kann die 
Elektronentemperatur und die Trägerkonzentration über den Radius des 
tohres ermittelt werden. Eine Messung nach der Zweisondenmethode liefert 
die Werte des Längsgradienten. Die Auswertung der Sondencharakteristiken 
wird nach zwei verschiedenen Methoden durchgeführt und zeigt, daß die Me- 
thode der Auswertung des Ionenstromteils nach Kagan im Rahmen der all- 
gemein bei Sondenmessungen zu erwartenden Genauigkeit im Vergleich mit 
der Langmuir-Methode befriedigende Übereinstimmung liefert. 

Die Trägerkonzentration läßt sich gut durch eine Besselfunktion an- 
nähern, und die von Fetz eingeführte Randbedingung, die die Aufladung 
der isolierenden Rohrwand berücksichtigt, gibt die Messungen im Rahmen der 
hier möglichen Genauigkeit richtig wieder. Aus der Beziehung J = e N b_ E, 
wird durch Integration der Trägerkonzentration über den Rohrquerschnitt 
die Beweglichkeit berechnet. Die erhaltenen Werte stimmen gut mit den bei 
dreißigmal kleineren Stromdichten von Wojaczek und Rademacher ge- 
fundenen Werten überein. Weder in der Trägerverteilung noch im Radial- 
potential gelingt es, bei den gemessenen Entladungsströmen das Eigenmagnet- 
feld nachzuweisen. Die in der Literatur auftretenden Formeln für den Ionen- 
sättigungsstrom, die sich durch einen konstanten Faktor unterscheiden, werden 
überprüft. Eine Entscheidung zwischen den erhaltenen Beziehungen nach 
Kagan, Bohm, Burhop und Massey und nach Langmuir ist auf Grund 
unserer Messungen nicht möglich. 


Sondenmessungen bei hohen Stromdichten, die der Bestimmung der 
Entladungsparameter in der positiven Säule einer stationären Niederdruck- 
entladung dienen, wurden bisher nur in einigen wenigen Fällen durchgeführt. 
Meist handelt es sich dabei um Entladungen in Hg-Dampf mit einer Queck- 
silbersumpfkathode. So untersucht z.B. Fetz!) eine Hg-Entladung im Druck- 
bereich 1.103 — 1,6 - 10 Torr bei Strömen bis zu 80 A und Stromdichten 


1) H. Fetz, Ann. Physik 40, 579 (1941). 5 
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bis 0,26 A/em?. Thonemann und Cowhig?) führen Versuche bis 100 A und 
Stromdichten von 4,5 A/cm? ebenfalls in Quecksilberdampf aus. In einer Arbeit 
von Danders?) werden Störeinflüsse eines äußeren Magnetfeldes auf eine 
Hg-Entladung gemessen. Die Entladungsstromstärken betragen bis zu 300 A 
in einem Druckbereich von 2,5 - 10% — 1. 10 Torr. Die Stromdichte liegt 
aber wegen des großen Rohrdurchmessers bei etwa 0,9 A/cm?. 

Entladungen hoher Stromdichte in Edelgasen sind bisher nur in Kapillaren 
untersucht worden oder dienten lediglich der Bestimmung der Stabilitäts- 
bereiche der Entladungssäule. Bei Kapillarentladungen sind jedoch die Ein- 
flüsse der Wand auf die Säule sehr groß und können Störungen der Entladung 
hervorrufen. Ferner wird das Messen mit Sonden problematisch, da das Ver- 
hältnis von Rohrdurchmesser zu Sondendurchmesser ungünstig wird und 
Störungen des Entladungskanals durch die Sonde zu erwarten sind. Sonden- 
messungen in Kapillarentladungen, verbunden mit spektroskopischen Unter- 
suchungen, haben Sacharowa und Kagant) durchgeführt. Die Autoren ver- 
wenden Kapillaren von 0,4 cm Innendurchmesser und untersuchen mit Zylin- 
dersonden von 0,02 em Durchmesser im Bereich von 1,2 - 10 — 1,2 - 10-1 Torr 
Entladungen in Hg, Ar, Xe und Kr. Die Stromdichte wird von 0,2—6,0 A/cm? 
geändert. Instabilitätsuntersuchungen mittels Sonden sind uns von Timo- 
fejewa und Granowski®) bekannt. In einem Druckbereich von 104 bis 
10 Torr werden hier zur Anzeige der Pinchinstabilitäten in stromstarken 
Entladungen von Hg, H,, Ar und Kr Sonden benutzt. Bei Rohrdurchmessern 
von 7 cm erreichen die Autoren ungefähr 45 A/cm?. Hier gleichfalls angeführte 
Entladungen bei höheren Stromstärken werden nicht stationär betrieben. 

Aufgabe unserer Untersuchungen ist es, für die stationäre Argonentladung 
bei Stromdichten zwischen 0,3 und 1,3 A/em? Sondenmessungen durchzu- 
führen, um den Einfluß der Stromdichte auf Elektronentemperatur, Elektronen- 
beweglichkeit und Trägerdichte festzustellen. —— 
pe 


Versuchsapparatur und Meßmethode 


Für die Untersuchungen verwenden wir ein aus Hartglas hergestelltes 
Entladungsrohr mit 8cm Innendurchmesser, dessen Länge etwa 100 cm be- 
trägt. Als Anode dient ein Zylinder aus Molybdänblech. Im Abstand 18,2 cm 
voneinander sind zwei Zylindersonden im mittleren Bereich des Rohres ange- 
bracht. Ihre Lage kann mit Hilfe eines Magneten in radialer Richtung ver- 
ändert werden. Die Sonden bestehen aus Wolframdraht von 0,06 em Durch- 
messer und haben eine freie Länge von 0,4cm. Eine Kapillare vom Außen- 
durchmesser 0,1 em bedeckt den übrigen Teil des Sondendrahtes. Für das Rohr 
wird eine Glühkathode vom Typ der Vorratskathoden verwendet. Gegen das 
lonenbombardement aus der Entladung haben wir keine besonderen Schutz- 
maßnahmen getroffen. Im untersuchten Stromstärkebereich treten keine 
Schichtungen in der Entladung auf. 

Das Entladungsrohr haben wir in der üblichen Weise gereinigt, gepumpt, 
ausgeheizt und formiert. Das Rohr wurde mit spektralreinem Argon vom 


2) P.C. Thonemann u. W.T.Cowhig, Proc. physic. Soc. 64, 345 (1951) 
3) M. Danders, Z. angew. Physik 9, 223 (1957). 

‘) B.M. 3Baxaposa, u W}W.M. Karan, &. Ont. u Cnextp. 1, 5, 627 (1956). 
5) 1.1. Tamogpeesa, B.!l. Tpanosernt, SKIT® 30, 477 (1956). 
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VEB Sauerstoffwerk Berlin gefüllt, wobei ein Gasdruck von 0,1 Torr ein- 
gestellt wurde. Wir betrieben das Rohr im abgeschmolzenen Zustand. 
> Wegen der thermischen Belastung der Rohrwandung und der Kathode 
bei höheren Stromstärken wird das Rohr impulsmäßig getastet. Die Impuls- 
zeiten sind mit etwa 0,6 sec so bemessen, daß sich in jedem Falle eine stationäre 
Entladung ausbilden kann. Mit Hilfe eines Oszillographen im Entladungskreis 
ist eine Kontrolle der Stationarität möglich. Die Sonde arbeitet ebenfalls im 
Impulsverfahren mit Zeiten zwischen 0,06 und 0,25 sec, einerseits um eine 
thermische Überlastung der Sonde und die dadurch hervorgerufenen Fehler 
auszuschließen, andererseits um Meßfehler durch Änderung der Arbeitsfunk- 
tion der Sonde, hervorgerufen durch Verunreinigungen, zu vermeiden. Eine 
eingehende Beschreibung der Arbeitsweise der ganzen Meßanordnung ist in 
6) angegeben. 

Vor jeder Messung wird die Sonde mittels eines Teslatrafos gereinigt. 
Die Charakteristiken werden oszillographisch aufgezeichnet und photogra- 
_ phiert. Die erhaltenen Bilder werden vergrößert und die Maßeinteilung der 
Achsen wird eingetragen. Durch Anwendung der Langmuirschen Sonden- 
theorie lassen sich aus unseren Messungen Elektronentemperatur, Plasma- 
potential und Elektronenkonzentration bestimmen. Da infolge hoher Ent- 
ladungsstromdichten Schwierigkeiten bei der Aufnahme des oberen Knickes 
der Sondencharakteristik auftreten können, wird zusätzlich zur Auswertung 
die Methode von Kagan und Mitarbeitern?) herangezogen. Dadurch lassen 
sich neben der Elektronentemperatur die anderen Parameter bereits aus dem 
_ lonenteil der Charakteristik berechnen und mit den nach der Langmuir- 
Methode gewonnenen Werten vergleichen. Die Längsfeldstärke E, wird mittels 
der Doppelsondenmethode bestimmt, wobei die Sonden einen Abstand von 
18,2em haben. Aus den gemessenen Kurven läßt sich für die stromlosen 
Sonden die Spannung zwischen diesen ablesen, woraus sich die Längsfeld- 


stärke E, ergibt. 


iis 


Die Sondentheorie und ihre Verfeinerung 


ACR 


In der Mehrzahl der Fälle wird bisher bei Sondenmessungen die von 
Langmuif und Mott-Smith®) gefundene Theorie der Sondencharakteristik 
als Grundlage angenommen. Nach dieser Theorie nimmt die Sonde im Bereich 
stark negativer Vorspannungen nur die Ionen aus dem Plasma auf, während 
die Elektronen nicht in die Nähe der Sonde gelangen. Um die Sonde bildet 
sich eine Raumladungsschicht aus, und man setzt voraus, daß das elektrische 
Feld von der Sonde außerhalb dieser Schicht Null ist. Wird die Sondenspan- 
nung positiver gemacht, so können die schnellen Elektronen gegen das Poten- 
tial anlaufen, was zu einer Verkleinerung der Dicke der Raumladungsschicht 
führt. Die Schicht verschwindet schließlich ganz, wenn der Wert des Plasma- 
potentials erreicht wird. Damit ist auch der Sättigungsbereich des Elektronen- 
stromes erreicht. Oberhalb des Plasmapotentials bildet sich dann, entsprechend 
zu den Vorgängen im lonenstromgebiet, eine Raumladungsschicht entgegen- 
gesetzter Ladung aus. Unter diesen Voraussetzungen ergibt sich, daß in den 


6) E.Hayess, Exp. Techn. Physik 8,4, 145 (1960) 
7) |}. M.Karan, Becrn. JIPY Cep. Ons. u Xum. No. 4, 63 (1957). 
8) I. Langmuir u. H. M. Mott-Smith, Physic. Rev. 28, 727 (1926). 
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Gebieten der Sättigungsströme die Ladungsträger in Form von Diffusions- 
strömen zur Sonde fließen. 

Theoretische Untersuchungen der Raumladungsschicht von Bohm, Bur- 
hop und Massey°), von Boyd!?), Wenzl!!), sowie Kagan und Perel?2), 
Schulz und Brown?) und Kamke und Rose”) zeigen, daß durch die ver- 
einfachende Annahme der völligen Abschirmung des Sondenfeldes an der 
Grenze der Raumladungsschicht der Ausdruck für den Ionenstrom besonders 
stark verfälscht wird und damit natürlich auch der erhaltene Wert für die 
Jonentemperatur. Die Langmuirschen Vorstellungen werden durch die 
Annahme einer Übergangsschicht, in der noch ein schwaches Restfeld wirkt, 
korrigiert. Es zeigt sich aber, daß noch eine gewisse Unsicherheit im Verlauf 
der Konzentration der positiven Ladungsträger in der Schicht verbleibt. Die 
Vereinfachungen der Langmuirschen Vorstellungen haben nur geringen Ein- 
fluß auf die Werte der Elektronenkonzentration, und die Bestimmung der 
Elektronentemperatur bleibt davon unbeeinflußt. Daher begnügt man sich 
gewöhnlich damit. T_ aus dem Anstieg des Elektronenstromes im Anlauf- 
gebiet nach der vereinfachten Theorie zu bestimmen, und man benötigt dann 
den Ionenstromteil der Charakteristik nur zur Berechnung des Anteils des 
Elektronenstromes am Gesamtstrom im Anlaufgebiet. Mit Hilfe der Elek- 
tronentemperatur läßt sich aus dem Elektronensättigungsstrom die Elektronen- 
konzentration n_ berechnen. 

Bei hohen Entladungsstromdichten ist die Verwendung einer von Kagan 
und seinen Mitarbeitern?) aus der strengen Theorie der Sonde im Ionenstrom- 
gebiet entwickelten Näherung nützlich. Diese erlaubt, alle Größen bereits aus 
dem Ionenstromgebiet zu bestimmen. Die zugehörige Theorie ist besser fun- 
diert als die Sondentheorie von Langmuir, und man erhält wegen der geringen 
Sondenbelastung unverfälschtere MeBergebnisse. Der Einfluß von Magnet- 
feldern auf die Sondenmeßwerte wird nahezu eliminiert, weil dieser sich haupt- 
sächlich auf den Elektronenteil der Charakteristik auswirkt. 

Für den Sättigungsstrom erhält man den Ausdruck: 


i, = (1) 
v* ist die Geschwindigkeit der Ionen an der Schichtoberfläche. Sie errechnet 
sich aus der Beziehung: 


S, ist die Oberfläche dieser Schicht und M die Ionenmasse. Für die Schicht- 
dicke erhält man nach ') die Bestimmungsgleichung 


e t m 
‘| Mr 

°»), D. Bohm, E.H.S. Burhop u. H.S. W. Massey, Characteristics of Electrical 
Discharges in Magnetic Fields; Guthrie and Wakerling (MeGraw-Hill Book Comp., Inc., 
1949). 

1) R. Boyd, Proc. Roy. Soc. London (A) 201, 329 (1950). 

N) F. Wenzl, Z. angew. Physik 2, 59 (1950). 

22) W.M.Karan u B.M. Mepeab, DAH CCCP 95, 765 (1945). 

1) G.J. Schulz u. S.C. Brown, Physic. Rev. 98, 1642 (1955). 

14) D. Kamke u. H. J. Rose, Z. Physik 145, 83 (1956). 

15) Kapzow, Elektrische Vorgänge in Gasen und im Vakuum (Deutscher Verlag der 
Wiss. 1955). 
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In Formel (3) ist 6? eine tabulierte Funktion von = 15), r, der Schichtradius, 
r der Sondenradius und / die wirksame Sondenlänge. Die Elektronentemperatur 


kann man aus einer Formel bestimmen, die in einem kleinen Bereich um den 
Nullpunkt des Sondenstromes gilt 


di, 
din (=) 
Praktisch verfährt man bei der Rechnung so, daß man die graphisch nach 
U, differenzierten Werte des Sondenstromes um i, = 0 halblogarithmisch 
aufträgt, um aus dem Anstieg der Geraden die Elektronentemperatur zu er- 
halten. Aus der Gleichung: 


Us, =—T_In (4] (5) 


am, 


die sich aus der Gleichheit von Ionen- und Elektroneneinstrom für stromlose 
Sonde (?, = 0) ergibt, läßt sich die Spannung der Sonde, bezogen auf das 
Plasmapotential für diesen Punkt, finden. Diesen Wert benötigt man zur 
Festlegung des Plasmapotentials U,, bezogen auf das Kathodenpotential. Der 
Wert U, in Gleichung (3) bezieht sich ebenfalls auf das Plasmapotential. Für 
einen Punkt im Sättigungsstromgebiet wird die Gleichung (3) angewendet. 
Da :; und U, bekannt sind, erhält man 5° und bestimmt mittels der Tabellen 
schließlich S,. Mit dem aus Gleichung (2) berechneten Wert für v. kann man 
nun in Gl. (1) eingehen und n. errechnen. 

Ähnliche Berechnungen des Verhaltens der Sonde imIonenstromgebiet 9) 13) 14) 
führen in der Ionenstromdichte zu Unterschieden um konstante Faktoren. Im 
weiteren sollen diese Diskrepanzen durch Vergleich mit den gemessenen Werten 
untersucht werden. 


Messungen 


i 


Für die Entladungsströme 16, 32, 48 und 64 A werden nach dem in ®) be- 
schriebenen Verfahren Sondencharakteristiken aufgenommen. Die Lage der 
Sonde wird dabei in Abständen von 0,3 zu 0,3cm von der Säulenmitte aus 


Abb. 2. Ionenstromteil, unter den 
gleichen Bedingungen aufgenom- 
men wie bei Abb. 1 


Abb. 1. Sondencharakteristik fiir J = 64 A, 
p = 1/10 Torr, Argon, Sondenstand in der 
Rohrmitte 
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geändert, und für jede Lage werten wir ungefähr fünf Charakteristiken aus. 
Neben den Gesamtcharakteristiken mit oberem Knick werden gesondert die 
Ionenstromteile in einem für die Auswertung günstigeren Maßstab registriert. 
Für einen Entladungsstrom von 64 A wird in Abb. 1 und 2 ein Beispiel gezeigt. 
Die halblogarithmise he Auftragung des Elektronenstroms ermöglicht die Be- 
stimmung der Konzentrationsv erteilung und der Elektronentemperaturver- 
teilung über den Radius nach der Langmuir-Methode. Ferner wird das 
Plasmapotential bestimmt. Die gleichen Ionenstromkurven werten wir schließ- 
lich nach der Methode von Kagan aus. 


Auswertung der Messungen und Diskussion der Ergebnisse = 
a) Die Elektronentemperatur ae 


In Abb. 3—6 sind die nach der Methode von Langmuir berechneten Mittel- 
werte der Elektronentemperatur als Funktion des Abstandes von der Rohr- 
mitte aufgetragen. Zu jedem Entladungsstrom gehört eine Gerade. Abb. 7—10 
zeigen die nach der Methode von Kagan berechneten Werte für T_ bei den- 
selben Entladungsstromstärken. Die Übereinstimmung der nach den beiden 
Methoden erhaltenen Werte für 7- ist gut. Die Elektronentemperatur ist 
annähernd bis auf etwa 0,5 em Abstand von der Wand für die jeweilige Strom- 
stärke konstant. Aus einer theoretisch gewonnenen Beziehung von der Form '*) 


= 1.16 107: (cp RP (6) 


deren Voraussetzungen fiir unsere Messungen nicht mehr in allen Teilen er- 
füllt sind, läßt sich in ausreichender Übereinstimmung mit den experimentellen 
Ergebnissen eine Elektronentemperatur von 7_= 18000 °K abschätzen. 


b) Die Trägerkonzentration und das Wandpotential 
Nach einem in 1?) beschriebenen Näherungsverfahren werden die Meßwerte 
durch Kurven der Form n_(r) = n_(0)- Jo (t Fr) angenähert, wobei J, die 


Besselfunktion nullter Ordnung und f eine Konstante ist. Das Näherungs- 
verfahren wird bis zur zweiten Ordnung durchgeführt und liefert die in den 
Abb. 11—14 eingetragenen Kurven für die Elektronenkonzentration als Funk- 
tion des Rohrradius. Die Werte nach Langmuir sind durch Kreise, die nach 
Kagan durch Kreuze gekennzeichnet. Sie unterscheiden sich im Mittel um 
den Faktor 1,9. Abweichungen der Meßwerte von einer Besselverteilung 
sind nicht nachweisbar. Die Werte f, die man aus der Näherung erhält, ent- 
sprechen aber nicht der ersten Nullstelle der Besselfunktion, wie das in der 
Theorie von Schottky!®?) auf Grund von vereinfachten Randbedingungen 
gefordert wird. Sie liegen etwa zwischen 2,2 und 2,3 anstatt bei f = 2,4. Das 
entspricht der Tatsache, daß die Ladungsträgerkonzentration an der Rohr- 


1) A,v. Engel u. M. Steenbeck, Elektrische Gasentladungen, Bd.1 u. 2 (J. 
Springer-Verlag, Berlin 1934). 

17) F.Kohlrausch, Praktische Physik, Bd. 1 (Teubner-Verlag 1950). 

8) W. Schottky, Physik. Z. 25, 342, 365 (1924). 
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wand noch von Null verschiedene Werte annimmt. Die von Fetz!) angestell- halt 
ten Überlegungen betreffs der Randbedingungen der Diffusionssäule ermög- _ 
fir d 
stron 
x10" 
12 
6 
6 
? 
air) 
x19" 
J2 
26 
20 
20 30 Sstd 10 20 30 Sstd 
[mm] fan] 6 
Abb. 6. Abb. 10. 
Abb. 3—10. Verlauf der Elektronentemperatur in Anhängigkeit vom Sondenstand für 
verschiedene Entladungsströme. Sstd. 0 entspricht der Rohrmitte. Abb. 3—6. Werte nach 
Langmuir. Abb. 7—10. Werte nach Kagan 
lichen eine Abschätzung der Abweichung des Wertes f. Fetz geht bekanntlich 
davon aus, daß sich eine Wandaufladung ausgebildet und die Ionen etwa die | Abb. 
Strecke einer freien Weglänge im freien Fall zur Wand zurücklegen. Sie er- | laduı 
en 
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halten dabei eine gerichtete Geschwindigkeit v,, die sich aus der Beziehung 
= 

i, (R) = n,(R)-e-v, Fwana (7) 


für den Ionenstrom zur Wand errechnen läßt. Da man den Ionensättigungs- 
strom durch Sonden messen kann, ist v, leicht zu bestimmen. Aus der Bedin- 


Abb. 12. Abb. 14. 


Abb. 11—14. Die Trägerkonzentration als Funktion des Sondenstandes für die Ent- 
ladungsströme. Abb. 11: 16 A. Abb. 12: 32 A. Abb. 13: 48 A. Abb. 14: 64 A. Sstd. 0 
entspricht der Rohrmitte. + = Werte nach Kagan, ¢ Werte nach Langmuir 
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gung, daß der lonenstrom und der Elektronenstrom zur Wand gleich sein 
müssen, kann man dann die Größe des Wandpotentials finden, 7 

In am a): (8) 
Experimentell erhalt man das Wandpotential mittels einer stromlosen Platten- 
sonde an der Wand. Wir erhalten fiir unsere Zylindersonde unter gleichen 
Bedingungen einen Wert, der mindestens näherungsweise mit diesem überein- 
stimmen muß. In Tabelle 1 sind in Abhängigkeit vom Entladungsstrom die 
nach (7) und (8) berechneten Werte für v, und Uwana sowie die mit der Sonde 
gemessenen U,x, zusammengestellt. Dabei wird für v, und somit auch für 
Uwana noch berücksichtigt, daß wegen der Bedingung der Quasineutralität der 
positiven Säule auch die Werte von n, (R), bestimmt nach Kagan und nach 
Langmuir, sich etwa um den Faktor 1,9 unterscheiden. 


w| _ 


e ‚ 
k T_ U wand 


Tabelle 1 


v, nach Uwana nach Uwana nach 
a vg Langmuir Kagan Langmuir 
4) [em/see] [Cgs] [Cgs] gs] 
3,3 - 105 1,5 - 105 21 10 24 10-3 21 - 
32 3,4 - 10° Et <3 20 - 10-3 26 - 10-3 26 - 10? 
64 3,0 - 105 1,4 - 105 - 24-10-3 28 - 10? 


Da für die Berechnung des Wandpotentials U wana. keine sehr große Genauigkeit 
erwartet werden kann, muß man die Übereinstimmung zwischen Rechnung 
und Experiment als ausreichend ansehen. 

Nach Fetz kann man zur Berechnung der Ladungsträgerverteilung so 
vorgehen, daß man aus der Beziehung 


gk TR ( Mvp 


(9) 
zunächst die Ionentemperatur berechnet. Die Formel wird aus dem Nähe- 
rungsansatz 


1 
> | eEir)dr (10) 
R 


erhalten, wobei die Energieverluste durch Stöße vernachlässigt werden und 
> = 1 gesetzt wird. Nach der Rechnung ergeben sich Werte für 7* zwischen 
5 8 8 
2460 °K und 2480 °K. Geht man mit diesen Werten in die Beziehung 
F= Ai 8k 

Y9aMT, 
@) i g + 
ein, wobei F = A ist und die freie Weglänge der Ionen /; aus 16) entnom- 

1 

men wird, so kann man f berechnen und mit unseren, aus den Experimenten 


gefundenen Werten vergleichen. Wir erhalten Werte für £ zwischen 2,0 und 
29 


Bedenkt man, daß nur im Bereich F > 0,15 eine größere Genauigkeit für 
f zu erwarten ist, da bei kleineren Werten die Kurve J,(f) asymptotisch zur 
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Achse verläuft, so stimmen die Werte mit dem Experiment recht gut überein 
Es ist zu beachten, daß die nach Formel (9) theoretisch gewonnene Ionentempe- 
ratur unsicher ist. Das ist aber für die Abschätzung der Abweichung von f 
unbedeutend, da T. in Formel (11) sehr unempfindlich eingeht. 


ec) Die Elektronenbeweglichkeit 


Wir gehen von der Beziehung 
J=eNb_E, (12) 
aus, wobei J der Entladungsstrom ist (der Ionenanteil ist zu vernachlässigen) 
und N die Zahl der negativen Ladungsträger pro Längeneinheit der positiven 
Säule darstellt. b_ bezeichnet die Beweglichkeit. Da wir J und E, gemessen 
haben und N durch Integration der gemessenen Trägerverteilung über den 
Säulenquerschnitt erhalten wird, läßt sich b_ leicht berechnen. Es gilt: 
R 
= f n_(r) r dr mit n_(r) = n_(0) r). (13) 
6 
Unter Benutzung der fiir Besselfunktionen der Ordnung x gültigen Formel 


= Ile) + FR) (14) 


läßt sich eine partielle Integration durchführen und man erhält 


, 22R-n(0),, by, 
N= J,(f). (15) 
Für b_ folgt aus (12) 10 
b_ t- J (16) 


22 R?2eE,n_(0) J, (€) 
Für die Trägerverteilung legen 
wir unsere Naherungskurven 


N Ww HAVE! 


Jy r) zugrunde. Endlich wer- 


den die Werte von b_ noch auf 
0 °C und 1 Torr reduziert. Dazu 
muß die Gastemperatur der Ent- 
ladung bekannt sein. 

In einer Arbeit von Wasser- 
rab!%) wird eine exakte Berech- 


un 
r 


nung für die mittlere Neutral- | pe 
gastemperatur fiir den Fall an- 12 % 16 18 20 22 24 26 26 or] 
gegeben, daß über den Ent- x10 


ladungsquerschnitt gleiche Gas- Abb. 15. Reduzierte Elektronenbeweglichkeit in 
Abhängigkeit von der Elektronentemperatur 
Gichte herrscht. Wiese Bedingung T_, Die Darstellung zeigt die Kurve von ®) 
ist bei uns sicherlich nicht genau unter Berücksichtigung der Korrektur 2!) mit 
erfüllt. Eine Berechnung unter den von uns gemessenen Werten 
Zugrundelegung der von uns 

gemessenen Rohrwandtemperatur von 112 °C ergibt eine mittlere Gas- 
temperatur von T, = 550 °K. 


g 


19) Th. Wasserrab, Z. Physik, 129, 190 (1951). 
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Da die Totvolumina des Entladungsrohres verhältnismäßig klein sind, 
erfährt die Teilchenkonzentration in der Säule durch deren Temperatur nur 
eine kleine Erniedrigung. Für die mittlere Gastemperatur 7, = 550 °K und 
einen bei 293 °K gemessenen Fülldruck erhält man einen Reduktionsfaktor 
Po = 0,88 + prin. (17) 
Das Produkt b_ - py wird schließlich in das Diagramm von Wojaczek und 
Rademacher”) eingetragen. Wie in Abb. 15 gezeigt wird, stimmen die 
Werte trotz der etwa dreißigmal größeren Stromdichte in unseren Messungen 
gut überein. Eine Entscheidung darüber, ob die Werte nach der Rechnung 
von Kagan oder nach Langmuir richtiger sind, wollen wir weiter unten 
durch Fehlerbetrachtungen herbeizuführen versuchen. 
a 


d) Das Radialpotential und die Radialfeldstärke by 


. . . r . 

Für eine Diffusionssäule errechnet sich unter Vernachlässigung von D. 
und b,, und wenn man annimmt, es wäre kein Magnetfeld vorhanden, die 
Potentialdifferenz zwischen Säulenmitte und Säulenrand nach 1) zu 

k n_(R) 
Uy(0) — Ur(R) = T In 18 
pi (0) = — n_(0) (18) 

Einen Vergleich der errechneten mit den gemessenen Werten des Plasma- 

potentials für 16 und 64 A Entladungsstrom zeigt Tabelle 2. 


Tabelle 2 


[A] n_(O) [Upi(0) — Upi(R)] [Cgs]  [Up1(0) — Upi( Jexp [Ügs] 
16 Kagan 0,138 98-10 | 84-1073 
Langmuir, 0,182 | 8,4 103 
64 Kagan 0,116 10,5 - 103 11,8 - 10-3 
Langmuir 0,138 2,0:30° 


Die berechneten Werte stimmen mit den MeBwerten befriedigend überein. Die 
mittels beider Meßmethoden experimentell bestimmten Werte weichen so 
wenig voneinander ab, daß sie nicht getrennt aufgeführt werden müssen. 

Die gemessenen Kurven von U,,(r) fiir J = 16 und 64 A werden graphisch 
differenziert, um den Radialgradienten zu erhalten. In dem Falle, daß das 
Eigenmagnetfeld merklichen Einfluß auf die Säule hätte, müßten Abweichun- 
gen im Verlauf der Radialgradienten und der berechneten Werte vorliegen. 
Derartige Vorgänge werden von Thonemann und Cowhig?) in Quecksilber- 
dampfentladungen gezeigt. Wie aus Abb. 16 und 17 zu ersehen ist, stimmen die 
berechneten Werte mit den nach beiden Methoden gemessenen recht gut über- 
ein. 

Schätzt man die maximal durch das Eigenfeld auftretende Feldstärke H 


2J 
H = f 


2) K.Wojaczek u. K. Rademacher, Ann. Physik 18, 237 (1956). 
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so kann H höchstens für 64 A Entladungsstrom am Rand der Säule 3,2 Oe 
werden. Die Kraft, die auf ein Elektron in radialer Richtung ausgeübt wird, ist 


“we 


e 
K=8 io (20) 
mit vg als axial gerichteter Geschwindigkeit. Wir wollen die elektrische Feld- 
‘berochnen, die diesem H in der Wirkung aquivelont wäre. 
Heb_E, Te: 


10 


wobei wir v, = b_ E, einsetzen. Unter Benutzung der Werte b_ = 3,5- 107 


2 
[aver |: E, = 0,115 [Volt/em] und H = 3,2 [Oe] erhält man für einen Ent- 


ladungsstrom von 64 A 
Emag = 4,3 - 10-4 [Cgs]. (22) 


Der Betrag ist so klein, daß ein Nachweis der Einwirkung des Eigenfeldes aus 
unseren Messungen noch nicht zu führen ist. 


muir-Methode) . 
Ann. Physik. 7. Folge, Bd. 8 
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Abb.16. Die Werte des Radialgradienten Abb. 17. Die Werte des Radialgradienten 
für einen Entladungsstrom von 16 A. für einen Entladungsstrom von 64.A. 
+ = Be, + Emag nach Kagan berechnet; Die Bezeichnung der Punkte erfolgt wie 
A = Ee. + Emag nach Langmuir berech- in Abb. 16 
= Be nach Messungen (Kagan- 
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oe e) Das Verhalten der Sonde im Ionensättigungsstromgebiet 


Eu Be Auf Grund von unterschiedlichen Annahmen und Näherungen ergeben sich 
bei der Berechnung des Ionenstromes durch die Autoren 9)12)13)14) Beziehungen, 
die sich nur durch konstante Faktoren unterscheiden. Es folgt eine Gegen- 
überstellung der von den verschiedenen Autoren gefundenen Formeln für den 
Ionenstrom: 


i, =m, € Bohm, Burhop und Massey?) (23) 


i, = Schulz und Brown") 
i, =n, ,e Kagan und Perel), Kagan’) (26) 
i, =n e Verne” -S Kamke und Rose) (27) 


Für y, wird als günstigster Wert 1/2 gefunden, während für € die Größe 
2,718 gilt. In den Formeln (23) bis (27) soll unter S die Sondenoberfläche und 
unter S, die Grenzfläche der Raumladungsschicht verstanden werden. 

Da T,/3 T_ ~ 1/100ist, kann man diesen Ausdruck in (23) vernachlässigen. 
Dann zeigt sich sofort, daß (23), (24) und (25) annähernd übereinstimmen. Wir 
setzen hier näherungsweise S,/S~ 1,2. Für die Beziehungen (26) und (27) 
liegen wesentlich größere Abweichungen vor. 

Wir wollen die drei unterschiedlichen Beziehungen für die Größe des Ionen- 
stromes, ausgedrückt durch die Gleichungen (23)—(25), (26) und (27), experi- 
mentell überprüfen. Der Ionenstrom geht in die Gleichung (5) für das Poten- 
tial U,, ein. Da U,, aus der Sondencharakteristik entnommen werden kann, 
ist ein Vergleich mit den errechneten Werten möglich. Die Beziehung (5) wird 
für den Punkt i, = 0 durch Gleichsetzung des Elektronen- und Ionenstromes 
für Gleichung (26) erhalten. Der Unterschied zwischen S, und S kann hier 
vernachlässigt werden. Für das Potential U,, ergeben sich durch Verwendung 
der unterschiedlichen Ansätze für den Ionenstrom folgende Beziehungen: 


“pt 


II k 1 4M ite 
Ua, =—T- In (aa Bohm, Burhop und Massey 
| 
Ug. =—T_In Kamke und Rose. (29) 
e m 
Durch Umrechnung läßt sich leicht die Relation finden: 
UL = Un = (30) 


Tabelle 3 zeigt eine Gegenüberstellung der mit Hilfe von (5), (28) und (29) er- 
rechneten Werte für U,, mit den experimentell ermittelten Werten un “4 


it = 0,4-n,e - 8 
M (24) 

7 Du 
We 
nis 
dal 
sch 
(26 
ten 
| deı 
Au 
| zu 
Tr 
die 
q M: 
da 
lie 
tic 

M 
sc 
4 ar 

be 
be 
d 
E 
Vi 

4 


E. Hayess u. K. Rademacher : Sondenmessungen in der positiven Säule 171 


Tabelle 3 
Sondenstand | I II um IV ry 
ul, [V] ul [Vv] unt [Vv] 
6 6,2 7,4 8,7 6,2 
12 6,2 7,4 8,6 6,6 
18 6,1 7,3 8,6 6,4 ‘ 
24 6,0 7,3 8,5 5,6 ) 
30 5,9 7,1 8,3 6,1 
36 5,9 7,0 8,2 6,9 
40 5,5 6,6 7,7 6,0 


Durch die Fehlerbreite der Werte von + 6—8% ist es schwierig, zwischen den 
Werten U 5 und un zu entscheiden. Die Werte U4)’ können durch das Ergeb- 
nis der Messungen als nicht gerechtfertigt angesehen werden. 

Die Untersuchung der Raumladungsschicht um die Sonde hat ergeben, 
daß Elektronentemperatur und Trägerkonzentration nach der Langmuir- 
schen Methode richtig sind. Ein Vergleich der aus den Formeln (23)—(25), 
(26) und (27) erhaltenen Konzentrationen mit den nach Langmuir berechne- 
ten kann eventuell weitere Aufschlüsse zur experimentellen Bestätigung einer 
der angeführten Beziehungen geben. Man erhält: 

2 
ni = 2 nil = (31) 

Aus Abb. 11—14 ersieht man weiterhin, daß die Werte nl um den Faktor 1,9 
zu klein sind gegen die aus dem oberen Knick der Charakteristik erhaltenen 
Trägerkonzentrationen n_. Die Gleichung (31) liefert dasselbe Ergebnis wie 
die Betrachtung der U,,. Zwischen den Werten nach Bohm, Burhop und 

Massey und nach Kagan bleibt noch eine Entscheidung offen. 

In einer weiter unten durchgeführten Fehlerbetrachtung wird noch gezeigt, 
daß die von uns bestimmten Trägerkonzentrationen etwa um 20%, zu hoch 
liegen können. Trägt man dem Rechnung, so verhalten sich die Konzentra- 
tionen 


ni ın_:n!! = 1:1,5:3. (32) 


Man kann also auf Grund dieses Ergebnisses keine sichere Entscheidung zwi- 
schen den von Bohm, Burhop und Massey und den von Kagan und Mit- 
arbeitern erhaltenen Beziehungen treffen. 


«) Es können keine merklichen Fehler durch Schwankungen der Betriebs- 
bedingungen der Säule entstehen. Die Brummspannung des Gleichstromnetzes 
beträgt nur 1%, der Betriebsspannung des Rohres. Schichtungen treten bei 
den verwendeten hohen Stromstärken nicht auf. Die Schichtungsfreiheit der 
Entladung wird oszillographisch kontrolliert. 

ß) Die Sonden werden vor jeder Messung gereinigt. Eine Impulsbelastung 
der Sonden erfolgt kurzzeitig. Unter diesen Voraussetzungen ergeben die Vor- 


versuche eine gute Reproduzierbarkeit der Meßergebnisse. on fe ; 
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y) Ein groBer Teil des Sondendrahtes ist durch eine Quarzkapillare abe: d 
deckt, deren Oberfläche sich im Plasma elektrisch auflädt. Um diese Kapillare 5 
bildet sich eine Raumladungsschicht aus. Dadurch ergibt sich die Abschir- er 
mung eines kleinen Gebietes in der Nähe des Überganges vom Quarzüberzug > . 
zur freien Sondenflache fiir den Einfluß der Ladungsträger. Eine Abschätzung 2 
nach der Methode von Kagan ergibt für die Größe der abgeschirmten Fläche 
den Wert F = 220d (o = Sondenradius, d = Schichtdicke ~ 0,018 cm). Die lad 
wirksame Sondenoberfläche wird nach dieser Abschätzung um etwa 4%, ver- ~ 
kleinert. | 
6) Im Ionenstromgebiet kann es besonders bei Edelgasen zur Verfälschung u 
der Meßwerte durch Emission von Sekundärelektronen aus der Sondenober- = 
fläche kommen. Man kann annehmen, daß dieser Effekt bei unseren Messungen We 
vernachlässigt werden kann. Der Sondendraht besteht aus Wolfram, der Be 
infolge des Impulsbetriebes nicht heiß genug wird, um zu emittieren. Photo- 
emission durch Metastabile ließ sich ebenfalls nicht beobachten. Außerdem Ku 
ist der Emissionseffekt in Ar gering im Gegensatz z. B. zu Neon”). ne 
e) Spiwak und Reichrudel*) zeigen, daß bei Messungen mit Zylinder- 
sonden die Elektronensättigungsströme bei Anwesenheit eines magnetischen 
Querfeldes von 10 Oe merklich geändert werden. Das Eigenmagnetfeld unserer 
Säule beträgt aber höchstens 3,2 Oe und ist noch zu gering, um merkbaren wei 
Einfluß auf die gemessene Trägerkonzentration zu haben. Fel 
Zusammenfassend können wir feststellen, daß nur der Fehler, der durch die ter: 
Aufladung der Sondenisolierung entsteht, von Bedeutung ist. Die Verkleine- xin 
rung der wirksamen Sondenfläche um 4%, muß bei der Fehlerabschätzung der wei 
Trägerkonzentration mit berücksichtigt werden. dus 
REN und Wojaczek) untersuchen den störenden Einfluß der Da 
Sonden selbst bei Gradientenmessungen. Sie stellen fest, daß sich dieser Effekt Lit 
besonders bei Driicken kleiner als 1,5 Torr bemerkbar macht und nach der ber 
Zweisondenmethode erhaltene MeBwerte um 5—10%, zu groß gemessen werden. als 
Diese systematische Abweichung muß auch unsere Ergebnisse beeinflussen. it 
Die Streuung der Gradientenwerte ergibt sich dadurch, daB die Ablese- ow 
genauigkeit des Schnittpunktes der Sondencharakteristik mit der U-Achse etwa 
bis auf 7%, bestimmt ist. Andererseits wirkt sich die Einstellgenauigkeit der Re 
Entladungsstromstärke auch als Abweichung in E, aus. Schätzt man diese mi 
mit 2%, ab, so kann man für den Fehler von E, insgesamt etwa 9%, angeben ' 
und der Gradient kann héchstens um 10% zu groB gemessen sein. q 4 12° 
e) Fehler der Elektronentemperatur 
Hier muß man die beiden Auswertungsmethoden getrennt untersuchen. Up 
Wir betrachten zunächst die Werte, die wir mittels der Langmuir-Methode als 
erhalten haben. Dabei wird 7'_ aus dem Anstieg des logarithmisch aufgetra: nu 
genen Elektronenstroms bestimmt. Eine systematische Abweichung ware nur deı 
-21) K. Wojaczek, Monatsber. DAW (Berlin) 2,1, 12 (1960). os 
22) C. Kenty, Physic. Rev. 48, 181 (1933). ' 
%) T.B. Cnusak an 3.M. Penxpyaeup, Bull. Acad. Sci. l’Urss 479 (1938). die 
* K. Rademacher u. K. Wojaczek, Ann. Physik 17, 293 (1956). i 
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durch eine nicht völlig parallel zum Oszillographenschirm ausgerichtete Film- 
ebene zu erwarten. Durch eine derartige Abweichung würde die Charakteristik 
steiler und damit 7_ verkleinert werden. Nimmt man die Abweichung von der 
Parallelität zu höchstens 2° an, so läßt sich für T_ = 17000 °K errechnen, daß 
der Fehler kleiner als 1%, werden muß. 

Statistische Fehler sind erklärbar aus der Einstellgenauigkeit des Ent- 
ladungsstromes und der Auswertegenauigkeit der aufgenommenen Kurven. 
Man kann den Gesamtfehler der beiden Einflüsse zu etwa 8%, angeben. 

In der Kaganschen Temperaturbestimmung wird der abklingende Ionen- 
strom vernachlässigt. Da aber das Verhältnis von Ionenstrom zu Elektronen- 
strom sehr klein ist, hat der Fehler dieser Vernachlässigung bestimmt einen 
Wert kleiner als 0,3%. Für die Einstellung der Filmebene gilt Entsprechendes, 
wie für die Langmuir-Werte gesagt wird. 

Für den statistischen Fehler, der sich auch aus Auswertegenauigkeit der 
Kurven und Meßfehler im Entladungsstrom zusammensetzt, kann man 7%, 
rechnen. 


d) Fehler der Trägerkonzentration 


Wenn die „Charakteristiken nach der Langmuir-Methode ausgewertet 
werden, so können Fehler in den Größen i_, T_ und F, eine Rolle spielen. Die 
Fehler in T_ haben wir soeben untersucht. Da i_ für den Punkt der Charak- 
teristik bestimmt wird, wo die Sonde den Wert Up, annimmt, und die Appro- 
ximation des oberen Knickes der Kurve auf verschiedene Weise vorgenommen 
werden kann, besteht für i” eine gewisse Unsicherheit. Wir bestimmen Up, 
durch lineare Extrapolation des Elektronensättigungsstromes; die Autoren?) 
wenden eine nichtlineare Vorschrift an. Unsere Werte müßten gegen die aus 
der nichtlinearen Extrapolation erhaltenen deshalb um etwa 20%, größer sein. 
Da aber der Wert als Up anzusehen ist, in dem zuerst Abweichungen von der 
Linearität der logarithmisch aufgetragenen Kurve des Elektronenstroms zu 
bemerken sind, müßte die nichtlineare Extrapolation genauer sein. Wir müßten 
also i_ mindestens als 20%, zu groß annehmen. In i_ ergibt sich durch Ablese- 
genauigkeit und Auswertegenauigkeit noch eine statistische Schwankung von 
etwa 6%. 

Einem systematischen Fehler in F, wird schon in Abschnitt a) y) 
Rechnung getragen (etwa 4%,). Ein Meßfehler in der Sondenfläche ist höchstens 
mit 1,5%, anzugeben. 

In n_ erwarten wir nach dieser Betrachtung einen Gesamtfehler von 11 bis 
12%, wobei n_ noch um 17—20%, zu groß gemessen sein kann. 

Bei der Rechnung nach Kagan untersuchen wir wieder die Vernach- 
lässigung in den Formeln. In (5) wird In r,/r vernachlässigt. Da r,/r von der 
Größe 1,2 ist, ergibt sich in Uy, ein Fehler von höchstens 2%, und damit in 
Up, von 0,4%. Mit etwa 0,7%, geht in unserem MeBbereich die Abweichung 
als Fehler von n_ ein. Weitere Fehler können durch i,, T_ und S, in die Rech- 
nung kommen. Der Fehler in i, ist höchstens ein Ablesefehler, den man aus 
der Kurve zu etwa 3%, abschätzen kann. Für 7'_ ist die Untersuchung bereits 
durchgeführt. Es bleibt noch eine Abschätzung der Abweichungen in S,. 
Diese Fläche bestimmt sich aus ß?, wobei in die Rechnung das Potential Ug, 
die Sondenlänge und der Sondenradius eingehen. Der systematische Fehler 
in U, ist bereits berücksichtigt, so daß man einen statistischen Fehler in Ug 
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von ungetihr 3%, berücksichtigen muß. Für ! und r kann man je einen Fehler 
von 0,75%, erwarten. Daher würde der Fehler in ß? etwa 9%, betragen und folg- 
lich $, mit 1% Fehler behaftet sein. 

Insgesamt erhält man also einen Fehler von ungefähr 7—8% in n_ und eine 
systematische Abweichung von 1—2%. 

Diese Fehlerbetrachtung zeigt, daß der Unterschied in den nach den beiden 
Methoden berechneten Trägerkonzentrationen nicht allein aus der Fehler- 
breite zu erklären ist. Es gilt hier das bereits über die Kagansche Näherungs- 
formel Gesagte. 

e) Fehler der Elektronenbeweglichkeit 


Abweichungen in b_ py haben ihren Ursprung in Fehlern von J, n_ und E,. 
Ferner spielt auch der Fehler des Reduktionsfaktors p, eine Rolle. Die Fehler 
in n_ und E, haben wir bereits diskutiert. Ablesefehler des Entladungsstromes 
betragen etwa 6%. Für p, ist ein Fehler von etwa 8% zu erwarten. Nach der 
Methode von Langmuir folgt so ein Fehler in b_ py von 34%, und die Werte 
können bis zu 28%, zu klein sein. 

Nach Kagan ist ein Fehler zu erwarten von 31%, und die Werte können 
nur bis 12%, zu klein sein. 

Wir ersehen schon aus der Betrachtung, daß die Streubr@ite der Werte 
nach Kagan geringer ist. Noch besser ist das aus den Meßwerten selbst zu 
ersehen. 

Aus dieser Betrachtung lassen sich bezüglich der Anwendung der ver- 
schiedenen in der Anwendung üblichen Beweglichkeitsformeln die Ergebnisse 
von Wojaczek und Rademacher bestätigen. Weitere Schlüsse zu ziehen 
verbietet, wie auch bei diesen Autoren, die Meßungenauigkeit. 

Dem Leiter des Instituts, Herrn Prof. Dr. Rompe, möchten wir für sein 
förderndes Interesse an der Arbeit danken. Ebenfalls danken wir Herrn Dr. 
Wojaczek für fruchtbringende Diskussionen. 


Berlin, Deutsche Akademie der Wissenschaften, Physikalisch-Technisches 
Institut, Bereich Strahlungsquellen. 
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Ae Bei der Redaktion eingegangen am 23. November 1960. ze ee 
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Zur Theorie der Stoßwellenausbreitung und der 2 
Kanalentwicklung bei starken Funkenentladungen _ 


Von H. Morgenroth, H. Heß, K. Kischel, K. Rademacher und K. Seliger 


Inhaltsübersicht 


Die von §S. I. Drabkina!) angegebene Formel für die Geschwindigkeit 


der bei starken Funkenentladungen auftretenden Stoßwellen wird auf beliebige a 
Symmetrieverhältnisse erweitert. In Anlehnung an das Experiment wird der = 
ebene Fall ausführlich diskutiert. Durch geeignete Spezialisierung der in der = 
allgemeinen Formel auftretenden Energiefunktion läßt sich das experimentell - - 
gefundene Maximum in der Stoßwellengeschwindigkeit theoretisch beschreiben. i, 
Es folgen einige Bemerkungen zur Theorie der Kanalentwicklung. 7 


Starke Funkenentladungen und ihre Busiie wurden bisher 
im wesentlichen von zwei Arbeitsgruppen untersucht, und zwar von S.L. 
Mandelstam und seinen Mitarbeitern!-*) und von A. C. Kolb5). Die Ziel- 
setzung der durchgeführten Arbeiten ist unterschiedlich. Während in den 
älteren Arbeiten von Mandelstam und Mitarbeitern der Funkenkanal selbst, 
insbesondere seine schnelle Ausdehnung und die damit verbundene Ausbreitung 
von Stoßwellen im Mittelpunkt der Untersuchungen steht, wird in den jüngeren 
Arbeiten von Kolb eine Anwendung der bei der Funkenkanalentwicklung 
auftretenden Stoßwellen untersucht, die zur Erzeugung hoher Temperaturen 
führen soll. Kolb stützt sich dabei auf die Arbeiten von Fowler und Mit- 
arbeitern®). In diesen Arbeiten wird erstmalig die T-Rohr-Anordnung benutzt. 
Aus den beobachteten Lumineszenzerscheinungen in Form leuchtender 
Störungsfronten, die sich mit großer Geschwindigkeit vorwärtsbewegen und 


1) §.I. Drabkina, KITO 21, 473 (1951); (Dt. Übers.: Abhandlungen aus der sowj. 
Physik, Folge III, 295 "(1953)). 

2) LS. Abramson u. N. M. Gegetschkori, KITO 21, 484 (1951); (Dt. Übers.: Abh. 
sowj. Phys., Folge III, 311 (1953)). 
an 3) N. M. Gegetschkori, HIT O 21, 493 (1951); (Dt. Übers.: Abh. sowj.Phys., Folge III, 

(1953)). 

4) I. S. Abramson, N. M. Gegetschkori, S.I. Drabkina u. S.L. Mandelstam, 
17, 862 (1947). 

5) A.C. Kolb, Physic. Rev. 107, 354 (1957). 

6) R. J. Lee u. R. G. Fowler, Physic. Rev. 81, 457 (1951); Fowler, Goldstein u. 
Clotfelder, Physic. Rev. 82, 966 (1952); Fowler, Atkinson u. Marks, Physic. Rev. 
87, 966 (1952). 
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über Zeiten existieren, die die zu erwartenden Abklingzeiten bei weitem über- 
treffen, schließen Fowler und Mitarbeiter auf Stoßwellen, deren Geschwindig. 
keiten zu selbsttätiger Anregung und Ionisation hinter der Stoßwellenfront 
ausreichen. Um die Wirksamkeit der genannten Anordnung noch zu erhöhen, 
benutzt Kolb ein senkrecht zur Entladungsachse und zum Seitenarm gerichte- 
tes Magnetfeld und erzielt damit höhere Stoßwellengeschwindigkeiten. 

Einen wesentlichen Beitrag zur Theorie der Kanalentwicklung und der 
Stoßwellenausbreitung bei starken Funkenentladungen stellt die Arbeit von 
S. I. Drabkina!) dar. Eine weitere, leider unzugängliche theoretische Unter- 
suchung von E. G. Harris wird in der Arbeit von Kolb5) erwähnt. Dieser 
Untersuchung entnimmt Kolb einen ähnlichen Ausdruck für die Stoßwellen- 
geschwindigkeit im Falle ebener Symmetrie, wie er bereits in !) für den zylinder- 
symmetrischen Fall angegeben wird. Schließlich führt Braginski’) zum 
gleichen Problem theoretische Untersuchungen durch. Der Autor geht dabei 
jedoch von einem Kolbenmodell für die Kanalausbreitung aus, das wir nicht 
weiter verfolgen wollen. 

Die Arbeit von Drabkina stützt sich auf die umfangreichen theoretischen 
Untersuchungen, die Sedow®)?) zur Frage der Ausbreitung starker Explo- 
sionswellen durchgeführt hat. Danach können Lösungen des jeweiligen hydro- 
dynamischen Gleichungssystems auf der Basis der von Sedow entwickelten 
Ähnlichkeits- und Dimensionsmethoden!P) gewonnen werden. Ein ähnliches 
Verfahren findet man bei Courant und Friedrichs!!) unter der Bezeichnung 
Progressionswellentheorie und sehr ausführlich bei C. F. v. Weizsäcker!2) in 
Form der in mehreren Arbeiten entwickelten und diskutierten Homologie- 
lösungen. 

S.I. Drabkina behandelt die freie Funkenentladung als zylindersym- 
metrische Explosion einer Ladung E, und gibt die Ausbreitungsgeschwindig- 
keit der Stoßwellenfront als Funktion der Zeit, der Energie E, und der Dichte 
0, des umgebenden Mediums an. Darüberhinaus wird ein entsprechender Aus- 
druck für den Fall angegeben, daß die zugeführte Energie nicht momentan, d. h. 
im Zeitintervall At = 0, sondern allmählich freigesetzt wird, so daß die im 
Stoßwellenbereich enthaltene Energie von t abhängt. Diese Formel gilt exakt 
nur für den Fall, daß dem betrachteten gasdynamischen Vorgang in jedem 
Moment die Eigenschaft der ‚Selbstähnlichkeit‘‘ zugeschrieben werden kann. 
Darunter versteht man die Möglichkeit einer stationären Darstellung der den 
gasdynamischen Zustand kennzeichnenden Größen p/p,, u/u, und o/e, als 
Funktionen von r/R. (p,, u, und o, stellen den Druck, die Strömungsge- 
schwindigkeit bzw. die Dichte hinter der Front der Störungswelle dar.) Die 
Homologieeigenschaft des Vorganges wird hier nur benutzt, um die Ge- 
schwindigkeit der StoBwellenfront als Funktion von ¢ zu gewinnen. Auf eine 
vollständige Darstellung von p,/p usw. als Funktionen r/R wird in der vor- 
liegenden Arbeit verzichtet. 

8.1. Braginski, 34, 1548 (1958). 
upnk.s. Matem. u Mex.(Z. angew. Math.u. Mech.)10, 241 (1946). 
9%) L. I. Sedow, JIAH CCCP 47, 94 (1945). 


2) L. I. Sedow, Dimensions- und ähnlichkeitstheoretische Methoden in der Mechanik, 


UdSSR, 1944. 

u) R. Courant u. K.O. Friedrichs, Supersonic flow and shock waves, New York 
1948. 

12) Zusammenfassende Arbeit: S. v. Hoerner, Fortschr. d. Phys. 6, 375 (1958). 
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Wir werden im Sinne einer Zusammenfassung die Ableitung des in Rede 
stehenden Ausdrucks im folgenden Abschnitt in etwas allgemeinerer Form 
gleichzeitig für alle drei Symmetriefälle wiederholen. Dabei gehen wir von einer 
genauen Formulierung des Problems der Stoßwellenausbreitung als Folge 
starker Funkenentladungen aus, wobei wir die Funkenentladung als zeitlich 
aufgelöste Explosion einer vorgegebenen Energie E, auffassen. Die Frage nach 
dem Grad der Verletzung der vorausgesetzten Homologie der Strömungsform 
während die Energie allmählich zugeführt wird, konnte nicht beantwortet 
werden. Der Erfolg rechtfertigt zunächst die Voraussetzung. Die mit bestimm- 
ten Energieverlusten verbundenen Anregungs- und lIonisationsprozesse 
hinter der Stoßwelle sollen hier nicht berücksichtigt werden. 

Die gewonnene Formel für die Stoßwellengeschwindigkeit D(t) bei zeitlich 
ausgedehnter Energiezufuhr H(t) wird dann im nächsten Abschnitt der theo- 
retischen Erklärung des bei unseren Versuchen beobachteten Geschwindig- 
keitsprofils D(t) zugrunde gelegt. Dazu wird E(t) an Hand experimentell ge- 
wonnener Anhaltspunkte für den Verlauf der Leistungskurve spezialisiert. 
Charakteristisch für das erwähnte Geschwindigkeitsprofil ist ein relatives 
Maximum in einem Abstand von einigen cm von der Entladungsachse. Ein 
solches Maximum beobachtet auch Kolb5), ohne jedoch eine plausible Er- 
klärung zu geben. Aus der Annahme, daß die in den Stoßwellenbereich in- 
vestierte Energie in bestimmter Weise von der Zeit abhängt, läßt sich der beob- 
achtete Effekt zwanglos erklären. 

Der 4. Abschnitt der vorliegenden Arbeit beschäftigt sich mit einer be- 
stimmten, von der bisherigen abweichenden Vorstellung über die Kanalaus- 
breitung unmittelbar nach dem Durchschlag. Dieser Bereich wird von der den 
Abschnitten 2 und 3 zugrunde liegenden Theorie nicht erfaßt, welche u. a. davon 
absieht, daß schon der Durchschlagsfaden eine endliche Ausdehnung besitzt. 
In der vorgeschlagenen Behandlungsweise wird dagegen derDurchschlagsfaden 
als die Hochdruckkammer eines ,,zylindersymmetrischen Stoßwellengefäßes‘“ 
angesehen. Die konsequente Durchführung erfolgt dann allerdings nur für 
den Fall des ‚ebenen‘ (klassischen) StoBwellenrohres, da hierfür bereits eine 
vollständige Theorie vorliegt (siehe z. B. ")). Dies ändert jedoch nichts an den 
sich ergebenden grundsätzlichen Vorstellungen über den Mechanismus von 
Stoßwellen- und Kanalausbreitung im Anfangsstadium. 


2. Aufgabenstellung und Ableitung der Geschwindigkeitsformel 


Gegenstand der folgenden Betrachtungen ist die auf einen Punkt, eine 
Gerade oder eine Ebene konzentrierte Explosion einer bestimmten Energie 
E,, die in der Form 


freigesetzt wird. Eine solche Explosion erzeugt eine Störung, deren Front sich 
je nach den vorliegenden Symmetrieverhältnissen kugelsymmetrisch, zylinder- 
symmetrisch oder eben ausbreitet. Zwischen Front und Explosionszentrum 
befindet sich ein Strömungsfeld {u (r, t), p(r, t), o(r, t)}, das der Kontinuitats- 
gleichung, der Bewegungsgleichung und der Adiabatengleichung unterworfen 
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ist. An der Front selbst sollen die Ubergangsbedingungen fiir eine starke 
StoBwelle (2, = 0) erfüllt sein: 


| 


y = ¢,/¢,. 


Die bis zur Zeit t(<t,) dem Energievorrat E, entnommene Energie soll sich 
in Form von thermischer und kinetischer Energie restlos im Strémungsfeld 


wiederfinden : 
t RW) 
R 
Für t > t, bleibt die Gesamtenergie im Strömungsfeld konstant. [ f(r) dr = 
0 


f V (R) ist das Volumen des Strömungsfeldes im kugelsymmetrischen, zylinder- 
7 symmetrischen oder ebenen Fall. 
; Von diesen grundsätzlichen Annahmen ausgehend läßt sich die Aufgaben- 
= stellung folgendermaBen formulieren: 
u Gegeben sind o, sowie dE/dt als Funktion von ¢ für das Zeitintervall 
0<t< &. Es gelten die Gleichungen (2.1) bis (2.3). Gesucht ist eine mit (2.1) 
bis (2.3) verträgliche Lösung der gasdynamischen Grundgleichungen sowie eine 
Beziehung R = R(t) als Bahngleichung für die Stoßwellenfront und schließlich 
D = dR/dt für deren Geschwindigkeit. 

Unter der Voraussetzung 


E(t)=K#, E(t)=E, (2.4) 


findet man bei Courant-Friedrichs!) eine vollständige Lösung des Problems 
auf der Basis der sogenannten Progressionswellentheorie. Die Energiefunk- 
tionen der Form (2.4) schränken den Anwendungsbereich dieser Theorie 
jedoch stark ein. Esergeben sich nur streng monotone, genauer, {(P) proportio- 
u nale Stoßwellengeschwindigkeiten (f(ß) ist eine rationale Konstante). Um 
° auch detailliertere Vorgänge erfassen zu können, wenden wir uns daher un- 
a mittelbar dem allgemeineren Standpunkt von Drabkina zu. Wir verzichten 
a dabei auf eine Darstellung des Strömungszustandes hinter der Front der 
u Störungswelle und beschäftigen uns nur mit der Geschwindigkeit der Wellen- 
front selbst. 

Ausgangspunkt der folgenden Überlegungen ist das Energieintegral (2.3). 
Wir schreiben den Integranden abkürzend in der Form 


(2.5) 
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Ist die Strömung homolog, so ist die Darstellung 
e(r) _ r 
stationär (Abb. 1). Führt man ER an Stelle von reinund setzt nach (2.7) 


e(ER) = g(&) e(R) Er) 
so ist E(R) 


1 
E(t) = J FER) Ras 


1 (2.8) 
e(R) a, £918) En dé, 
da f(r) ~ r„ ist (n = 0,1,2 entsprechend 
dem ebenen, dem zylindersymmetrischen 
bzw. dem kugelsymmetrischen Fall). An- Bw 
dererseits existiert ein Zwischenwert r’ derart, daß a” Tah 


R 
E(t) = e(r’) f f(r)dr =e V =eb, R™, m = 1, 2, 3. 
6 


Hierbei ist ¢ die mittlere Energiedichte und V das Volumen des Störungsfeldes 
zur Zeit ¢. Wir vergleichen (2.8) und (2.9) und erhalten 


1 
Rm = e(R)a, Rt (2.10) 
oder 
é = Ae(R), A= const, (2.11) 


da das Integral auf der rechten Seite von (2.10) zeitunabhangig und m = n + 1 

ist. Diese Darstellung gilt für jedes t. Also ist zu jeder Zeit ¢ die mittlere 

Energiedichte der Energiedichte am Rand des Störfeldes proportional. Wir 
weisen noch auf 

e(R)— e(R)— = 


hin. Danach ist &’ = const bzw. r’(t) = R(t) oder 7’ ~ R. Dies bedeutet, 
daß bei homologen Strömungsformen die Geschwindigkeit des Punktes, 
in dem wertmäßig die mittlere Energiedichte angenommen wird, der Ge- 
schwindigkeit der Randfläche der Störung proportional ist. Im einzelnen 
gilt dies für den mittleren Druck und für die mittlere kinetische Energie. 

Gleichung (2.11) stellt die wichtigste Beziehung der durchzuführenden 
Ableitung dar. Sie ermöglicht es, von (2.9) ausgehend mit Hilfe der Rand- 
bedingungen (2.2) wie folgt weiterzurechnen: 

E(t) =Ae(R)V =AV [Au =AV 

Wir setzen wieder V = b,, R™, m = 1, ¥ 3, lösen nach D = dR/dt auf, multi- | 
plizieren mit dt und mapa auf beiden Seiten: 
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> Daraus folgt die Bahngleichung des Stoßwellenkopfes: Pre 

ibn 


Die Differentiation nach t führt unmittelbar zu der gesuchten Formel für die 
Stoßwellengeschwindigkeit: 


m+2\ 4 


D(t) = 3* (m, y) (bn Qo) = al (2.13) 


Dieser Ausdruck ist für den zylindersymmetrischen Fall (m = 2, b,, = xh) 
mit dem von Drabkina angegebenen identisch. Wir werden uns im folgenden 
Abschnitt unserer experimentellen Anordnung entsprechend mit dem Fall 
ebener Symmetrie mit m = 1 und b,, = q beschäftigen, wobei g der Quer- 
schnitt des Rohres ist, in dem sich die StoBwelle reameemanel Zuvor jeden 
noch einige Hinweise. 

Wir spezialisieren die Energiefunktion auf den Fall 


E(t) = E, = const für alle t. 
Dann folgt aus (2.13) 


D(t) = A*(m, y) nas 


Diese Formel geht auch aus der Progressionswellentheorie nach Courant 
und Friedrichs für den Spezialfall 6 = 0 (siehe (2.4)) hervor. Für m = 1 
wurde sie von Harris berechnet und von Kolb in seiner Arbeit®) erwähnt. 
Sie gibt die Verhältnisse bei momentanem Umsatz der zur Verfügung stehen- 
den Energie E, wieder. Zugleich stellt (2.14) einen asymptotischen Ausdruck 
für die Stoßwellengeschwindigkeit in einem Zeitintervall tf) <t < tyrit dar. 
(txrit ist die Zeit, oberhalb der die Stoßwelle nicht mehr als stark angesehen 
werden kann und die Voraussetzung p, > po verletzt ist.) Dies geht aus folgen- 
der Überlegung hervor: Gemäß unserer Aufgabenstellung wird dem Störfeld 
bis zum Zeitpunkt ¢, Energie zugeführt. Von t, an aber ist H(t) = E, = const, 
und nach einiger Zeit stellt sich der gleiche Strömungszustand ein, der sich 
eingestellt hätte, wenn von vornherein der Fall momentanen Energieumsatzes 
vorgelegen hätte. Nach dieser Zeit verhält sich die Stoßwelle :ngeschwindigkeit 
also — (2.14). Wir werden dies im nächsten Abschnitt noch präzisieren. 


3. Spezialisierung der Energiefunktion E(t) 


_ Zur Erklärung des von uns beobachteten Maximums der StoBwellenge- 
schwindigkeit (Abb. 2) wollen wir E(t) spezialisieren. Einen Anhaltspunkt liefert 


dazu der von uns gemessene Leistungsverlauf L(t) = = bis zum Ende der 
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ersten Stromhalbwelle. Da wir eine geschlossene mathematische Darstellung 
versuchen, approximieren wir L(t) durch eine Parabel der Form er 


Lit) = At? 0<t<2r, (3.1) 


wobei t die Lage des Maximums von L(t) angibt. Durch die gewählte Naherung 
gelingt einerseits die erforderliche Integration; andererseits vermag die damit 


+ 


nar Cmax to”? 
Abb.2. Stoßwelle windig- Abb.3. Idealisierter Ver- 
keit im 7-Rohr. lauf der Leistungskurve. 


verbundene Idealisierung an den grundsätzlichen Verhältnissen nichts zu 
ändern. Wir werden am Schluß dieses Abschnittes dazu noch eine Bemerkung 
machen. 

Die Integration der Leistungskurve ergibt 


t 
E()'= At(t-5) .., 


Die in jedem Moment im Stoßwellenfeld enthaltene Energie Bra somit 


der Darstellung (siehe Abb. 4) ; 


AR (7-5) für 0<t< (8.31) 
E(t) = 
A g(r = E, für t > ty. (3.32) 


Der Zeitpunkt t = t,, indem die Energiezufuhr aussetzt, muß nicht unbedingt 


mit dem Zeitpunkt übereinstimmen, in dem die Ableitung der Energiefunktion 
verschwindet, d. h., die Energiezufuhr kann 1 E(t) 
schon früher aussetzen (etwa bei to Abb. 4). |; 
Allerdings wird die StoBwellengeschwindig- |’ 
keit dann mathematisch gesehen unstetig. 
Die Darstellung der Stoßwellengeschwin- 
digkeit gelingt in geschlossener Form mit 
Hilfe von (3.3) natürlich nur für das 
Intervall 0<t<.t, Für größere Zeiten 
muß zunächst numerisch gerechnet werden, ru 
und zwar so lange, bis die Stoßwellenge- TE ato 
schwindigkeit in Abhängigkeit von der Zeit = 
hinreichend genau der Formel (2.14) ge- 4 
Die allgemeine Beziehung (2.13) fiir die 
einer Energiefunktion E(t) zugeordneten 
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StoBwellengeschwindigkeit lautet fiir den Fall ebener Symmetrie, (die wir 
gemäß den geometrischen Verhältnissen unserer experimentellen Anord. 


nung wählen), 
-1/3 


Di) =( —) VEO al (3.4) 


Dabei wants A*(m, y) = & gesetzt; b,, ist gleich dem Rohrquerschnitt g. Das 


fain fat = Ait. +27) - =) ] 


3 
Damit ist 
= A [3 t 


= Funktion besitzt ein Maximum und gibt bei geeigneter Wahl der Kon- 
_ stanten den in Abb. 2 skizzierten Verlauf der StoBwellengeschwindigkeit bis 

zur Zeit {, prinzipiell richtig wieder. 
Für £>t, muß zunächst auf (3.4) zurückgegangen werden. Beachtet 


a ER (3.32), so ergibt sich a 
.. (t) le > 1/3° = (3.7) 
Vi VE(t) dt + | 


a af J VE dt lautet für den Spezialfall (3.31) 
6 


Wir schreiben 


ty 
J VE (dt = 0< <1, 
0 
und erhalten 
Ey’? 
[¢ —t)(1 — 


; Für große ¢ geht dieser Ausdruck in die Form (2.14) für m = 1 über. Be- 
an merkenswert ist dabei, daß nicht £ >‘), sondern nur t>t,(1—n) zu sein 
braucht, um mit der Formel (2.14) fiir den momentanen Energieumsatz 
weiterrechnen zu können. Für die spezielle Energiefunktion (3.21) ist = 
0,64 » 1/2, wenn ft, = 2r gesetzt wird. Somit gilt (2.14) für t >t. 

Die Gesamtdarstellung des in Abb. 2 skizzierten Verlaufes der Stoß- 
wellengeschwindigkeit wird somit durch die drei Beziehungen (3.6), (3.8) 
und (2.14) (für m = 1) entsprechend 0<t< t, t> t, und ¢t 4/2 wieder- 
gegeben. 

Wir interessieren uns noch für die Lage des Maximums in Abhängigkeit 
von der Zeitkonstanten rt. Dazu hätten wir (3.6) zu differenzieren. Zur Er- 
leichterung differenzieren wir (3.4) und erhalten 
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als Bedingung dafür, daß D verschwindet. Mit (3.2) und (3.5) ergibt dies 


Von hier führt der Ansatz 


zu der Gleichung 


9(2 — x) — 2(9 — 22) (1— z/3)P = 0, 


deren Wurzel mit ausreichender Genauigkeit bei x = 4/3 liegt. Nach (3.10) 
ist damit 

tmax 


wobei r nach (3. 1) der Zeitpunkt ist, in dem die m Leistungskurve 
ihr Maximum durchläuft. Zwischen dem Leistungsmaximum und dem Maxi- 
mum der Stoßwellengeschwindigkeit liegt also eine gewisse Zeitdifferenz mit 
konstantem Verhältnis der entsprechenden Zeitwerte. 


Für > 0 geht D—0. Der Beweis gelingt am besten, wenn man den 
Zähler sowie den Ausdruck (7 — =)" im Nenner von (3.6) nach Taylor ent- 
wickelt und nach der ersten Ableitung abbricht. Es zeigt sich dann, daB 
D(t) mit t wie {2/3 verschwindet. Die l’Hospitalsche Regel versagt hier. 

Wir interessieren uns schließlich noch für das Verhalten von (3.6) für r — 0. 
Dabei schrumpft der variable Bereich der Energiefunktion zwischen t = 0 
und ¢ = é auf Null zusammen und der Endwert E, der zugeführten Energie 


wird immer früher erreicht. Wir verlangen nun 
El.-o = E, für t> 0, (3.12) 


d.h., wir diskutieren die Annäherung an den Fall des momentanen Energie- 
umsatzes. 

Um (3.12) zu gewährleisten, muß die Konstante A in geeigneter Weise 
von t als Parameter abhängen. Nach (3.3) gilt für {;= 2r 


5 4 
=at”. 
Damit hängt E fiir ¢ > 2r nicht mehr von r ab: ai 
und wir kénnen 
a = 3/4 E, 
setzen. Fiir die Leistungskurve hat dies Late 


=ar’t=art 


zur Folge. Das ee muß also mit t—> 0 wie 1/r ansteigen. 
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Aus (3.11) folgt nun unmittelbar, daß sich die Lage des Maximums 
mit t—> 0 zum Zeitnullpunkt hin verschiebt. Das Verhalten der Höhe des 
Maximums bestimmen wir an Hand von (3.6). Es ergibt sich 


1/3 

| 

3 


Ein Vergleich mit (2.15) für den Fall m = 1 lehrt, daß das Maximum der Stoß- 
wellengeschwindigkeit für r— 0 mit der gleichen Potenz ansteigt, wie im 
Falle des momentanen Energieumsatzes für t > 0 die StoBwellengeschwindig- 
keit selbst. Diese Verhältnisse sind in 

| Dit) = der Abb. 5 dargestellt. 
Zum Schluß dieses Abschnittes soll 
Br noch darauf hingewiesen werden, daß die 
Beziehung (3.11) mitunter scheinbar ver- 
letzt ist, nämlich dann, wenn (experi- 
mentell) das Leistungsmaximum zeitlich 
später auftritt, als das Maximum der 
ifs StoBwellengeschwindigkeit. In diesem 
, Falle kann dasder ersten Stromhalbwelle 
7 entsprechende Hauptmaximum der Lei- 
7  stungskurve nicht für das Geschwindig- 
Abb. 5. D= f(t) für verschiedene keitsmaximum der Stoßwelle verantwort- 
t-Werte ,<r,<t, lich gemacht werden. Hier gibt eine Arbeit 
von B. Kühn?) einen Hinweis. Danach 
steigt die Leistungskurve bei höheren Drucken nicht monoton an, sondern 
durchläuft ein ,, Vormaximum“, das nach Kühn mit dem Zeitpunkt maximaler 
Energiedichte im Funkenkanal zusammenfällt. Diese maximale Energiedichte 
wird dem Durchschlagsfaden unmittelbar vor dessen Ausdehnung zugeschrieben. 
Sie entspricht der höchstmöglichen Stromdichte bei maximalem Ionisations- 
grad. Von diesem Zeitpunkt an wird der Widerstand des Funkenkanals nur 
noch durch seinen Querschnitt bestimmt. Er verringert sich durch die ein- 
setzende Ausdehnung so schnell, daß der Strom wegen der zu überwindenden 
Induktivität nicht im gleichem Maße ansteigen kann. Dadurch fällt die 
Leistung kurzzeitig ab und durchläuft ein Minimum. Der Zeitpunkt dieses 
Minimums entspricht etwa dem Übergang in das Bogenstadium, in dem 
die Stromstärke wieder leistungsbestimmend wird. Die experimentelle Beobach- 
tung dieses Vormaximums bedarf einer außerordentlich leistungsfähigen und 

hochauflösenden Meßanordnung. 


Aus den angeführten Überlegungen geht hervor, daß das Hauptmaximum 
in der Leistungskurve nur dann für das Maximum der Stoßwellengeschwindig- 
keit verantwortlich gemacht werden kann, wenn es mit dem Energiedichte- 
maximum zusammenfällt. Andernfalls kann über r ohne genauere Messungen 
keine Aussage gemacht werden. Die Beziehung tmax = 4/3 r bleibt natürlich 
gültig, da auch das Vormaximum eines solchen Leistungsverlaufs durch eine 
Parabel approximiert werden kann. 


18) B.Kühn, Ann. Physik 8, 241 
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In der Arbeit von H. Heß u. a.!4) werden die hier dargestellten Ergebnisse 
im einzelnen mit dem Experiment verglichen. Dabei wird auch der bei uns 
vorliegende Zusammenhang zwischen dem Leistungsverlauf und dem Maxi- 
mum der Stoßwellengeschwindigkeit diskutiert. 


zn 


4. Zur Frage des Kanalmechanismus 

Die in den Abschnitten 2 und 3 entwickelten Vorstellungen können offenbar 
nur für nicht zu kleine Zeiten t gelten. Dies geht schon daraus hervor, daß 
D=0 als Anfangswert der Stoßwellengeschwindigkeit nicht vorkommen 
kann, da die untere Grenze für die Geschwindigkeit einer Stoßwelle die Schall- 
geschwindigkeit c, des Mediums vor der Stoßwelle ist. Im vorliegenden Falle 
wurde darüberhinaus mit der Näherung für starke StoBwellen gerechnet,d. h. 
mit M = D/e, >1. Außerdem wurde in der Näherung der Abschnitte 2 und 3 
die räumliche Ausdehnung des Durchschlagsfadens vernachlässigt. Die spe- 
ziellen Vorgänge im Durchschlagsfaden sowie die Herausbildung des eigent- 
lichen Funkenkanals können daher in dieser Näherung nicht erfaßt werden. 
Diese Überlegungen zwingen dazu, einer Betrachtung des Anfangsstadiums 
der Kanal- und Stoßwellenausbreitung realere Annahmen zugrunde zu legen. 

Dazu gehen wir von der Vorstellung aus, daß bei einer Funkenentladung 
unmittelbar nach dem Durchschlag ein fadenförmiger Bereich sehr hoher 
Stromdichte und damit sehr hoher innerer Energie entsteht. Der Durch- 
schlagsfaden, in dem gegenüber dem angrenzenden Medium eine hohe Tempe- 
ratur und somit ein hoher Druck herrschen, gibt Anlaß zu einem Ausgleichs- 
vorgang, der das Auftreten einer Stoßwelle zur Folge hat. Betrachtet man nur 
sehr kleine Zeiten, so läßt sich dieser Ausgleichsvorgang mit dem Vorgang in 
einem klassischen Stoßwellenrohr vergleichen. Dabei ist der Durchschlags- 
faden mit der Hochdruckkammer, das an den Durchschlagsfaden angrenzende 
Medium mit der Niederdruckkammer eines solchen Rohres zu identifizieren. 
Bei diesem Vergleich muß man jedoch zwei wesentliche Unterschiede im Auge 
behalten: 

1. Im klassischen Stoßwellenrohr ist die Temperatur der Hochdruckkam- 
mer gleich der Temperatur der Niederdruckkammer, und der Druckunter- 
schied in den beiden Kammern vor dem Durchschlagen der Membran wird 
durch eine Dichteänderung in der einen Kammer erzeugt. Im Gegensatz 
dazu wird der hohe Druck im Durchschlagsfaden durch Temperaturerhöhung 
hervorgerufen, während die Dichte gleich der Dichte des umgebenden Mediums 
ist. Dieser Umstand schließt jedoch die Anwendung des beim Stoßwellenrohr 
üblichen Rechenverfahrens noch nicht aus. 

2. Im Gegensatz zur ebenen Symmetrie der Vorgänge im Stoßwellenrohr 
ist der Ausgleichsvorgang für den Durchschlagsfaden zylindersymmetrisch. 
An Stelle des klassischen Stoßwellenrohres müßte also folgendes Problem 
untersucht werden: Eine Membran in Form eines Hohlzylinders trennt eine 
bestimmte Gasmenge hohen Druckes und hoher Temperatur von dem um- 
gebenden Medium niedrigen Druckes, niedriger Temperatur und gleicher Dichte. 
Wie vollzieht sich der Ausgleichsvorgang, wenn die Membran plötzlich und 
überall gleichzeitig entfernt wird ? Offenbar werden die Geschwindigkeiten 


4) H.Heß, K. Kischel, H. Morgenroth, K. Rademacher u. K. Seliger, Ann. 
Physik. 
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der Stoßwellenfront, der Kontaktfläche und der Front der auf die Achse des 
Hohlzylinders zulaufenden Verdünnungswelle von der Zeit abhängen. Damit 
entfällt eine zum ebenen Fall des Stoß- 

| | % % wellenrohres analoge Beschreibung, 


at en 6a denn hier beruht das Lésungsverfahren 
VAX\e gerade darauf, daB sich fiir die Fort- 
\ K 


pflanzung der Verdiinnungswelle, der 


= 6b Kontaktfläche sowie der StoBwellen- 


front konstante Geschwindigkeiten ein- 
Jo stellen. 


* Das _,,zylindersymmetrische Stoß- 
wellenrohr‘“‘ ist mathematisch noch 
nicht untersucht worden. Einige grund- 

geg 
sätzliche Anhaltspunkte zur Klärung 
der Vorgänge bei der Ausdehnung des 
Funkenkanals lassen sich jedoch schon 
am Modell des klassischen ebenen 
Stoßwellenrohres (Abb. 6a) gewinnen. 
tine ausführliche Theorie des ebe- 
6e nen Stoßwellenrohres findet man bei 
ut Courant-Friedrichs!!). Danach las- 
rl | sen sich aus den vorgegebenen An- 
x oF fangszustinden 0 = {pp, 09, Uy = 0)} 
Abb. 6 und 2 = {p,, 0, Ug = 0)} mit py > pp, 
a) Stoßwellenrohr; Rn 09 > 0,, alle Größen berechnen, die den 

tle: Charakteristik dr/dt gasdynamischen Zustand des Systems 
welle; 1arakteristik : AX r 

nach dem Zerstören der Membran be- 
taktfläche K: dz/dt = u,, Stoßwelle schreiben. Die Ergebnisse dieser Rech- 
S:dx/dt =D. nungen sind schematisch in den 
c) Dichteverteilung zur Zeitt=t. Die Abbn. 6b bis 6f wiedergegeben. Die Be- 
gestrichelte Linie deutet das rasche Zer- wegungsvorgänge verlaufen durchweg 
fließen des Dichtesprunges an der Kon- 5 
Pa am mit konstanter Geschwindigkeit. Nach 
d) Temperaturverteilung zur Zeit t= t,. links in die Hochdruckkammer hinein 
e) Druckverteilung zur Zeit t= t,. Kein pflanzt sich eine sogenannte Verdün- 
„ Drucksprung an der Kontaktflache.  nungswelle fort. Der Kopf dieser Ver- 
f) Verteilung der Strömungsgeschwindig- hü 
keit zur Zeit t = A Kein Sprung an der “ünnungsweie bewegt sich aul einer 


u x 


Kontaktfläche. Charakteristik 02), deren Neigung 
— = _ = — —( 
dt Us V’ 05 ’Ug ( = 


konstant ist. Entsprechend bewegt sich der RR der Verdünnungswelle 
mit der Ceschwindigkeit 

dx Pı 
dt = uy Cy = uU | 

wobei u, die Strömungsgeschwindigkeit und c, die Schallgeschwindigke it in 
dem zu ermittelnden Zwischenbereich 1 angeben. Nach rechts in die Nieder- 
druckkammer hinein pflanzt sich mit der Geschwindigkeit D = const der 
Kopf einer ‚„‚konstanten‘‘ StoBwelle fort, nach deren Passieren das Gas die 
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konstante Strömungsgeschwindigkeit vu, annimmt. Mit eben dieser Geschwin- 
digkeit bewegt sich auch eine ‚‚Kontaktunstetigkeitsfläche‘, die das Gas aus 
der ehemaligen Hochdruckkammer von dem der Niederdruckkammer trennt. 
An dieser Fläche erleidet die Dichte einen Sprung. Der Druck bleibt jedoch 


stetig; wegen T = wird mit der Dichte auch die Temperatur unstetig. 


p 
oR 
Die Abbn. 6c, 6d und 6e geben den Dichte-, den Temperatur- und den Druck- 
verlauf als Funktion vom Ort x zur Zeit t = t, (siehe Abb. 6b) wieder. 

An Hand der vorweggenommenen Ergebnisse soll hier noch kurz das Lö- 
sungsverfahren selbst erläutert werden. 

Man geht davon aus, daß p und w an der kritischen Stelle x, (Kontakt- 
fläche) stetig sind. Druck und Strömungsgeschwindigkeit sind hinter der Stoß- 
welle identisch mit den entsprechenden Größen hinter der Verdünnungswelle. 
Es sind also zunächst p, und u, aus den vorgegebenen Anfangszustanden 
2 und 0 zu ermitteln. Dies geschieht dadurch, daß man zunächst alle Zustände 
1’ = fi (pi, ui) sucht, die mit 0 durch eine Stoßwelle gekoppelt sein können. 


Die Gesamtheit dieser Zustände läßt sich durch u 
eine Kurve pj = f5(w}. 0) darstellen. Sodann dp 
hat man die Gesamtheit aller Zustände 1” dar- ß “3 Tt 
zustellen, die mit 2 durch Verdünnungswelle 
gekoppelt sein können: 1 2). Der f3(u,£) 
die gesuchten Werte p, und u, (Abb. 7). Aus ne 
den Stoßwellenbeziehungen folgen nun noch die f2(u,.0) 
Stoßwellengeschwindigkeit D und die Dichte im fy 

Bereich x, < x < x, auf der StoBwellenseite der f rn 
Kontaktflache. Aus p, folgt die Temperatur U; 
T,,. Aus den Gleichungen für eine Verdiin- Abb.7 
nungswelle ergibt sich andererseits die im 


Bereich 2,< x < ag, daraus nach cj = 7 Pı die Geschwindigkeit Uy — Cy, des 
Av 


Auslaufs der Verdünnungswelle sowie nach = p ‘tie Temperatur 
0. 


Die Funktionen f2(u}, 0) bzw. fs (wy, 2) folgen aus Rankine-Hugoniot- 
schen StoBwellenbeziehungen bzw. fiir die V erdünnungswelle aus der Existenz 
einer Rieman-Invarianten bei ‚einfachen Wellen‘, zu denen auch die zen- 
trierte Verdünnungswelle gehört (siehe "!)). Damit ist der Vorgang im ebenen 
Stoßwellenrohr vollständig beschrieben. 


Welche Folgerungen ergeben sich nun, wenn man den mit der Ausbreitung 
eines Funkenkanals verbundenen Ausgleichsvorgang durch ein solches Modell 
beschreibt ? 

Es zeigt sich zunächst, daß die Kanalbegrenzung in ihrem Antangsstadium 
im Gegensatz zu der bisher üblichen Annahme nicht mit der Stoßwellen- 
front, sondern mit der Kontaktfläche zusammenfällt. Beide Flächen pflanzen 
sich von vornherein mit unterschiedlicher Geschwindigkeit fort (Abb. 6b), 
und zwar bewegt sich die Kontaktfläche etwa um 20%, langsamer als die 


Stoßwellenfront. (Für starke Stoßwellen gilt u, = „ : | P). Ferner ist die Zahl 


der Teilchen im Kanal während der Expansion zunächst unveränderlich, da 
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der Massenstrom durch eine Kontaktflache verschwindet. Dies gilt jedoch 
nur so lange, bis die in der Achse des Funkenkanals reflektierte Verdiinnungs- 
welle die Kontaktfläche einholt und an dieser nur zum Teil reflektiert wird. 
Der Teil der Verdiinnungswelle, der die Kontaktfliche passiert, reiBt einen 
Teilchenstrom mit sich, wodurch sich die Zahl der Teilchen innerhalb der Kanal- 
begrenzung stark verringert. 

Zu diesem Zeitpunkt sind die im Realfall vorliegenden Verhältnisse mit 
diesem Modell jedoch nur noch sehr bedingt deutbar. Erstens ist eine Kon- 
taktfläche ihrer Natur nach nicht stabil, sondern zerfließt, wie dies in den 
Abbn. 6ce und 6d durch die gestrichelten Linien angedeutet ist. Ferner kann 
man die durch den Stromfluß bedingten Vorgänge nicht mehr außer acht 
lassen, da durch sie die Kanalausbreitung auf einen bestimmten Radius be- 
grenzt ist. Dies führt dazu, daß die Kontaktfläche aufspaltet: Ein Teil pflanzt 
sich als ein Bereich zunächst starken, mit der Zeit jedoch immer breiter werden- 
den Dichtegefälles mit der Geschwindigkeit der alten Kontaktfläche fort. 
Dieses Gebiet wird beobachtet und ist von Drabkina!) als „Kanalhülle“ 
bezeichnet worden. Der andere Teil bildet die endgültige Kanalbegrenzung 
und stellt im wesentlichen ein Gebiet starken Temperaturgefilles dar. Im 
Gegensatz zu der hydrodynamisch bestimmten Bewegung der ‚Hülle‘ wird 
die Bewegung der Kanalbegrenzung vom Zeitpunkt ihrer Herausbildung an 
von anderen Gesetzen bestimmt, in denen der Stromverlauf die wesentliche 
Rolle spielt. Darauf soll hier jedoch nicht näher eingegangen werden (siehe 
auch Braginski’) und Kühn®)). 

Es muß noch erwähnt werden, daß beim ebenen Stoßwellenrohr die Ab- 
hängigkeit sowohl der Stoßwellengeschwindigkeit als auch der Geschwindig- 
keit der Kontaktfläche (= u,) von der Ausgangsdichte o, im Bereich 


10-2? < p< 10 Torr (H,) annähernd durch 

uy [m/s] wiedergegeben werden kann. Dies 
x entspricht der im Abschnitt 2 auf anderem 

Wege gewonnenen Abhängigkeit (vgl. 2.14 

und 2.15 für m = 1). Da eine Auflösung 

der maßgeblichen Gleichung nicht mög- 

lieh ist, wurde u, (09) bzw. D(o,) bei festen 

jo? 107 1 0 3, Werten für 7, und damit p, graphisch 


Me aufgetragen (Abb. 8). Zur Wahl von T, 
gestrichelte Gerade entspricht 9-1/3 liegen Anhaltspunkte aus der Literatur 
vor’). Die für die Stoßwellengeschwindig- 
keit berechneten Werte entsprechen etwa den von uns gemessenen mittleren 
Geschwindigkeiten. 
Damit soll die Diskussion des vorgeschlagenen Modells fiir das Anfangs- 
stadium der Funkenkanalentwicklung bzw. der StoBwellenausbreitung abge- 
schlossen werden. Die Anlehnung an das klassische Stoßwellenrohr liefert 
eine klare Vorstellung vom Beginn des Ausgleichsvorgangs nach erfolgtem 
Durchschlag. Es sollte möglich sein, in dieses Modell auch die Plasmaeigen- 
schaften des Durchschlagsfadens und des aus dem Kanal herausgezogenen 
Gases einzufügen. 


Berlin, Phys.-Techn. Institut der Deutschen Akademie der Wissenschaften, 
Bereich Strahlungsquellen. 
Bei der Redaktion eingegangen am 23. November 19(0. 
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Untersuchungen an Funkenentladungen 
in T-förmigen Entladungsgefäßen 
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> mit 
Kon- 


1 den 


kann H.Heß, K. Kischel, H. Morgenroth, K. Rademacher und K. Seliger 
acht 

s be- 2 Mit 11 Abbildungen 7 

lanzt 

den- con Inhaltsübersicht 

fort. Es wird über oszillographische Messungen von Strom und Spannung bei 

ülle“ einer Funkenentladung in einem T-förmigen Entladungsgefäß, das zuerst von 

zung Fowler beschrieben wurde!), sowie über Aufnahmen der Stoßwellenausbrei- 
Im tung mit einer Drehkamera berichtet. 

wird Als Füllgas diente Argon; der Anfangsdruck im Entladungsgefäß wurde 

gon zwischen 0.2 und 500 Torr variiert. 

liche Aus den Oszillogrammen werden Leistung und Energie berechnet; die 


iehe Aufnahmen der Drehkamera ergeben die Stoßgeschwindigkeit in Abhängig- 
keit von der Zeit und vom Ort. 


Ab- Es wird gefunden, daß die bis zu einer bestimmten Zeit im Funken ver- 
dig- brauchte Energie druckabhängig ist und bei etwa 5 Torr einen minimalen Wert 
eich annimmt. 
ırch Der zeitliche Verlauf der Stoßgeschwindigkeit, die Beziehung zwischen 
Dies dem zeitlichen Auftreten der Maxima in Leistung und Stoßgeschwindigkeit 
- sowie die Abhängigkeit der maximalen Stoßgeschwindigkeit vom Druck werden 
2.14 mit Arbeiten von Kolb?) und Morgenroth u. a.%) verglichen. 
En Es zeigt sich, daß die in ?) verwendeten Modellvorstellungen eine Inter- 
a pretation unserer experimentellen Ergebnisse gestatten. 
isch 
lig. 1. Experimentelle Anordnung 
ren Die Kapazität des Kondensators, der die Entladung speiste, betrug 10 u F, 
die Ladespannung stets 6 kV (Abb. 1). 

1gs- Das Entladungsgefäß war 7'-formig; das Stoßrohr aus Hartglas hatte eine 
ge- Länge von 15 cm und einen Innendurchmesser von 1 em. Die Elektroden be-  _ 
fert standen aus Wolfram (2 0,45 em) mit halbkugeliger Oberfläche; der Elek- 
em trodenabstand betrug 1 cm. 
en- Der Anfangsdruck des als Füllgas verwendeten Argon wurde zwischen 
1en 0,2 und 500 Torr variiert. ol 

4 


1) Fowler, Goldstein u. Clotfelder, Physic. Rev. 82, 879 (1951). 
en, 2) A.C. Kolb, Physic. Rev. 107, 345 (1957). 

3) H. Morgenroth, H. Heß, K.Kischel, K. Rademacher u. K. Seliger, Ann. 
Physik .. 
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Der Strom wurde mit Hilfe eines induktionsarmen Widerstandes oszillo- 
graphisch gemessen. 

Zur Spannungsmessung wurde ein Widerstandsteiler verwendet. Die Un- 
abhängigkeit des Teilerverhältnisses von der Frequenz wurde durch eine Prüfung 
mit Rechteckwellen hinreichend gut bestätigt*)°)®). Ein koaxialer Aufbau des 
Teilers, sowie eine sehr klein gehaltene Leiterschleife beim Anschluß an die 
Funkenstrecke verringerten 
die durch Induktionswirkung 
auftretende Spannung. 

Diese induzierte Spannung 
\ konnte auf 10% der Meß- 
spannung herabgedrückt und 
durch eine zusätzliche Mes- 
sung berücksichtigt werden. 

Für die oszillographischen 


AL 


c Messungen wurde der Im- 

- & pulsoszillograph OG 1—8 vom 

=: VEB Funkwerk Köpenick 
Abb. 1. Schematische Darstellung der Versuchs- benutzt. 


, anordnung. a Entladungsgefäß, 5 Wolframelektro- 
den, e Schaltfunkenstrecke, d Strommeßwiderstand, 
e Spannungsteiler, zum Ozillographen,gzur Pumpe, 


Die Stoßwellenausbreitung 
wurde mit einer Drehkamera 


a h Koaxialkabel beobachtet. 
2. Ergebnisse der oszillographischen Messungen 
2.1. Strom 


Der Strom als Funktion der Zeit ist eine schwach gedämpfte Schwingung. 
Im allgemeinen erlischt die Entladung nach 5 Schwingungen, da dann die 
verbliebene Spannung nicht mehr zum Zünden ausreicht (Abb. 2). 

Die Frequenz der Schwingung beträgt etwa 30 kHz. Daraus folgt für die 


d 


za 


| Selbstinduktivität des Entladungskreises in guter Näherung Z = 2,75 wH. 
# | . Der Stromwert im ersten 
r . . 
u Meßgenauigkeit isteine Druck- 
0 abhangigkeit nicht festzu- 
tiysec) stellen. 
u | Nach der Theorie für den 
no. | ungedämpften Kreis gilt: 


ie ot Abb. 2. Entladungsstrom als Funktion der Zeit Jin=U, Vi = 11.4kA. 


Erwartungsgemäß ist der experimentelle Wert infolge der Dämpfung insbe- 
sondere durch den Funkenwiderstand kleiner als der für den ungedämpften 
Kreis berechnete. 

1) J.S.Abramson u. N.M. Gegetschkori, J. exp. theor. Phys. 21, 484 (1951). 

5) J. Müller, FTZ 4, 211 (1951). 

6) J.Czech, Oszillographenmeßtechnik Berlin-Borsigwalde 1959, Verlag für Foto- 
 Kino-Radio-Technik. 
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Da der Funkenwiderstand infolge seines im zeitlichen Mittel recht kleinen 
Wertes nicht sehr stark in den Ausdruck fiir das Strommaximum beim ge- 
dämpften Kreis eingeht, ist von vornherein keine bedeutende Druckabhängig- 
keit des Stromes zu erwarten, obwohl der Widerstand nach unseren Messungen 
eine solehe Abhängigkeit zeigt. RER, © 


2.2. Spannung an der Funkenstrecke 


Nach der Zündung des Schaltfunkens steigt die Spannung an der Meß- 
funkenstrecke fast bis auf die Ladespannung an. Dieser Wert wird während 
der von Druck und Stoßfaktor abhängigen Verzugszeit beibehalten, bis der 
Spannungszusammenbruch erfolgt. Die Zeit dieses Spannungsabfalls ist eben- 


190: 
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Abb. 3. Spannung an der Funkenstrecke Abb. 4. Spannung beim 1. Strommaximum 
als Funktion der Zeit als Funktion des Druckes 


falls vom Druck und vom Stoßfaktor abhängig. Die in Abb. 3 gezeigte Kurve 
setzt erst im sogenannten Bogenstadium des Funkens ein; hier ist die Span- 
nung nahezu unabhängig vom Strom. 

Mit abnehmendem Strom erlischt die Entladung schließlich, die Spannung 
steigt sprunghaft in entgegengesetzter Richtung bis zum Wert für den Wieder- 
einsatz der Entladung an, und das Bogenstadium beginnt erneut. 

(Verwendet man statt Wolfram Messing als Elektrodenmaterial, so findet 
man regelrechte „Zündspitzen‘“ in der Spannung zu Beginn jeder Halb- 
welle *) 8) 9). 

Der zeitliche Verlauf der Spannung ist vom Anfangsdruck im Entladungs- 
gefäß abhängig. Abb. 4 zeigt die Druckabhängigkeit der Spannung zur Zeit 
des ersten Strommaximums. 

Da beim Strom keine Druckabhängigkeit gefunden wurde, zeigen auch 
Widerstand und elektrische Leistung der Funkenstrecke zur Zeit des Strom- 
maximums als Druckfunktionen einen ähnlichen Verlauf wie die Spannung. 


2.3. Funkenenergie 


Die bis zum Auftreten des ersten Strommaximums im Funken verbrauchte 
elektrische Energie ergibt sich aus den Oszillogrammen für Strom und Span- 


*) H. Kaiser u. A. Wallraff, Ann. Physik 34, 297 (1939). 7 

8) A. Wallraff, Z. techn. Physik 17, 44 (1936). . “8 

®) G. Glaser, Z. Physik 148, 44 (1955). u 
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nung nach 


_ wobei 7 die Schwingungsdauer ist. 
Abb. 5 zeigt Strom, Spannung, Leistung, Energie und Widerstand der 
Funkenstrecke als Funktionen der Zeit. 


Die bis t = z in der Funkenstrecke freiwerdende Energie durchläuft in 
Abhängigkeit vom Anfangsdruck 


im EntladungsgefiB bei etwa 
5-10"! ft ei 5 Torr ein Minimum (Abb. 6). 
Dieses Verhalten läßt an Hand 

Druckabhängigkeit des Stoß- 
aktors (Verhältnis von Stoßspan- 

100 ; nung, bei uns etwa gleich der 
Ladespannung, zu statischer Durch- 


schlagsspannung) und seines Ein- 
flusses auf die Zusammenbruchszeit 
der Spannung an der Funkenstrecke 


erkliren, 


E (Wsec) 
m: 
10 
9: 
8: N 
7: 
6; 
Abb. 5. Strom, Spannung, Leistung, Ener- Abb. 6. Funkenenergie bis zum 1. Strom- 
gie und Widerstand der Funkenstrecke als maximum als Funktion des Druckes 
Zeitfunktion 
ies 2.4. Energiebilanz zur Zeit des ersten Strommaximums 
7 Die gesamte im Kondensator zu Beginn der Entladung gespeicherte 


‘ Energie ergibt sich aus Eye, = C Up zu 180 Wsec. 


= 

Zur Zeit t = 7 ist diese Energie folgendermaßen verteilt: 270. Fo 
Fr 

= 


Im Kondensator ist keine Energie gespeichert. 
1 
> Im Magnetfeld des Entladungskreises befinden sich nach E, = — LP | 
etwa 134 Wsec. 
- Im Shunt sind rund 2 Wsee verbraucht worden, während in der Funken- 
strecke .in Abhängigkeit vom Druck zwischen 6 und 12 Wsec umgesetzt 
worden sind. 
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Der Rest zur Gesamtenergie von rund 35 Wsec entfällt auf den Verbrauch 
in der Schaltfunkenstrecke, auf den Widerstand der Zuleitung und auf dielek- 
trische Verluste im Kondensator. 

Es werden also z. B. bei einem Anfangsdruck von 500 Torr bis zum ersten 
Strommaximum 46 Wsec im ganzen Kreis verbraucht — das sind 25%, der 
Gesamtenergie. Davon entfallen 12 Wsec auf die Funkenstrecke — das sind 


26% der bis t = > verbrauchten Energie oder 7%, der Gesamtenergie. 


3. Ergebnisse der Aufnahmen mit der Drehkamera 

In Abb. 7 ist die Stoßgeschwindigkeit als Funktion des Abstandes der 
Stoßfront vom Funkenkanal dargestellt. 

Die Geschwindigkeit wächst zunächst mit 
zunehmendem Abstand an, durchläuft ein Maxi- 
mum und nimmt dann monoton ab. Dieser Ver- 
lauf wird auch in ?) angegeben. Bei unseren Unter- 
suchungen trat das Maximum sowohl mit magne- 
tischer Treibung als auch ohne diese auf. 

Man kann bei der Betrachtung der Aufnahmen 
mit der Drehkamera nicht von vornherein die 
Front des leuchtenden Plasmas mit der Stoß- 
wellenfront gleichsetzen. 


5 ! Dikm/sec) - 


i ; = 7Torr 
Das ist nur dann gestattet, wenn die Leucht- 2\ ch 
front das Ende des Stoßrohres erreicht. 0 5 rkml 
Nur in diesem Falle kann man annehmen, Abb. 7. Stoßgeschwindigkeit 


daß die Leuchtfront längs des ganzen Weges mit ls Funktion des Abstandes 
der Stoßfront identisch ist. Kommt die Be- 
wegung der Leuchtfront jedoch schon vor Erreichen des Rohrendes zum 
Stillstand, soist das ein sicheres Zeichen dafiir, daB sich die StoBwelle schon 
eine gewisse Zeit vorher von der Leuchtfront gelést hat und dieser voraus- 
geeilt ist. Die reflektierte StoBwelle hindert dann die leuchtende Plasma- 
front am Erreichen des Rohrendes. 

Dieses Verhalten findet man bei langen 
Rohren oder hohen Drucken. Bei unseren Mes- 
sungen ist die Stoßfront bis über das Maximum 
der Stoßgeschwindigkeit hinaus mit der Leucht- 
front identisch (Abb. 8). Die Meßwerte für die 
Stoßgeschwindigkeit streuen in einem Bereich 
von — 10%, um den Mittelwert. 

Die Stoßgeschwindigkeit erreicht bei niedri- 
gen Drucken ihr Maximum bei einem Abstand 
von etwa 2cm vom Funkenkanal. Zwischen 2 
und 5Torr verringert sich dieser Abstand sprung- 
haft auf 1 cm und nimmt dann mit zunehmen- apb.8. ‘Aufnahme des Stoß- 
dem Druck bis auf 0,5 cm ab (Abb. 9). rohres mit der Drehkamera; 

Trägt man die Stoßgeschwindigkeit über der Po = 1Torr. Untere Begren- 
Zeit auf, die seit Beginn der Entladung ver- Bode der Leuchterscheinung: 

oden der Entladungskammer ; 
gangen ist, so erhält man die ausgezogene — obere Begrenzung: Stoßrohr- 
Kurve in Abb. 10. 


ende 


7 
| = 
4 


darzustellen, 
klären. 


ergibt sich nach 


ladungsgefäß, E die in der 
eine Konstante. 
Für die Energiefunktion macht 


(-5) 


Abb. 9 
keitsmaximums vom Funkenkanal 


Der in ?) angegebene Ausdruck 


.. 


enthalten. u 
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Über die in 2) angeführte Beziehung für die Stoßgeschwindigkeit hinnsia 


die nur in der Lage ist, den zu größeren Zeiten hin abfallenden Teil der Kurve 
3) die Möglichkeit, auch das Maximum zu er- 


gefundenen Verlauf näherungsweise wiedergibt: 


Abstand des Stoßgeschwindig- 


D(t)~ ( =A) 3 
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Berücksichtigt man, daß die Energie im Funkenkanal nicht momentan 
umgesetzt wird, sondern eine Funktion der Zeit ist, so erhält man im Falle 
ebener Symmetrie nach ?) aus hydrodynamischen Betrachtungen: 


Dabei ist qg der Querschnitt des Stoßrohre: 


t 
VE(t) | | E(t) dt 


8, 0) die Anfangsdichte im Ent- 


Entladungsstrecke verbrauchte Energie und a 


man einen Ansatz, der den experimentell 


=2r 


kit) = 4 (2) 
Hier ist T die Zeit, zu der die Leistung 4 O/km/sec) 
ihren maximalen Wert annimmt. 6! 
Der Zusammenhang zwischen Tr . 
und der Zeit, zu der das Maximum | 
der StoBgeschwindigkeit auftritt, | 
lautet: 54 
4 | 
an 3 bi (3) | / 
Die experimentellen Werte genügen / 
etwa dieser Beziehung, eine Druck- yi 
abhängigkeit ist nicht zu erkennen. i. 
3 
7 | 
| 
9 5 10 3 8 
erperımenteller Verlau tly sec) 
Morr) „ — —- Verlauf nach (7) 


Verlauf nach (4) 


Abb. 10. Stoßgeschwindigkeit als Funk- 
tion der Zeit; p = 1 Torr 
6 


(4) 


% 


ist in der Gleichung (1) als Spezialfall für die momentane Energieumsetzung 
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Paßt man die Funktion (1) unter Verwendung des in (2) angegebenen 
Energieverlaufs im Maximum an die experimentelle Kurve an, so erhält man 
die gestrichelte Kurve der Abb. 10. 

Die Beziehung (4) führt auf die strichpunktierte Kurve. 

Man sieht, daß die Gleichung (1) mit der Annahme (2) eine wesentlich 
vollständigere Annäherung an den experimentellen Verlauf ergibt als die bisher 
ausschließlich verwendete Beziehung (4). 

Aus (1) folgt für ¢ = const eine Abhängigkeit der Stoßgeschwindigkeit 
vom Druck in der Form 

D, -const 0% (5) 
Abb. 11 zeigt die maximale StoBgeschwindigkeit (die etwa stets zur gleichen 
Zeit nach Beginn der Entladung auftritt) in Abhängigkeit vom Druck. Die 


estrichelte Gerade entspricht dem 
8 Dna,(km/sec) 


1 
Exponenten — — 


Das ist auch zu erwarten, da die 
Druckabhängigkeit der Energie nicht | 
— 
allzu groß ist und A nur mit dem 07 ] 0 100 


Abb. 11. Maximale Stoßgeschwindigkeit 
als Funktion des Druckes 


Obwohl man nach unseren Mes- | 
sungen annehmen muß, daß auch A : 
druckabhangig ist, gibt (5) die ex- ; 
perimentellen Werte gut wieder. 
| 
| 


plTorr) 


Exponenten + in (1) eingeht. 


Die Messungen ergaben, daB die bis zum ersten Strommaximum im Funken 
verbrauchte Energie in Abhängigkeit vom Anfangsdruck im Entladungsgefäß 
bei etwa 5 Torr einen minimalen Wert annimmt. Dieses Verhalten läßt sich aus 
der Abhängigkeit der Zusammenbruchszeit der Spannung von Druck und 
Stoßfaktor erklären. 

Der experimentelle Verlauf der Stoßgeschwindigkeit wurde mit der in *) 
angeführten Beziehung (1) verglichen, die eine wesentliche Erweiterung des 
bisher verwendeten Ausdrucks (4) darstellt. Es zeigte sich, daß die Funktion 
(1) den experimentellen Verlauf bei entsprechender Anpassung recht gut be- 
schreibt. 


Berlin, Physikalisch-Technisches Institut der Deutschen Akademie der 
Wissenschaften, Bereich Strahlungsquellen. 


Bei der Redaktion eingegangen am 23. November 1960. 
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Zur Theorie der Feldemission aus Halbleitern 
Von R. Fischer, H. Neumann und Ch. Kleint 
Mit 2 Abbildungen 
Die Theorie der Feldemission aus Halbleitern wird auf p-Halbleiter für 
beliebige Temperaturen erweitert. Formeln für die Feldeindringung werden 
hergeleitet. Oberflächenzustände wurden nicht berücksichtigt. 
= 
1. Einleitung 
Zur Theorie der Feldemission aus Halbleitern sind Arbeiten von Morgu- 
: list), Busch?), Stratton®), Wassiljew*) und Elinson®) bekannt, die die 
N Feldemissionsstromdichte in Abhängigkeit von der Feldstärke, der Tempe- 
3 ratur und einigen charakteristischen Größen des Halbleiters angeben. In allen 
3 Fallen beschrankte sich die Rechnung auf n-Halbleiter. In dieser Arbeit werden 
; Gleichungen für die Stromdichte aus p-Halbleitern gegeben und einige Un- 
5 richtigkeiten in der Arbeit Wassiljews*) korrigiert. Oberflächenzustände 
werden nicht berücksichtigt, eine Emission aus Oberflächenzuständen, Stör- 
. termen oder dem Valenzband wird vernachlässigt. Abb 
x Die durchgeführten Rechnungen und die angenommenen Näherungen sind 
k . . . 
N an die entsprechenden von Murphy und Good®) angelehnt, die die Feld- Nine 
: emission aus Metallen untersuchen. In die Grundgleichung | Felde 
Absta 
werden die Durchdringungswahrscheinlichkeit D(E,, F) der Potentialschwelle (e Di 


vor der Halbleiteroberfläche (s. Abb. 1) und die Zahl N(E,, T)dE, der pro (x) 
Sekunde auf die Flächeneinheit der Oberfläche treffenden Elektronen mit 
einer Gesamtenergie zwischen EZ, und E, + dE, eingesetzt, und die Integration 
von E, = —y, der unteren Kante des Leitungsbandes, bis E, = oo liefert 
dann unter den angegebenen Voraussetzungen die Stromdichte j. (E, Gesamt- 
energie der Elektronen, F Feldstärke, e Absolutwert der Elementarladung, mit | 
T Temperatur, y Elektronenaffinität; s. Abb. 1.) a 


1) N. D. Morgulis, 17, 983 (1947). 


2) G. Busch u. T. Fischer, Brown Boveri Mitt. 45, 532 (1958). _ vu), 
3) R. Stratton, Proc. Phys. Soc. London (B) 68, 746 (1955). j gegel 
') G. F. Wassiljew, Paanorexunka n Iektponnka 8, 962 (1958), liert. 
>, M.I. Elinson, u Jaektponnka 4, 140 (1959). 
6) E.L. Murphy u. R. H. Good jr., Physic. Rev. 102, 1464 (1956). ) 
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D(E,,F) kann mit Hilfe der WBK-Methode berechnet werden und ist 
(Lit. in ®)) 
2i 
D(E,, F) | i+ exp | p(&) asl : (2) 
4 
p(x), die Impulskomponente der Elektronen in x-Richtung, ergibt sich aus 
Gesamtenergie E, und potentieller Energie zu : 
P(x) = 2m [Bet tere] 
Haibleiter vakuum 
den 00 
D 
_ 
gu- 2 = +A 
die 
1pe- 4x 
len - @° 
‘den 
Un- 
nde 
tör- 
Abb. 1. Bändermodell mit Feldeindringung (nichtentartetes Elektronengas). E, Ge- 
sind samtenergie der Elektronen (potentielle Energie infolge Makropotential + Kristallener- 
eld gie), untere Kante des Leitungsbandes, Ey, obere Kante des Valenzbandes, Ey Fermi- 
is Niveau, A Akzeptoren, E , Ionisierungsenergie der Akzeptoren, AR Bandverbiegung infolge 
Feldeindringung, y Elektronenaffinität, ® Austrittsarbeit, ©, ¢, Betrag des energetischen 
Abstandes des Fermi-Niveaus von der unteren Kante des Leitungsbandes bzw. der oberen 
(1) Kante des Valenzbandes im Halbleiter-Inneren 


wie | € Dielektrizitätskonstante des Halbleiters). x, und x, sind die Nullstellen von 
p(x) (s. Abb. 1), die von dem gewählten E, abhängen; der zweite bzw. dritte 


u Term rechts berücksichtigen Bildkraft und anliegendes Feld. u 6 
tion Integration im Exponenten liefert für D(E,.F) u 
ix, F) = |1+ ex :) y 
mt- | D(E,, F) (ase) (3) 
ing, mit den Abkiirzungen 
) 2 e-F 
e+ 1 
v(u), die sog. Nordheimfunktion, ist bei Murphy und Good®) genau an- 
gegeben und für 0 < u < 1 bei Burgess, Kroemer und Houston’) tabu- 
liert. 
7) B. E. Burgess, H. Kroemer u. J.M. Houston, Physic. Rev. 90, 515 (1953). 


| 
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Gleichung (3) unterscheidet sich von der entsprechenden bei Wassiljew 
(4), Gl.(1; 3)]. Wie die Rechnung ergibt, darf die Größe x nur im Argument 
Nordheimfunktion v(u) auftreten. 


2. Nichtentartetes Elektronengas 


Bei nichtentartetem Elektronengas gilt die Maxwell- Boltzmann -Stati- 


 stik und N(E,. T)dE, ergibt sich dann zu 
bs 
n | E, + 


u N(E,. T) dE, = dE,. (4) 


2 
(2amkT)'” \ kT | 
Hier und im folgenden wird fiir die effektive Masse stets die Masse des freien 
Elektrons eingesetzt. n ist die Elektronenkonzentration an der Oberfläche im 
Halbleiterinneren, m die Elektronenmasse und k die Boltzmann -Konstante. 


Eine exakte Lösung des Integrals (1) nach Einsetzen von (3) und (4) ist 
nieht möglich. Der Integrand wird deshalb an dem Punkt E, = n, wo er ein 
rer besitzt. entwickelt. Diese Entwicklung hat aber nur solange Sinn, 
wie |y |) S = y ist. Bei || > y würde zwar der Integrand sein Maximum formal 
bei i = | > y besitzen, da aber im verbotenen Band, also bei |E, > y 
laut ph PT A keine Elektronen emittiert werden können, wird das tatsäch- 
liche Maximum bei |E,| = y liegen (Emission aus dem Valenzband wird auch 
ohne die obige Annahme i. a. keine Rolle spielen, da dort D(E,,F) zu klein 
ist). Man wird also für 4) > y an der Stelle E, = — y entwickeln. 


Für n ergibt eine einfache Rechnung 


> 
- 6 1/2 (H. ) (9) 
Ta 8(k T)? ) - 
mit 
Solange = y, also 
F< 2)2y'2k T t(x F E (H. E.) (6) 
wird an der Stelle 4 entwickelt, und die an Murphy und Good®) angelehnte 


a Rechnung ergibt fiir die Stromdichte 


- n-F 1 1 
= 1/21 x Giz 
(4) 
I, 
; >. 8) In allen mit (H. E.) gekennzeichneten Gleichungen sind sämtliche Größen in 


Pant, Hartree -Kinheiten einzusetzen. Näheres dazu und Definitionen in ®). 


Al 


R. 


Dal 


muß 
erhä 


Für 


d.h. 
artuı 
eben 
Taı 


SOWO 


Feld 


dring 


Ng ist 


oll 
1h 
= 
| 
| 
Fan 
| 
| 
t(u) 
wer 
dort 
zum 
| Bei 
|| 
| 
Es se 


R. Fischer, H. Neumann u. Ch. Kleint, Zur Theorie der Feldemission aus Halbleitern 199 


ljew | Dabei ist h das Plancksche Wirkungsquantum und @(u) ist gegeben durch 


ment = 
a u) = — 2 v(u) (wu) 


2kT 4 


G (2) = J exp {— 27} dz = (y +) (H. E.) 


Stati- 


t(u) und @(u) sind von v(u) abgeleitet und können ®) oder 7) entnommen 
(4) werden. 


Die Bedingung (6) ist bei Wassiljew [4), Gl. (18)] in anderer Form gegeben, 
dort ist auch 7) > y möglich. Solche 7-Werte müssen jedoch für Halbleiter, 
freien | zumindest unter den hier gemachten Voraussetzungen, ausgeschieden werden. 


he im (Bei Stratton*) wird an der Stelle E, = — y entwickelt.) 
tante. 
Sobald 
Er ), (HE) (8) 
Sinn, 
formal ) muß also an der Stelle E, = — y entwickelt werden. Für die Stromdichte 
>y! erhält man dann 
1 auch n t( y ) 
1 klein QakT)2 {kT F (9) 
4 | exp (« (H. E.) 
(5) FS 2) 2 yee rex ; ), (H. E.) 


d.h. für große F oder kleine 7, geht Gl. (9) in die bei Stratton fiir Nichtent- 
artung angegebene Formel iiber|*), G. (29)]. Das war zu erwarten, da Stratton 
| ebenfalls an der Stelle 2, — y entwickelt, die Rechnung aber nur für kleine 


T ausführt. 


(6) | Die Elektronenkonzentration » an der Oberfläche im Halbleiterinnern, die 
sowohl in Gl. (7) als auch in (9) vorkommt, kann unter Berücksichtigung der 
Feldeindringung (Abschn. 4) berechnet werden. Ist AR die durch Feldein- 

elehnte | dringung verursachte Bandverbiegung (Abb. 1), dann ist 


(7) | % ist die Elektronenkonzentration im ungestörten Halbleiter. 
4 3. Entartetes Elektronengas 
Sen is Für diesen Fall gilt die schon bei Wassiljew*) durchgeführte Rechnung- 


Es sei deshalb nur das Ergebnis unter Berücksichtigung der gegenüber #) ver. 


| 
IR 
4 


4 ‘ 


as 


änderten Gl. (3) für D(E,, F) gegeben: 


er 1 ackT 
ap {- 4} Far 


(ack 1+ckT (ck | 1-+ck7T „| 
— exp ı— + exp ı— ¢* 
12 | er | l+ckT | “3 | 
r 
mit den Abkürzungen 
4 2 y — 2 F!: y — 2 \ 
b= =); c= 2/2 tl 
3 F y — ¢* F 
- « <* ist der Abstand der unteren Kante des Leitungsbandes vom Fermi-Niveau 
= : bei Entartung (Abb. 2). Die Kriterien für den Gültigkeitsbereich der Gl. (11) 
sind bei Wassiljew*) angegeben. 
he 
Halbleiter X Yokuum 

P 

al 

Abb.2. Bändermodell bei Entartung. [* Ab- 

nn stand des Fermi-Niveaus von der unteren Kante 
= 


des Leitungsbandes 


In den meisten Fällen genügt es, auch bei mittleren Temperaturen mit einer 
= Zl Näherung fiir 7’ — 0 zu rechnen. Aus (11) wird dann 


1 
exp {— 6} [1 — (1 + H. E. 12 
Date P )expı ( ) (12) 


(b und ¢ wie oben). Setzt man den Wert von [* ein, der für n-Halbleiter be- 
rechnet wurde, so geht Gl. (12) in die bei St ratton für Entartung angegebene 
Formel über [?, Gl. (30)]. io 


g 


gilt 


Ak 


a 
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¢* kann, wie aus Abb. 2 zu ersehen, aus AR und {, berechnet werden; es 


AR ergibt sich aus der Feldeindringung (Abschn. 4), 


Pa | 


4. Theorie der Feldeindringung 


Bei starken elektrischen Feldern dringen die Feldlinien weit in das Innere 
des Halbleiters ein und führen zu Änderungen der Ladungsträgerkonzen- 
tration an der Oberfläche. Dieser Effekt kann beschrieben werden durch eine 
Bandverbiegung AR (Abb. 1), die für n-Halbleiter in Abhängigkeit von der 
Feldstärke F bereits von Morgulis®) und Stratton?) berechnet worden ist. 
Ihre Ergebnisse gelten auch noch im Fall der Eigenleitung. Bei einem p-Halb- 
leiter müssen aber andere Beiträge in der Ladungsbilanz berücksichtigt werden. 
Bei der folgenden Rechnung wird die Konstanz der Fermi-Grenze voraus- 
gesetzt (vgl. 3)). 


Nimmt man an, daß die Donatoren-Konzentration wesentlich kleiner ist 
als die Akzeptoren-Konzentration, dann lautet die Poisson -Gleichung 


207 2 


da? € jE,—¢,+ U(zx)\ 

| k j 
(14) 

i 
U(x) ist die zusätzliche potentielle Energie, die die Teilchen durch das Makro- 
potential V erhalten (U(x) = —e:-V(x)), £, bzw. ¢, der Betrag des energe- 
tischen Abstandes des Fermi-Niveaus von der oberen Kante des Valenzbandes 
bzw. von der unteren Kante des Leitungsbandes im Innern des Halbleiters 
und E, die Ionisierungsenergie der Akzeptoren (Abb. 1). Weiter ist 


E VB LB h2 ge ar 
xi? dx 


Der erste Term in Gl. (14) gibt die Löcherkonzentration im Valenzband!P), 
der zweite die Zahl der ionisierten Akzeptoren pro cm? 1°) und der dritte die 
Elektronenkonzentration im Leitungsband!P) bei einem bestimmten Wert von 
U(x). Die Randbedingungen für Gl. (14) sind 


U(x«=0)=—AR Ua>-o)=0 


aU _ e|F(0)| dU 


N.D. Morgulis, 16, 959 (1946). 
10) O. Madelung, Handb. d. Physik (herausgegeb. von S. Fliigge) Bd. 20, 1957. 
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F (0) V/cm ist die Feldstärke unmittelbar an der Spitzenoberfläche im Vakuum. 
Gl. (14) läßt sich noch etwas vereinfachen, da bei größeren Feldstärken, die 
für Feldemission ja ohnehin nur in Frage kommen, die Löcherkonzentration 
in der Nähe der Oberfläche klein sein wird und da die Annahme gemacht 
werden kann, daß alle Akzeptoren ionisiert sind. Dann kann auf die Löcher 
im Valenzband die Maxwell-Boltzmann -Statistik angewendet werden und 
die Zahl der ionisierten Akzeptoren ist konstant gleich N 4. Für die Elektronen- 
konzentration im Leitungsband läßt sich keine Vereinfachung angeben, da 
auch der Fall der Entartung untersucht werden soll. Durch Vereinfachung 
von (14) folgt somit 


de {Po exp ler |— NiaFır( (16) 


(?, Löcherkonzentration im Valenzband des ungestörten Halbleiters). 


Die erste Integration der Gl. (16) unter Berücksichtigung der Randbedin- 


2 15) liefert: 
® gungen (15) liefer 
(eF(0))? 42e _ AR 
— exp ir) - 


Diese Gleichung ist unter den angegebenen Voraussetzungen allgemein gültig 


und kann für die Grenzfälle der Entartung und Nichtentartung des Elektronen- 
gases noch weiter vereinfacht werden durch geeignete Approximation von Fj). 


Nichtentartung bedeutet 0, das Fermi-Integral Fj). kann 


dann nach !°) gut argenähert werden durch 


Der Fehler ist für y< 1,5 kleiner als 4%. Damit wird aus Gl. (17) | y 


(e F(0))? = 8aee pok T(1-exp[-7 + Na AR (18) 
4k {log (1 + exp {- - log (1 +4 exp {= 


~ 

» 


Entartung hingegen bedeutet —"— > 0. Fie (y) wird i in diese 


10) gut angenähert durch 
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Der Fehler liegt für y > 1,5 unterhalb 1,5%. Einsetzen in Gl. (17) und Aus- 
führen der Integration liefert a 


= Po k T (1—exp {- +N, 


2 5/2 a? 2 

Da die fiir Fy). (y) verwendeten Näherungen recht gut sind, wird auch das ; Ei: 
Übergangsgebiet von Nichtentartung zu Entartung durch die Gl. (18) und a 
(19) erfaBt. 

Bei Kenntnis der notwendigen Daten kann man mit Hilfe der Gl. (18) und 
(19) zusammengehörige Werte von AR und F berechnen. AusAR ist dann n 4 
mittels Gl. (10) und ¢* mittels Gl. (13) leicht zu ermitteln. Daraus kann schließ- ‘So 
lich die Stromdichte der kalten Emission eines gegebenen p-Halbleiters in E 
Abhängigkeit von der Feldstärke berechnet werden. L 


Leipzig, Physikalisches Institut der Karl-Marx-Universität. cree a 


= 


Bei der Redaktion eingegangen am 25. 1960. 
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Feldemission aus Silizium stro1 
Von Ch. Kleint, H. Neumann und R. Fischer Ge ( 
> 
Mit 18 Abbildungen Fort 
Lost 
Es wurde die Feldemission von p-leitenden Silizium-Einkristallen durch 
Aufnahme von Stromspannungs-Kennlinien bei verschiedenen Temperaturen 
m untersucht. Die Ergebnisse werden mit der in der voranstehenden Arbeit 
3 weitergefiihrten Theorie verglichen, indem eine Anpassung der experimentellen 
R Kennlinien an die theoretischen Kurven durchgeführt wurde. Sie ergab für Wa 
2 die Temperaturabhängigkeit der Kennlinien (gerader Teil) Übereinstimmung zücl 
7 mit der Theorie bei einer aus dem Temperaturverhalten der Feldemission ım \ 
’ erschlossenen Austrittsarbeit von ® = 5,2 eV. Da Mikrogeometrie der Spitze, lang 
: Oberflächenzustände und -bedeckung nicht bekannt waren, mußten bei dem 
3 Vergleich zwischen Theorie und Experiment einige Annahmen eee werden. Kri 
= d 
3 Wahrend die kalte Emission aus Metallen recht eingehend idl 
d wurde, sind verhältnismäßig wenige Untersuchungen an Halbleitern und Wie 
a Isolatoren veröffentlicht worden. Bekannt sind Arbeiten über CdS- und CdSe pir 
ro von Apker und Taft?) und Klimin?), über LaB, von Elinson und Wassil- stel 
- jew®)4), über ZnS von Bertoldi und Kleint?°), über SiC von Busch und 
Ik Fischer®), über Ge von Elinson und Wassiljew’) und Allen®), über Te 
= von Kleint und Fischer?) und über Si von D’ Asaro”), Allen!!)12), Kleint ng \ 
und Fischer?) und Fischer!?). Die Kennlinien, die bei einem Teil der ge- 
nannten Arbeiten erhalten wurden, konnten qualitativ mit der Theorie von 
Stratton") und Wassiljew!5) in Übereinstimmung gebracht werden. Ein 
1) L.Apker u. E. Taft, Physic. Rev. 88, 1037 (1952). 
2) A.I. Klimin, *T® 27, 719 (1957). (Po 
3) M.I. Elinson u. G.F. Wassiljew, Paınotexnnka m Jaextpounka 2, 348 (1957). lad 
4) M.I. Elinson u. G.F. Wassiljew, Paınorexnnka mn Iıektponnka 3, 945 (1958). 
is 5) W. Bertoldi u. Ch. Kleint, Ann. Physik 7, 388 (1959). 
6) G. Busch u. T. Fischer, Brown Boveri Mitt. 45, 532 (1958). Sté 
7) M.I. Elinson u. G.F. Wassiljew, Paınorexunka n Iıertponnka 4, 728 (1959). sch 
a 8) F. G. Allen, Havard University Technical Report No. 237. dat 
®) Ch. Kleint u. R. Fischer, Z. Naturf. 14a, 753 (1959). 7 
4 ») L.A. D’Asaro, J. Appl. Phys. 29, 33 (1958). Sti 
1) F.G. Allen, J. Eisinger, H.D. Hagstrum u. J.T. Law, J. Appl. Phys. 30, die 
1563 (1959). 
#2) F.G. Allen, J. Phys. Chem. Solids 8, 119 (1959). 
: 13) T. Fischer, Freundl. priv. Mitt. sowie Helv. phys. Acta, 33, 961 (1960). 
0...) R. Stratton, Proc. physic. Soc. London (B) 68, 746 (1955). 
15) M.I. Elinson u. G. F. Wassiljew, Pagnorexunka 8, 962 (1958). 
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quantitativer Vergleich zwischen Theorie und Experiment ist bisher noch nicht 
durchgeführt worden. Das ist bedingt durch mangelnde Kenntnisse der ge- 
nauen Gestalt der emittierenden Spitze und der Daten über die Oberflachen- 
eigenschaften der Halbleiter. Die Resultate der bisher veröffentlichten Arbei- 
ten #)?)8)12) stimmen bezüglich der Temperaturabhängigkeit des Feldemissions- 
stromes nicht überein. Während Elinson und Wassiljew®)’) bei LaB, und 
Ge eine Temperaturabhängigkeit feststellen, wurde von Allen )12) weder bei 
Ge (bis 620 °K) noch bei Si (bis 1200 °K) eine Abhängigkeit gefunden. Durch 
Fortsetzung früherer Messungen?) wurde deshalb versucht, einen Beitrag zur 
Lösung der beiden Probleme, nämlich quantitativer Vergleich mit der Theorie 
und genaue Messung und Deutung der Temperaturabhängigkeit ¢ des ———g 
stromes, am speziellen Beispiel des Siliziums zu liefern. 2 Beer 

2. Kristalle | 


Die für die Versuche verwendeten Kristalle wurden nach einer bei H. v. 
Wartenberg!) und Johnson und Amick!”) angegebenen Methode ge- 
züchtet. Dabei wird SiCl, durch Zn reduziert; das entstehende Si scheidet sich 
im wesentlichen in Form von Nadelkristallen ab. Diese waren bis zu 15 mm 
lang und 10 bis 40 um dick; Wachstumsrichtung ist die 111-Richtung!”). 


Wesentlich für die Auswertung sind die elektrischen Eigenschaften der 
Kristalle. Der Leitungstyp wurde durch Thermokraftmessungen bestimmt. 
Alle untersuchten Kristalle waren p-leitend. Unterstützt wird dieses Ergebnis 
durch die im Funkenspektrum einiger Kristalle gefundenen Zn-Linien; Zn 
wird nach Carlson}§) in Si als Akzeptor eingebaut. 


Die sehr großen Werte von o0(7',) = 10% bis 10° 2 cm für den spezifischen 
Widerstand o bei Zimmertemperatur 7’, und der Verlauf von o = o(T) deuten 
auf große Reinheit hin. Bei Zimmertemperatur zeigten alle Kristalle Stör- 
stellenleitung und wurden bei höheren Temperaturen eigenleitend. 


Die in den Feldemissionsformeln vorkommende Elektronenkonzentration 


n, wurde ermittelt aus o und den bekannten Beziehungen: | aoe 
1 


(p, Löcherkonzentration, n; Inversionsdichte, e Absolutbetrag der Elementar- 
ladung, u, bzw. u, Beweglichkeiten für Löcher bzw. Elektronen). 


Werte für die Beweglichkeiten wurden der Literatur entnommen). Die 
Störstellenkonzentrationen konnten aus dem Temperaturverlauf von o abge- 
schätzt werden und lagen zwischen 1012 und 1014 em-*. Wenn man annimmt, 
daß im wesentlichen nur eine Sorte von Störstellen vorhanden ist, kann aus 
Störstellen- und Löcherkonzentration mit Hilfe des Massenwirkungsgesetzes 
die Ionisierungsenergie der Störstellen bestimmt werden. Sie ergab sich zu 


1) H.v. Wartenberg, Z. anorg. allg. Chem. 265, 186 (1951). : 
”) E.R. Johnson u. J.A. Amick, J. Appl. Phys. 25, 1204 (1954), 
13) R.O. Carlson, Physic. Rev. 108, 1390 (1957). { 

1) F. J. Morin u. J. P. Maita, Physic. Rev. 96, 28 (19 


urch | 
uren 
beit 
allen 
für 
Jung 
sion | 
itze, 
dem 
ıcht 
und 
dSe 
sil- 
und 4 
Te 
int | 
ge- 
von a 
Ein 
57). 
8). 
9). | 
30, 7 


pee Annalen der Physik. 7. Folge. Band 8. 1961 $a ek 


etwa 0,30 eV, was gut mit den von Carlson®) für Zn in Si ermittelten Werten 
übereinstimmt. 

Eine Abhängigkeit des Widerstandes von der Beleuchtung konnte in keinem 
Falle festgestellt werden. 


3. Experimentelle Technik 


Um die für Feldemission notwendigen hohen Feldstärken zu erzielen, 
wurden die Kristalle durch elektrolytisches Ätzen in Flußsäure mit einer feinen 
Spitze versehan. Die Spitzenradien lagen zwischen 1000 und 9000 A. Sie wur- 
den in einem hochauflösenden Lichtmikroskop gemessen; die Oberflächen- 

on konnte in einem Elektronenmikroskop untersucht werden. Die Be- 
festigung und Kontaktierung der Kristalle 
in der Versuchsröhre erfolgte meist wie bei 
Bertoldi und Kleint?). 

Eine schematische Darstellung der 
verwendeten Versuchsröhren zeigt Abb. 1. 
Der leitende Belag auf der Innenwand der 
Röhre diente zur Vermeidung von Auf- 
ladungen der Glaswand; er bestand aus 
Zinnoxyd”°) und war für sichtbares Licht 
gut durchlässig?!), so daß auch Messungen 
an beleuchteten Kristallen gemacht wer- 
den konnten. Der eingeschmolzene Metall- 
ring verhinderte bei geeigneter Polung 
eine Verfälschung der Strommessung durch 
Isolationsströme über das Glas. 

Die Temperaturabhängigkeit des Feld- 
emissionstromes wurde im Bereich von 
— 170 °C bis 400 °C gemessen. Um die 
Spitze auf die gewünschte Temperatur zu 
atl iol bringen, wurde die gesamte Röhre gekühlt 

ae oder erwärmt. Die Temperatur des Kri- 
+ stalls wurde gemessen mit Hilfe der 
Widerstandsänderung einer Wolfram- 
spirale, die sich im Wärmekontakt mit 


+. tou 


wid 


—(8) 


_ Abb.1. Schematische Darstellung einer 
_ Versuchsröhre. 1 leitender Belag, 


2 Leuchtschirm, 3 Kristall, 4 Pump- 
stutzen, 5 Anodeneinschmelzung, 
6 Glas, 7 eingeschmolzener Metallring, 
8 Zuführungsdraht zur Feldkathode 


der Spitze befand. Wie eine Abschätzung 
ergab, braucht bei den auftretenden 
Strömen eine Erhitzung der Spitze durch 
Joulesche Wärme noch nicht berück- 


sichtigt zu werden. 

Der Hochvakuumpumpstand war ähnlich dem in 5) beschriebenen. Es wurde 
bis zu 100 Stunden evakuiert bei gleichzeitigem Ausheizen der Versuchsröhre 
bei 400 bis 450 °C. Der dann noch vorhandene Restgasdruck wurde aus Messungen 
mit einem Bayard-Alpert-Manometer oder aus der Bedeckungszeit einer reinen 

 Wolframspitze22) ermittelt; es ergaben sich in allen Fällen Werte kleiner als 
10 Torr. 


2%) R. ieteen, Rev. sci. Instrum. 24, 993 (1953). 
wk. Ishiguro, T. Sasaki, T. Araiu. J. Imai, J. Phys. Soc. Japan 18, 296 (1958). 
22) M. Drechsler u. E. Heß, Physik. Verh. 4, 211 (1953). 
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Es wurde versucht, die Spitzenoberfläche durch Wärmeleitung (ähnlich 
wie bei Wolfram oder anderen Metallen), Wärmestrahlung, Elektronenbeschuß 
oder Felddesorption zu reinigen. Die Kristalle konnten bis zu schwacher Rot- 
glut erhitzt werden, jedoch konnte kein Bild einer saubsren Spitze erhalten 
werden wie etwa bei Allen!2) und Fischer). Die aus der Wachstumsrichtung 
der Kristalle (Abschn. 2) zu erwartende Bildsymmetrie war jedoch gut zu er- 
kennen. Da die Messungen sich über mehrere Stunden hinzogen, werden die 
Ströme stets von bedeckten Spitzen stammen. 

Bei der Auswertung von Messungen an Si muß außerdem berücksichtigt 
werden, daß sich beim Ausheizen der Versuchsröhre durch Anlagerung des 
vom Glas abgegebenen Bors auf der Si-Oberfläche eine p-leitende Schicht 
ausbildet®). Da aber die Kristalle p-leitend sind (Abschn. 2) und die Stör- 
stellenkonzentrationen im Innern und an der Oberfläche etwa gleich sind®), 
wird dieser Effekt vernachlässigt. 


4. Meßergebnisse 


Die hier angegebenen Meßwerte waren während der gesamten Versuchs- 
dauer reproduzierbar; auch nach Temperaturänderungen wurde bei einer 
bestimmten Temperatur stets dieselbe Kennlinie gemessen. 
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reziproke Spannung U"(kV") reziproke Spannung U'(kV™') 
Abb.2. Kennlinien mit Anlaufstrom. I: Abb.3. Gerade Kennlinien. I: Kristall 6; 
Kristall 37/3: Spitzenradius ry = 7-10-cm; II: Kristall 38:19 = 
II: Kristall 36/2: ry = 4-10 cm; III: Kristall 7 
UI: Kristall 36/1: = 5:10 cm 


2) F.G. Allen, T. M. Buck u. J.T. Law, J. Appl. Phys. 81, 979 (1960). es 
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4.1 Kennlinien bei Zimmertemperatur 


Die Messungen ergaben drei verschiedene Formen von Strom-Spannungs- 
Kennlinien. Abb. 2, 3 und 4 zeigen einige Beispiele. In Abb. 2 sind gebogene 
Kennlinien dargestellt, die bei großem 1/U einen relativ wenig von der Span- 
nung U abhängigen Strom zeigen. Dieses Verhalten wurde auch schon von 
Kleint und Fischer?) beobachtet. Die in Abb. 3 dargestellten Kennlinien 
sind über den gesamten gemessenen Strombereich Geraden. Diese Form wurde 

schon von Allen’) und Fischer?) bei Si gefunden. Abb. 4 zeigt einige 
Kennlinien, die in ihrem oberen Teil einen scharfen Knick zeigen. Über eine 
solche Form ist bisher bei Si noch nicht berichtet worden. 


107 107 o 100°K 
3 2 @ 300°K 
10°8 
‘5 
3 10? 10° 
§ \ 
un-10 
10 10°” 
= 10 \ 
I z\ 
-1, 
0 06 08 10 04 06 
reziproke Spannung U"(kV") reziproke Spannung 
Abb. 4. Kennlinien mit And der Stei- - “1 
im Teil. Kristall 40: U" (kV) 
10- em; II: Kristall37/1:r, =1-10-cem; Abb.5. Kennlinien für 100 °K und 300 °K 
III: Kristall 21: rg =3 - 105 em (Kristall 6) 


oY 4.2 Temperaturabhängigkeit des Feldemissionsstromes 

Die Abb. 5, 6 und 7 zeigen Beispiele fiir die Temperaturabhangigkeit des 
Emissionsstromes. Wie ersichtlich, fallen die Kennlinien für Zimmertemperatur, 
200 °C und 400 °C zusammen. Das ist in Übereinstimmung mit Allen’), 
der von Zimmertemperatur bis 1200 °K keine Temperaturabhängigkeit erhalten 
hatte. Bei — 170 °C hingegen ergeben sich bei gleichen Spannungen kleinere 
Ströme, und die Steigung der Geraden ist größer. 
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4.3 Beleuchtungsabhängigkeit des Feldemissionsstromes 


nungs- Bei allen durchgeführten Versuchen konnte keine Abhängigkeit des Stromes 
bogene f von der Beleuchtung beobachtet werden. Das war zu erwarten, da sich auch der 
Span. § Widerstand der Kristalle bei Beleuchtung nicht meßbar änderte (Abschn. 2). 
m von § Dieses Ergebnis steht ebenfalls in Übereinstimmung mit den Resultaten 
nlinien Allens?2). 
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(Kristall24, r= 1-10-5cm;bei 77 =300°K: und 473°K (Kristall 25: r, = 3:10>® cm; 
o=6-10°Qem, N, = 3-10" cm) bei 7 = 300 °K: 9 = 10° 2-cm, N, = 


Abb. 6. Kennlinien für 300 °K und 673 °K Abb. 7. Kennlinien für 100 °K, 300 °K 
06 
3 + 1013 cm?) 


5. Vergleich der Messungen mit der Theorie 
0 °K 5.1 Stromdichte-Formeln 


Da die verwendeten Kristalle p-leitend sind, müssen die in *) für die 
Feldemissions-Stromdichte hergeleiteten Gleichungen benutzt werden. Setzt 
man die bekannten Werte für die Konstanten h, k, m und e ein und nimmt die 

des Dielektrizitätskonstante des Siliziums zu e = 12 an), so ergeben sich mit der 
tur, Elektronenkonzentration n in em-®, 7 in °K, Fin Vem" und y, n, ¢* in eV 
112), die folgenden Stromdichte-Formeln (4), (6) und (7) Die Bedeutung der einzel- 
Iten nen Größen ist %4) zu entnehmen. 

nere 


#4) R. Fischer, H. Neumann u. Ch. Kleint, Ann. Physik (7) 8, 196 (1961). 
35) E.M.Conwell, Proc. I. R. E. 40, 1327 (1952). 
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: Aus Gl. (7) in %) wird für 
F < 0,88 - 104 7 


(3) 


(4 
+ 4,946 - 109 — 
| 


7 

ire = — = 3,75 - 10° 


y  4,67-10- 1,30 - 107 - 


3 Ge ) und ¢(u) hängen mit der Nordheim-Funktion v(u) zusammen und können 


Aus Gl. (9) in *4) wird für 
F > 0,88 - 104 T ¢(6) 
(6) 


2 Aus Gl. (12) in %) für den Fall der en des Elektronengases wird 


5.2 Berechnung von n und “hae : 
Sowohl die Micktouneniensentention n als auch ¢* [*4), Abb. 2] werden 


_ durch das Eindringen des Feldes in den Halbleiter und die damit verbundene 
_ Bandverbiegung AR [*), Abb. 1] bestimmt. Zur Berechnung von AR werden 


die in %) hergeleiteten Formeln (18) und (19) verwendet. n und £* sind gegeben 
durch 

kT 


Betrag der Energiedifferenz ,,.Ferminiveau — untere Kante des Leitungs- 


__ bandes“ im Halbleiterinneren, n, Elektronenkonzentration im Inneren). 


Um die folgenden Rechnungen übersichtlich zu halten, wird der Vergleich 

Theorie — Experiment an einem speziellen Beispiel durchgeführt (Kristall 25). 

% Ähnliche Ergebnisse erhält man bei einer analogen ee für die anderen 
Kristalle. 
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Tabelle 1 enthält die benutzten Daten für Kristall 25. 
Einsetzen dieser Werte in Gl. (18) aus %) ergibt AR in Abhängigkeit von 
der Feldstärke (Abb. 8). £* folgt aus Gl. (19) von 2%); in Abb. 9 sind die be- 
rechneten Werte aufgetragen. Mit Hilfe von AR kann jetzt auch die Elektronen- 
konzentration n in Oberflächennähe berechnet werden (Abb. 10). = 5 


Tabelle 1 (Kr. 25) 


Temperatur T | *K 100 300 473 | 673 | 1200 
Löcher- 
konzentra- |em 108 1,7- 10% | 1,2-1014 2,1 - 1016 5 1018 
tion 
Elektronen- 
konzentra- M | em [2,3 - 10% 16-107 | 1,2-10% 2,1 - 1018 5 - 1018 
tion | 
Lage der | 
Fermikante a eV 0,92 0,70 | 0,605 0,605 | 0,605 
Stérstellenkonzentration N, em? 3 - 1013 
Breite der verbotenen Zone av AE eV 1,21, nach 1) 
Elektronenaffinität wer % y eV angen. zu 3,5; 4,0; 
4,5 und 5,0 
Dielektrizitätskonstante _ 12, nach ®) 


1,0 
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Abb. 8. Bandverbiegung AR in Abhängigkeit von der Feldstärke F für Kristall 25 für 
verschiedene Temperaturen 


5.3 Die Elektronenaftinität 


Zahlenwerte für die Elektronenaffinität y für Si werden in der Literatur 
nicht angegeben. Es sind aber Messungen der Austrittsarbeit ® veröffentlicht 
worden, und da die Beziehung y = ® — €,, besteht (£, Abstand Fermikante- 
Leitungsband, s. 5.2), kann man mit den angegebenen Daten (Tab. 2) abschät- 
zen, welche Werte für y in Frage kommen. Mit Ausnahme von !l) und 12) wurden 
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n-leitende Kristalle untersucht. y dürfte damit zwischen 3,5 und 5 eV liegen, 
da die Breite des verbotenen Bandes 1,21 eV betragt. 


Fi 0,4 Tabelle 2 
2 E 03 (eV) Literatur 
~ 02 45 
gs 4,02 

38 4,6---4,9 
By 
3 T 5 
q 


Feldstarke F (10’Vcrm ) 


Abb. 9. Abstand ¢* des Ferminiveaus von der unteren >. 

Kante des Leitungsbandes bei Entartung in Abhangig- 
phe keit von der Feldstärke F für Kristall 25 ee 


5.4 Berechnung der Kennlinien 


Mit den in Tab. 1 gegebenen Werten kann j in Abhängigkeit von F mit y 
und T als Parametern berechnet werden, wobei y innerhalb der in 5.3 angege- 
benen Grenzen variiert wird. Abb. 11 zeigt den unteren Teil einer für Kr. 25 
mit 7 = 1200 °K und y = 3,5eV berechneten Kennlinie. Die ausgezogene 
Kurve wurde nach Gl. (4), die gestrichelte nach Gl. (6) ermittelt. Wie ersicht- 
lich, liegt der untere gebogene Teil bei Stromdichten, die experimentell nicht 
mehr nachweisbar sind. Für Si kann also das in Abb. 2 dargestellte Abbiegen 
der Kennlinien im unteren Teil offenbar nicht nach der benutzten Theorie 
erklärt werden, da bei kleineren 7’ und größeren y die Stromdichten noch kleiner 
werden. (Anders dagegen bei Ge und LaB,. Bei Ge ist die Breite des verbotenen 
Bandes kleiner, so daß die Elektronenkonzentration bei höheren Temperaturen 
größer als bei Si ist, und bei LaB, ist die Austrittsarbeit und damit die Elek- 
tronenaffinität bedeutend kleiner, so daß die von Elinson und Wassilje w*)’) 
gemessenen Kennlinien wahrscheinlich deutbar sind.) 


Der Übergang zur Entartung des Elektronengases äußert sich in einer 
Änderung der Steigung der Kennlinie. Jedoch liegt der theoretisch zu erwar- 
tende Knickpunkt bei Feldstärken oberhalb 108 Vem-! (Abb. 13); dieser Wert 
wurde bei unseren Untersuchungen wahrscheinlich nicht erreicht. Der Knick 
in den experimentell gemessenen Kurven (Abb. 4) muß deshalb vermutlich 
anders erklärt werden. 

Für die Anpassung theoretischer an experimentelle Kennlinien wurde 
deshalb der mittlere geradlinige Teil für nichtentartetes Elektronengas benutzt. 
Die Abb. 12, 13 und 14 stellen Berechnungen in diesem Bereich für verschiedene 
Werte von y und 7 dar. 


26) W. E. Meyerhof, Physic. Rev. 71, 727 (1947). 
27) A.H. Smith, Physic. Rev. 75, 953 (1949). 

®) L. Esaki, J. Phys. Soc. Japan 8, 347 (1953). Re oy Atay 


2) J. A. Dillon u. H.E. Farnsworth, J. Appl. nee 29, 1195 (1958). 
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Wie ersichtlich, unterscheiden sich die Kurven fiir 300 °K bis 1200 °K nur 
geringfügig, so daß bei den Messungen zwischen 300 °K und 673 °K keine oder 
nur geringe Temperaturabhängigkeit zu erwarten ist. Dagegen ist der Unter- 
schied der Kennlinie für 100 °K gegenüber den anderen deutlich zu erkennen. 
Das steht qualitativ in Übereinstimmung mit dem experimentellen Ergebnis 


(Abb. 5 bis 7). 
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Abb. 10. Elektronenkonzentration Abb. 11. Berechnete Kennlinie für Kri- 
nim Leitungsband an der Oberfläche stall25 bei nichtentartetem Elektro- 
in Abhängigkeit von der Feldstärke nengas für y = 3,5eV und T = 1200 K® 


F fiir verschiedene Temperaturen 


5.5 Durchführung der Anpassung experimenteller an theoretisch Kurven 


Um einen Vergleich zwischen Theorie und Experiment durchführen zu 
können, muß die Feldstärke F in Spitzennähe in Abhängigkeit von der Span- 
nung U und die emittierende Fläche bekannt sein. F wird durch U und die 
geometrische Gestalt der Spitze bestimmt. Nach Müller und Bahadur®®) 


%) E.W. Müller u. K. Bahadur, Physic. Rev. 102, 624 (1956). er 
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BT gilt etwa (unter der Annahme, daß die Gestalt der Spitze zwischen Kugel und 


. 
_ (Fy (Vem-!) Makrofeldstärke vor der Spitze im Abstand r (cm) vom Krüm- 
 mungsmittelpunkt der Spitze; r, (em) Krümmungsradius). 


10° 107 
Ze I 100°K I 100% 
> 105 IV 673°K ant I 673° 
3 V 1200°K & IV 1200% 
= 10° x 
— 
w 10? 
I 
8 \ S 7 
10? 
10 
| & \ 
2 10° N \ Ww 
uj IN § 


Wor 02 03 04 as 2 03 
reziproke Feldstarke F’  reziproke Feldstärke 


7, m 1 

Abb. 12. y = 4,0eV Abb. 13. y=4,5eV. Gestrichelt dieVerlänge- 
we. ye rung der Kennlinie fiir Entartung 


Allen und Mitarbeiter haben festgestellt!!), daß sich auf frisch geätztem 
Silizium Mikrospitzen befinden, die Durchmesser von 10 bis 100 A haben. 
Da unsere Kristalle nicht stark geglüht werden konnten, werden diese Mikro- 
spitzen sicher noch nicht vollständig beseitigt worden sein, die Spitzenober- 
fläche wird noch Rauhigkeiten und Unebenheiten aufweisen. Diese Mikro- 
geometrie bestimmt die tatsächlich in Spitzennähe vorhandene Feldstärke 
und wird durch Einführung eines sog. Feldverstärkungs- oder Grobfeinfaktors 
ß berücksichtigt. Die wirksame Feldstärke F ergibt sich aus der Makrofeld- 
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el und Nach Lewis!) ist ß keine Konstante, wie früher vielfach angenommen 
worden war, sondern eine Funktion des Abstandes x von der Oberfläche. Seine 
Rechnungen zeigen, daß 6 an der Oberfläche sehr groß sein kann, dann aber 
(8 mit wachsendem Abstand sehr rasch auf den Wert eins abfallt. In der Feld- 
emissionsformel ist ß(x,) einzusetzen, wo x) den Abstand von der Oberfläche 


Kr angibt, in welchem der Potentialwall durch 
das äußere Feld bis zur Höhe der Fermi- 107 
kante (bei Entartung) bzw. der unteren I 100% 
Kante des Leitungsbandes (bei Nichtent- 
artung) heruntergebogen ist. Im Gegen- 10° 300°K 
satz dazu ist fiir die Berechnung von "  473°K 
°K AR B(0) zu nehmen, da ja die unmittelbar 105 I 673°K 
°K an der Oberfläche vorhandene Feldstärke 
wirksam wird. Der Abstand x, ist bei IV 1200°K 
°K || Feldemission etwa 10 A. Da B(0) und f(a») 10° 
J%K |} in unserem Fall nicht sehr differieren, wird I 


— If der Einfachheit halber mit einem kon- 
stanten mittleren Feldverstärkungsfaktor 
gerechnet. 


Über die emittierende Fläche $ kann 
keine genaue Aussage gemacht werden. 
Sie ist abhängig von der Feldstarke*). 


Emissionsstromdichte j (Acm’?) 


Definitionsgemäß ist (I Strom; j Strom- 10’ 
dichte): \ 
J=j-8, 
bei der Anpassung der Kurven kann nur \ 
_ If darauf geachtet werden, daß S nicht zu ; 
groB oder zu klein wird. . 10° 
Die Anpassung der experimentellen an 04 
Zu die theoretische Kurve verläuft demnach reziproke Feldstärke u 
a) Unter Verwendung des gemessenen Abb. 14. y = 5,0 eV 
Spitzenradius für Kr. 25 (ry = 3-10-5cm) Abb.12—14. Berechnete Kennlinien 


. wird nach Gl. (8) für r= r, das Makro- für Kristall 25 bei Nichtentartung 
länge feld F, berechnet und die Strom-Feldstärke- [nach Gl. (6)] für ee Tem- 
Keanlinie gezeichnet. peraturen 


b) Durch Variation von 8 wird erreicht, daß berechnete und experimentell 
gemessene Kennlinie für 7’ = 300 °K dieselbe Steigung haben. 


fikro. c) Die emittierende Fläche ergibt sich nach (10). Bei ry = 3- 10-5 cm 
ober. | Wird sie etwa zwischen 10-12 und 10-9 cm? liegen. 


fikro- d) Unter Verwendung derselben Werte für # und S werden die für 100 = 5 


stärke | gemessenen Kennlinien in das j-\/F-Diagramm eingezeichnet. 


ei Die Abb. 15—18 zeigen das Ergebnis der durchgeführten Anpassung fir 
verschiedene y. Die beste Ubereinstimmung ist fiir y = 4,5 eV erreicht worden | 


()) T. J. Lewis, J. Appl. Phys. 26, 1405 (1955). 
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Die Kennlinien fiir 300 °K und 473 °K fallen nach der Theorie zusammen, so 
daß das experimentelle Ergebnis bestätigt wird. Der Elektronenaffinität y 
_ von 4,5eV entspricht bei Zimmertemperatur eine Austrittsarbeit von ® = 
5,2 eV. Das paßt zu den Messungen von Dillon und Farnsworth”), wenn 
4 man bedenkt, daß keine sauberen Oberflächen vorliegen und adsorbierte Ver. 


_ unreinigungen, vor allem Sauerstoff, eine Erhöhung der Austrittsarbeit be. 
wirken 32). 


\ a 100°K 10° a 100% 
© 300°K | \ © 300% 


@ 473°K 
10? 


© 

> 

<= 

2 

5 10 \ 
\ 
s 

& 


Emissionsstromdichte j (Acm-2) 
S, 
S 


10’ N 10° 
v\v 2 \\ 
Ä re 
reziproke Feldstärke F"' reziproke Feldstärke 
(10’cmV") 
Abb. 15. p=3,5eV; B=2,55; $=8-10-!!cm?; Abb.16. p=4,0eV; B=2,75; S=6-10-" cm’; 
theoretisch: I: 100 °K, II: 300 °K, III: theoretisch: I: 100 °K, II: 300 °K, II: 
. 473 °K, IV: 673 °K, V: 1200 °K 473 °K, IV: 673 °K, V: 1200 °K 
u 6 Diskussion des durchgeführten Vergleichs zwischen Theorie und Experiment 
6.1 Oberflächenzustände 
nt, Die in dem Vergleich benutzten Formeln wurden ohne Beriicksichtigung 
von Oberflächenzuständen abgeleitet. Letztere dürfen bei genauerer Betrach- 
_ tung jedoch nicht vernachlässigt werden, da sie bei Si über das ganze verbotene 
Band mit relativ großer Dichte verteilt sind®)%). Eine Vorhersage über Lage 


und Dichte der Zustände kann in keinem Falle gemacht werden #). Das zeigt 
sehon die Schwierigkeit, sie in die Auswertung einzubeziehen. Denn, sollen 


82) J. Bardeen, Halbleiter und Phosphore, Garmisch-Partenkirchen (1958), S. 81. 
3) H. Statz, G. A.de Mars, L. Davis u. A.Adams, Semiconductor Surface Phy- 
sics, Philadelphia (1956), S. 139. eee 
34) W. G. Litowtschenko u. O. W. Snitko, ®TT 2, 591 (1960) 


Obeı 
gen | 


(Fa 
in st 
stan 


Emissionsstromdichte j (Acm?) 


Abb 
theo 


| 
| 
\ 
| Abb 
nich 
| recl 
Obe 
tati 
eine 


Ch. Kleint, H. Neumann u. R. Fischer: Feldemission aus Silizium 217 


Oberflächenzustände berücksichtigt werden, so ist in 4) für die Randbedingun- 
gen der Poisson -Gleichung zu setzen ®) 

(F äußeres Feld, Ao Ladung/Fläche in den Oberflächenzuständen). Ao ist aber 

in starkem Maße abhängig vom Feld F sowie von Lage und Dichte der Zu- 

stände. Schon bei Berücksichtigung nur einer Art von Niveaus ist die Theorie 


„FA | 
a 100°K a 100°K 
300% © 300% 
10 \ 10° 
\ 
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10 \ 10 IT 
3 \ 
10° 10' 
S 
I a 
10° 

01 02 Q3 0% 0 01 02 0 
reziproke Feldstärke F”  reziproke Feldstärke 
07%") 

Abb. 17. y=4,5eV; B=3,2; S=5-10-"cm?; Abb. 18. = 5,0 eV; B= 4,0; S= 
theoretisch: I: 100 °K, II: 300 °K und 10-11 cm?; theoretisch: I: 100 °K, II: 

473 °K, III: 673 °K, IV: 1200 °K 300 °K und 473 °K, IH: 673 °K, 

IV: 1200 °K 

Abb. 15—18. Anpassung der Meßpunkte an die theoretisch berechneten Kennlinien für 


Kristall 25 


nicht einfach!4), so daß bei mehreren Arten die mathematische Behandlung 
recht schwierig werden dürfte. Solange keine genauen Kenntnisse über die 
Oberflächenzustände vorhanden sind, kann eine Einbeziehung in den quanti- 


tativen Vergleich zwischen Theorie und Experiment nicht erfolgen. = yo 


6.2 Berechnung des Feldes und der emittierenden Fliche = 
Bei der in 5.5 geschilderten Anpassung wurde die Mikrogeometrie durch 
einen Feldverstärkungsfaktor berücksichtigt. Dieser und die emittierende 
Fläche wurde dann so gewählt, daß die Anpassung möglich war. Das gibt dem 


™ Physik. 7. Folge, Bd. 8 
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Verfahren eine gewisse Unsicherheit. Allerdings schränkt die Benutzung des 


Temperaturverhaltens eine willkürliche Parameterwahl stark ein. Bei 


deshall 
Voll befriedigend wäre die Methode erst dann, wenn aus Makro- und Mikro- §, die 


geometrie der Spitze das Feld und die emittierende Fläche berechnet werden § param 
können, oder wenn es gelingt, die Spitze durch Glühen oder Felddesorption f jert w. 


so zu verrunden, daß Gl. (8) mit nur geringem Fehler anwendbar ist. die St 
Die 
6.3 Die Abweichung von der Linearität im unteren Teil der Kennlinie die de 


Der bei niedrigen Spannungen teilweise aufgetretene, relativ flach ver- Feldei 
laufende Teil der Kennlinie (Abb. 2) kann mit den verwendeten Annahmen § Wendi 
nicht quantitativ erklärt werden. Wie verschiedene Abschätzungen ergeben f versch 
haben, können wahrscheinlich auch Oberflächenzustände nicht einen solchen | freie | 
Kennlinienverlauf verursachen. Subenko, Klimin und Sokolskaja®) Fü 
nehmen an, daß Aufladungserscheinungen zu solchen Abweichungen vom fi verwe! 
geradlinigen Verlauf führen können. Da bei unseren Versuchen die Glaswand | Bestät 
ein festes Potential hatte (Abschn. 3), dürften jedoch Aufladungen weitgehend f und e: 
vermieden worden sein. Außerdem war dieser Teil der Kennlinie stets gut | lage : 
reproduzierbar. Fischer!) vermutet, daß eine kapazitive Einstreuung von f Dabei 


Störspannungen für diesen Effekt verantwortlich zu machen ist. angab 
Bu Fi 6.4 Der Knick im oberen Teil der Kennlinie Le 


Bei einem qualitativen Vergleich mit der Theorie wird man zunächst an- 
nehmen, daß die gemessene Änderung der Steigung (Abb. 4) durch einsetzende 
Entartung des Elektronengases hervorgerufen wird. Wie aber Abb. 13 zeigt, 
gehen die Kurven für Entartung und Nichtentartung bei Feldstärken inein- 
ander über, die größer als 10% V/cm sind. Aus der Anpassung ergeben sich aber 
experimentell erreichte Feldstärken von maximal 5 bis 6 - 10° V/em. Demnach 
muß vermutlich eine andere Erklärung für den Kennlinienknick gesucht 
werden. 


Wie in Abschn. 2 angegeben, ist der spezifische Widerstand o der Kristalle 
sehr hoch. Das Abknicken der Kurven könnte also eventuell durch Spannungs- 
abfall am Kristall erklärt werden. Der Widerstand von Kr. 21 (Abb. 4) müßte 
dann 5 - 10°Q2, der von Kr. 37/1 (Abb. 4) 101% Q und der von Kr. 40 (Abb. 4) 

10° sein. Entsprechend den Versuchsbedingungen (o ~ 1022 cm, Kri- 
stallänge 3 bis 4 mm, Kristalldurchmesser 20 um) ergibt sich ein Widerstand 
von 10° 2, so daß die obige Erklärung für die Kennlinienform der Abb. 4 nicht 
abwegig erscheint. 


. Zusammenfassung 


Nach den in *4) BE Beziehungen kann die Feldemissions- Stem: 
dichte j bei Kenntnis der Feldstärke F an der Spitze, der Elektronenaffinität 
y und der Elektronenkonzentration n an der Halbleiteroberfläche berechnet 
werden. In den Versuchen wurde dagegen die Abhängigkeit des Emissions- 
stromes I von der Spannung U bestimmt (bei verschiedenen Temperaturen T). 
Die Löcherkonzentration p, im Halbleiterinneren sowie py = p,(T) waren 
bekannt, desgleichen der makroskopische Spitzenradius. 


3) O. W. Subenko, A.J. Klimin u. I. L. Sokolskaja, ®TT 1, 1845 (1959). 
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Beim Vergleich Theorie-Experiment für eine bestimmte Temperatur müssen 
deshalb Annahmen gemacht werden über die Größe der emittierenden Fläche 
§, die Mikrogeometrie der Spitze (Feldverstärkungsfaktor 6) und y. Die 
Parameter S, 6 und y können jedoch nur innerhalb angebbarer Grenzen vari- — 
iert werden (5.3, 5.5), die Auswahl wird zusätzlich stark eingeschränkt, weil 
die Stromdichte von y und ß exponentiell abhängt. 

_ Die eee des Emissionsstromes wird bestimmt durch | 


freie Parameterwahl weiter ein. | 
Für die Anpassung wurde der mittlere geradlinige Teil der Kennlinien 
verwendet. Es gelang, der angegebenen Grenzen eine 


und experimentelle j-1/F- Kennlinien konnten hinsichtlich Anstieg und sheolatt r 
lage für verschiedene in gebracht 1 werde 


angaben verträglich ist. 


Leipzig, Physikalisches Institut der Karl-Marx-Universität. 
Bir fern Bei der Redaktion eingegangen am 25. Nov ember 1960. 
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= Hochaufgelöste paramagnetische Protonenresonanz- 
spektren einiger Metallalkylverbindungen') 


Von @. Klose 


Mit 6 Abbildungen 1. Th 
Inhaltsiibersicht 

Verhä 


Es werden die hochaufgelösten paramagnetischen Protonenresonanzspek- || proto 

tren von Diäthylcadmium, Tetraäthylzinn, Diäthylselenid und der ent-| ;chiet 
Be Methylverbindungen bei 25 MHz untersucht. Die Spektren von | kann 
_ Diäthylcadmium und Tetraäthylzinn können als Überlagerung der Spektral- | werde 
typen A,B,, A,B,X, A,B,Y (X, Y-Metallisotöpe mit dem Spin !/,) und das äthylı 
Spektrum von Diet ylealenid als A,B,-System interpretiert werden. Die Spek- gieniv 
_ tren der entsprechenden Methylverbindungen entstehen durch Überlagerung } Jct, y 
der Spektraltypen A,, A,X, A,Y und A,Z (Z — Kohlenstoffisotop C1). Die D 
Elektronegativität von Cadmium, Zinn und Selen wird erstens aus der inneren 
chemischen Verschiebung der Athylverbindungen und zweitens aus den OR. 
Kopplungen der Protonen der Methylverbindungen bestimmt. Die exakte 


Lösung der Spin-Hamiltonfunktion für Protonenresonanzspektren des A,B,- mit » 

Typs wird durchgeführt. geme: 

Hierk 
Ma 

Einführung gn 

(CH,- 


Wenige Arbeiten befaßten sich bisher mit der Untersuchung der weit- CH,- 
_reichenden indirekten Spin-Spin-Kopplung zwischen Protonen und anderen entsp 
Kernen. In denerschienenen Arbeiten wurden fastausschließlich Kopplungen zwi- # ;ußer 
schen Protonen und Fluorkernen bzw. Protonen und Phosphorkernen untersucht | zwise 

” a z. B. ?2)?)). Die einzigen uns bekannten Arbeiten, in denen Kopplungskon- A 
iF stanten zwischen Protonen und Kernen mit großen Massenzahlen gemessen f §pjn. 
wurden, sind die Arbeiten von Baker‘), von Raymond, Flautt, Jaffe, | aufse 

_ Reynolds®) und von Narasimhan, Rogers®). Um die weitreichende in- | pejde 
4 = Spin-Spin-Kopplung über zwei bzw. drei Bindungslängen zwischen | Rech 


Protonen und den Isotopen des Cadmiums und Zinns mit dem Spin !/, bzw. } der A 
dem Selenisotop mit dem Spin !/, zu untersuchen, führten wir ein intensives J (.ter 


Er ) Über einige Ergebnisse dieser Arbeit berichtete der Verfasser auf dem Colloque für d 
Ampite 1960 in Pisa. 
2) H.S. Gutowsky, L.H.Meyer und D. W.McCall, J. chem. Physics 28, 98 | Kurz 


(1965). 
2 3) H.S. Gutowsky, C.H. Holm, A.Saika u. G.A. Williams, J. Amer. chem. ?) 
Soc. 79, 4596 (1957). 8) 
a E. B. Baker, J. chem. Physics 26, 960 (1957). (1955) 
be 5) E.D. Raymond, T. J. Flautt, H.H.Jaffe u. G.F. Reynolds, J. chem. 9) 
Physics 80, 1422 (1959). ten K 
8% P. T. Narasimhan u. M.T. Rogers, J. Amer. chem. Soc. 82, 34 (1960). 10 
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_entsprechen der Wechselwirkung der beiden Protonengruppen mit dem 
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Studium der Spektren von Diäthylcadmium, Tetraäthylzinn, Diäthylselenid 
und der entsprechenden Methylverbindungen durch. 

Die Deutung der Spektren der Äthylverbindung des Cadmiums und Zinns 

wird vor allen Dingen durch die ‚geringe natürliche Häufigkeit der Cadmium- 
und Zinnisotope mit dem Spin :/, erschwert. Um die Kopplungskonstanten 
und die chemischen Verschiebungen der untersuchten Äthylverbindungen zu 
erhalten, mußten die Spektren der Äthylgruppen theoretisch berechnet werden. 
Wir führen deshalb im ersten Teil dieser Arbeit die exakte Lösung der Spin- 
Hamiltonfunktion für Protonenresonanzspektren des A,B,-Typs durch. 


1. Theoretische Berechnung eines Protonenresonanzspektrums vom AgBs>-Typ 


Bekanntlich hängt das Aussehen des Spektrums einer Äthylgruppe vom 
Verhältnis der indirekten Spin-Spin-Kopplungskonstanten J zwischen den 
Protonen der CH,- und CH,-Gruppe zur relativen (inneren) chemischen Ver- 
schiebung 6 dieser beiden Gruppen ab. Ist dieses Verhältnis klein gegen 1, so 
kann die einfache Störungstheorie zur Berechnung des Spektrums benutzt 
werden”). Ist diese Bedingung nicht erfüllt, wie es bei den Spektren von Di- 
äthylcadmium und Tetraäthylzinn der Fall ist, muß zum Auffinden der Ener- 
gieniveaus und relativen Intensitäten die zugeordnete Säkulargleichung ge- 
lost werden. 


X =, 1,(A) 1,(B) B) 


Y 


yH H ‘ 
4 und 9, = 3 = geschrieben werden, wenn die Energie in Hertz 


2n 
gemessen wird und das konstante Magnetfeld in negative z-Richtung zeigt. 
Hierbei ist y das gyromagnetische Verhältnis der Protonen, H, bzw. H, das 
Magnetfeld am Ort der Kerngruppe A (CH,-Gruppe) bzw. der Kerngruppe B 
(CH,-Gruppe), und J,(A) und J,(B) sind die Spinoperatoren für die CH,- und 
CH,-Gruppe in z-Richtung. Die ersten beiden Terme der Hamiltonfunktion 


mit = 


äußeren Magnetfeld, während der dritte Term die Spin-Spin-Wechselwirkung 
zwischen den beiden Gruppen berücksichtigt. 

Als Eigenfunktionen 0-ter Ordnung könnten die 32 einfachen Produkt- 
Spin-Eigenfunktionen benutzt werden, die man für 5 Kerne mit dem Spin !/, 
aufschreiben kann, und die bei verschwindendem Wechselwirkungsterm der 
beiden Protonengruppen die stationären Eigenfunktionen darstellen. Die 
Rechnung kann jedoch wesentlich vereinfacht werden durch Berücksichtigung 
der Äquivalenz der Protonen der CH, bzw. CH,-Gruppe in den Spinfunktionen 
0-ter Ordnung. Zum Aufschreiben dieser Wellenfunktionen verwenden wir 
die von McConnell, Mc Lean und Reilly®) abgeleiteten Spinfunktionen 
für die symmetrisch) äquivalenten Kerne der CH,- und CH,-Gruppe _ der 
Kurzbezeichnung von Whitman”) (Tab. 1). a: 


?) J.T. Arnold, Physic. Rev. 102, 136 (1956). ne Oe 
Mi. H.M. McConnell, A.D.McLean u. C.A. Reilly, J. chem. Physics 23, 1152 
(1955). 
®) Symmetrisch heißt, daß der Hamiltonoperator bei V ertauschung der äquivalen- 
ten Kerne der CH,- bzw. CH,-Gruppe ungeändert bleibt. u 
10) D.R. Whitman, Dissertation, Jale Universität 1957. 
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Mit den 8 Funktionen für die CH,-Gruppe und den 4 Funktionen der CH,. 
Gruppe erhalten wir die 32 Produkt. Spin- Funktionen 0-ter Ordnung für én 
A,B,-System. Diese Funktionen wurden in der Tab. 2 mit den berechneten 
diagonalen und nichtdiagonalen Matrixelementen aufgeschrieben. Wie aus 
| der Tab. 2 unmittelbar abgelesen werden kann, zerfällt die Säkulargleichung 


Tabelle 1 
ie ar Spinfunktionen 0-ter Ordnung für die symme- 
us trisch äquivalenten Kerne der CH,- und CH,- 
“ere Gruppe 
Bezeich- 
Funktion nung 


CH,-Gruppe axa 
(aa B+ apa + Bax) 


(«aa B+aBba —2B ax) 


(BBa+BaB+a 88) 
er 
ze 2a BB) 


CH,-Gruppe «« 


1 
5 + Ba) 


T_, 
in zwei 3-3, sechs 2-2 und zwölf 1-1-Determinanten. Davon sind die Deter- 
minanten, die den D’T- und D’S,-Funktionen entsprechen, gleich den Deter- 
minanten für die DT- und DS, „Funktionen. Sie brauchen also bei der Berech- 
nung der Übergangsfrequenzen nicht berücksichtigt zu werden, während sie 
bei den Intensitäten der DT- und DS,-Ubergange durch einen Faktor 2 ihre 
Berücksichtigung finden. 

In der Tab. 3 sind die durch Lösung der Säkulargleichungen erhaltenen 
Energieniveaus und die exakten Wellenfunktionen zusammengestellt worden. 
Zur Vereinfachung der Schreibweise wurden dabei folgende Abkürzungen ein- 


[Fo +o +1] s(t w+o+1]", 


1/2 


j 99 
We 
n 
1 
| 
3 
_ 
= 
16 
17 
49 
| 
A 20 
21 
i) 
22 
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Tabelle 2 
Wellenfunktionen 0-ter Ordnung und Matrixelemente für den Fall A sB, 
diagonales nichtdiagonales 
» Funktion I, Matrixelement Matrixelement 
5 3 3 
= 
1 Q, T. | +, +4 5 J 
| | 3 
2 Q% T, 3 "4 3 
= 1 (®, = —— 
3 7; | v6 
2 2 
4 Q-ı - Zu tn -4J 
‘ 1 1 (®,|X|®,) = 
9 Q, Ty, | „ra V6 
3 (®; |X| ,) 
6 T-, | | 37 
3 
7 | - Zn tm 
1 1 (®, |X| ®,) = J 
8 | Q47, ‘tee | V4 
1 
10 Q-, 3 - — Vz 3 
3 (Po |X| = J 
11 Q-,T, V6 
3 
12 QT. | -z 
3 3 
| 
nd sie 3 1 
i DT — + +—J 
| 1 1 
ltenen 18 D_, T, 1 trp — 
orden 7 (Py, |X|.) = — 
n ein 1 To 
1) | 
9% 3 
2 | | 


2 


> 
> 
ag 
{ a 
f 
2 


Tabelle 2 (Fortsetzung) 
Funktion I, Matrixelement Matrixelement 
1 
23 D, Sy 4 
vy - 
1 
24 D_ So 7 V4 
| 
3 1 1 
25 D; T, Va Vp J 
1 1 
26 1 
27 | mm 2 4, 2 
1 
28 D,T, Fr 1 
: = 1 (Pog |X| Dog) = — > J 
29 Di T., 2 


a Die eingeführten Winkel werden durch die Beziehungen 


3 11 13 15 19 1 15 
u +(20 (jo +75) % 7? 0 
und 


3 11 13 15 19 1 15 
+(20 +7)e- (Fe 0 


wobei u, < u, < u, und w, < w, < wz gilt. Die Koeffizienten a,, b,, C,,; Gn On 
und ¢, mit n = 1, 2, 3 werden aus den Gleichungssystemen Daa ge 


= 0, 


G. 
Ta 
1 
1 
1 
1 
Gis 
gegebe und w, mit ind die Wurzeln der Gleichungen a. 
18 
2¢ 
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Tabelle 3. Energieniveaus und exakte Wellenfunktionen für den Fall A,B, 
— n Energie | Wellenfunktion 
5 3 
3 
2 +s Q; T, sin 6_ + Q, 7, cos 
3’ Q; T, cos 6. — Q, T, sin 
2 6 
4’ a, Q-, Tı + To +, T-ı 
1 é 
2 
A 1 ö 
1 ö 
1 ö 
1 ö = 
9’ +7+%9 Gs Q,T_, + 6; Q-, To, + C3 Ty 
ö 1 as 
11 6 Q_, T_, cos 6, — Q_, Ty sin 64 
12’ Q., 
13’ Q3 So 
14’ Q Sp 
15’ Q-1 So 
16’ Q-, So 
17’ > DT, 
0 18’ D_, T, cos ys — D, T, sin y+ 
19’ D_, sin + D, T, cos y+ 
1 6 1 
20’ ——v, 
b 2 2 4 
ans n ö 
21’ = 
22” 
23’ 


a 
| 
| 1 
24’ | D_, So 


E- 226 Annalen der Physik. 7. Folge. Band 8. 1961 6. K 
Tabelle 4 
Übergangsfrequenzen und relative Intensitäten der Absorptionslinien für den iy 
ABs sie) 
Ubergangsfrequenz 
Linie Ubergang rel. zu v, in Einheiten relative Intensität a ue 
von 6 
93 | 11’ 
16’ — 15’ (QS) 2 94 | 11’- 
24’ > 23° (DS) 0 9 95 | 11". 
1’ (QT) -i+ +C {3 sin 6_ — V2 cos 6_}? 96 | 
| 
1 QT) ->+70 {V2 sin 6_ + V3 cos 6_}? alr 
1 5 2 = 
PM) | -—-ze+t {V3sin 6, + V2 cos 0,}? og | 7 
12’ > 11’ (QT) -4-20-0 {V2 sin 0; — V3 cos | 8° 
18’ > 17’ (DT) 2 {/2siny, + cos ys}? g 
| 1 3 2 = 
2 {sin yp, — V2 cos 8) 
| 29° 20’ (DT = 1 3 9 {si ya 
| 22’—21’(DT) -4-2o+F | 2 {/2siny- + cosy_}? 9 
| | | 
0>18(DN| —1—-F'+F 2 {V2 cos (y_ + yx) + cos y_ sin y,}* 
20’>19 (DT)| —1+F+4+F 2 {/2sin (y_ + yx) + cos y_ cos y,}* 
| 21’ +18’ (DT)| —-ı-F-F 2 {/2sin (y_ + y+) + sin y_ sin y,}? 
| 21’+19(DT)| —1+F-—F 2 {V2 cos (y_ + y+) + sin y_ cos u 
| 
| 3(Q7) {sin 6 (2a, + b,) — cos (V3 + 
(QT)| — + + —u,—C | {sin (2a, + /2 — cos (V3 5, + V2 
| | {sind (2a, + — cos (VB by + V2 
2’ (QT) {cos (2a, + b,) + sin (V8 5, + 2 «)}? 
2° (QT) {cos 6_ (2 a, + V2b,) + sin 6_ (V3 b, + V2 
| 0 
2’ (QT) {cos 6_ (2a, + + sin 6_(Y3 + V2 «)} 
10° > 9 (QT) | — + _ + wo +w,— C’| {sin 0, (2a, + Y2b,) + cos 0, + 2 
. 
1’> 8(QT) | {sin 8, (2a, + + cos 6, (V3B, + 
10°> 7(QT) | — + jo {sin 0, (2a, + /26,) + cos 0, (V3 b, + y22)} 
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Tabelle 4 (Fortsetzung) 
den } Übergangsfrequenz 
wil inie bergang rel. zu v, in Einheiten a relative Intensität Ba BE 
sg 
—3 | {cos 6; (2, + — sin 0, (V35, + 
| 1 1 Jak > 2 
8 (QT) | -—— {cos 0, (2a, + V26,) — sin 0, (//3 6, + V/2<,)}? 
%  1l'> 7 (QT) | + _ + w, + | {cos 6, (2a, + V25,) — sin 0, (//3 6, + 
1 
4 (QT) {V3q@at+e a) + (a,b, + a,) + 20, b,}? 
| | 
| — + Uy — {V3 (a3 + + V2 (a,b, +b, a) + 2b, 
| 1 2 
B(QT)| — + Us — {V3 (ast, + + V2 (a,b, + + 2 by b,}? 
1 9 
| {V3 (a, & + Gq) + V2 (a, b, + + 2b, bg} 
{V3 a) + V2 (a, by + by + 2b, bg}? 
0 | 5' (QT) t% ws VB (ag Cy + Cy + V2 (ag by + by ay) + 2b, dg 
| 1 
a|8> 6 QT) | - {VB (ay Ce + Gy) + V2 (ay by + bs + 2by b,}* 
| | 
32 | 4° (QT) - | {VB (ay + + V2 (a, + +26, 
| | 
| 1 y le 
}33 | 9> | — {V3 (gt + 6303) + V2 (ay by + by ag) + 2 by bg}? 
1 - 2 
34 | I> {V3 (ag Cy + Cy + V2 (ag by + by + 2b, 


q 
erhalten. 
Bei Beachtung der Auswahlregeln, daß Übergänge nur zwischen Energie- 
niveaus, die zur gleichen Symmetriegruppe gehören, und für AF, = + 1 mög- 
lich sind, erhalten wir die Übergangsfrequenzen. In Tab. 4 sind die Übergangs- 
frequenzen relativ zu v, in Einheiten der chemischen Verschiebung 6 mit den 
entsprechenden relativen Intensitäten zusammengestellt worden. Von den 34 
Linien sind die Linien 12, 16, 17, 22, 23, 27, 28, 31 und 32 so intensitätsschwach, 
daß sie nicht beobachtet werden können. Sie werden daher bei der Berechnung 
der Spektren von Diäthylcadmium, Tetraäthylzinn und Diäthylselenid nicht 
berücksichtigt. 
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2. Protonenresonanzspektren von Diäthyleadmium und Tetraäthylzinn grup 

B Die eg re von Diäthylcadmium (Abb. 1, oben) und und 

Tetraäthylzinn (Abb. 2, oben) bestehen aus einigen intensiven Mittellinien nr 

Pr. mit mehreren schwachen Satellitenlinien. Die Satellitenlinien sind das Ergebnis ui 

der ungleichen indirekten Spin-Spin-Kopplung der Cadmium-Isotope Cd! or d 
(Spin natürliche Häufigkeit 12,75%) und (Spin !/,, natürliche Häufig- 


> 


Tu m, 


Abb. 1. Experimentelles und berechnetes Spektrum von 
Diäthylcadmium bei 25 MHz 


nic 

der 

keit 12,26%) bzw. der Zinnisotope Sn!!? (Spin !/,, natürliche Häufigkeit 7,5%, liel 

und $nlıs (Spin !/,, natürliche Häufigkeit 8,6%,) mit den Protonen der CH,- der 

und CH,-Gruppe. Nach der Störungstheorie 1. Ordnung sind für die Kopplung abr 
5 eines Isotops mit den Protonen der Athylgruppe entsprechend eines A,B,X- Ve 
_ Systems ein Quadruplett und ein Triplett auf jeder Seite des Spektrums des Ke 

_ A,B,-Systems zu erwarten®). Da sowohl Cadmium als auch Zinn zwei Isotope Lit 
2 mit dem Spin !/, und ausreichender Häufigkeit besitzen, müssen also bei ZW 
genügend großen Kopplungen der Metallisotope mit den Protonen der Äthyl- lit 


4 
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gruppe entsprechend der Überlagerung der Linien der Spektraltypen A,B,X 
und A,B,Y (X und Y kennzeichnen die beiden Cadmium- bzw. Zinn-Isotope) 
zwei Tripletts und zwei Quadrupletts auf jeder Seite auftreten. | 

Das experimentell erhaltene Spektrum von Diäthylcadmium zeigt statt 
der erwarteten Zahl von Multipletts nur drei Linien auf der Seite höheren 
Feldes, die wohl die für Tripletts charakteristischen gleichen Abstände aber 


Abb. 2. Experimentelles und berechnetes Spektrum von 
Tetraäthylzinn bei 25 MHz 


nicht das einfache 1:2:1-Intensitätsverhältnis besitzen, und drei Linien auf 
der Seite schwächeren Feldes. Beim Spektrum von Tetraäthylzinn sind deut- 
lich Satellitenlinien nur auf der Seite höheren Feldes zu erkennen, die wohl 
den für ein Quadruplett charakteristischen gleichen Abstand voneinander 
aber nicht das 1:3:3:1-Intensitätsverhältnis zeigen. Diese unübersichtlichen 
Verhältnisse kommen dadurch zustande, daß erstens infolge der „kleinen“ 
Kopplungskonstanten zwischen den Metallisotopen und den Protonen einige 
Linien der Multipletts in den intensiven Mittelteil des A,B,-Spektrums fallen; 
zweitens durch die geringe natürliche Häufigkeit der Isotope werden die Satel- 
litenlinien von den schwachen Seitenlinien des A,B, EEE überdeckt; 
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drittens infolge des nahe bei 1 liegenden Verhältnisses der Kernmomente der 
beiden Cadmium- bzw. Zinnisotope erscheinen die Linien des A,B,X-Spektrums 
nicht von den Linien des A,B,Y-Spektrums getrennt, und daß viertens die 
Linien der Tripletts und Quadrupletts übereinander fallen. Bei beiden Spektren 
wurde angenommen, daß die beobachteten Satellitengruppen auf der Seite 
höheren Feldes durch Überlagerung der erwarteten Quadrupletts und Tri- 
pletts entstehen. Diese Annahme wird durch das gemessene 4:9:6-Intensitäts- 
verhältnis beim Diäthylcadmium und das beim Tetraäthylzinn gemessene 
4:6:4:3-Intensitätsverhältnis bestätigt. 

Eine entsprechende einfache Bestimmung der Multipletts auf der Seite 
schwächeren Feldes ist infolge der Kompliziertheit der Spektren unmöglich. 
Um in diesen Fällen die Kopplungskonstanten der Metallisotope mit den 
Protonen zu erhalten, müssen die Protonenresonanzfrequenzen der CH,- und 
CH,-Gruppe v, und v, des A,B,-Spektrums bestimmt werden, daß die an die 
Metallisotope mit dem Spin 0 bunden Äthylgruppen ergeben. Diese Spek- 
tren wurden mit Hilfe der Formel der Tab. 4 berechnet. Die Kopplungskon- 
stanten zwischen den Metallisotopen und den Protonen der CH,-Gruppe J , , 
und J ,y ergeben sich dann als der doppelte Abstand der Mittellinien der Satel- 
litentripletts auf der Seite höheren Feldes von »,, und entsprechend sind die 
Kopplungskonstanten zwischen den Metallisotopen.und den Protonen der CH,- 
Gruppe J,y und J,y gleich die doppelte Differenz zwischen dem Zentrum 
der Quadrupletts und »,. Die bei den Mittellinien der Satellitentripletts vor- 
handene Aufspaltung (s. Abb. 1 und 2) ermöglicht die Bestimmung der Kopp- 
lungen beider Cadmium- bzw. beider Zinnisotope mit den Protonen der CH,- 
bzw. CH,-Gruppe. 

Die als Übe.lagerung der Spektraltypen A,B,, A,B,X und A,B,Y inter- 
pretierten und berechneten Spektren von Diäthylcadmium und Tetraäthyl- 
zinn wurden unter die entsprechenden experimentellen gezeichnet (s. Abb. 1 
und 2, unten). Die zu den verschiedenen Multipletts gehörenden Satelitten- 
linien wurden durch gestrichelte Linien verbunden und die Kopplungskon- 
stanten eingezeichnet, um die Übersicht zu erleichtern. Die recht gute Über- 
einstimmung zwischen den experimentell erhaltenen und den berechneten 
Spektren ist offensichtlich. Die für die Kopplungskonstanten und chemischen 
Verschiebungen erhaltenen Werte wurden in der Tab. 5 zusammengestellt. 


3. Protonenresonanzspektrum von Diäthylselenid 


Das Spektrum von Diäthylselenid (Abb. 3, oben) zeigt ein wesentlich 
anderes Aussehen als die Spektren von Diäthylcadmium und Tetraäthylzinn. 
Das ist darauf zurückzuführen, daß das J/ö-Verhältnis im Vergleich zu den 
Werten bei den Äthylverbindungen des Cadmiums und Zinns klein ist und 
daß keine Satellitenmultipletts auftreten, obwohl infolge des Vorhandenseins 
des Selenisotops Se’? mit dem Spin !/, und der natürlichen Häufigkeit von 
8,3%, welche erwartet werden. 

Die Berechnung des Spektrums als A,B,-Typ ergab, daß alle beobachteten 
Linien zugeordnet werden konnten mit Ausnahme der intensitätsschwächsten 
Linie im ,,Triplett*‘ der CH,-Gruppe. Es wird vermutet, daß diese Linie auf 
die Kopplung des Se’’-Isotops mit den Protonen der CH,-Gruppe zurückzu- 
führen ist. 
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von Dimethylcadmium, Tetramethylzinn und Dimethylselenid 


Die Spektren von Dimethylcadmium (Abb. 4), Tetramethylzinn (Abb. 5) 
und Dimethylselenid (Abb. 6) bestehen aus einer intensiven Mittellinie, die Ets 


| | 
M i Il | 
Abb. 3. Experimentelles und berech- 
netes Spektrum von Diäthylselenid 


bei 25 MHz 


=» 


von den Protonen der an die Metallisotope mit dem Spin 0 gebundenen CH,- 
Gruppen herrühren (A,-System), und zwei symmetrisch zur Mittellinie gelege- 
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Abb. 4. Spektrum von Tetramethylzinn bei 25MHZ 
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nen schwachen Satellitenlinien, die durch die indirekte Spin-Spin-Wechsel- 
wirkung der Metallisotope mit dem Spin !/, zustande kommen (A,X-System). 
Beim Tetramethylzinn sind diese Satellitenlinien in zwei Linien aufgespalten 
entsprechend dem Vorhandensein der zwei Zinnisotope Sn!!? und Sn}, 
während die beim Spektrum von Dimethyleadmium entsprechend dem Vor- 
handensein der Cadmiumisotope Cd!!! und Cd!!3 erwartete Aufspaltung der 
Satellitenlinien nicht beobachtet wurde. Außer diesen Satellitenlinien werden 
bei allen drei Spektren zwei weitere sehr schwache Linien beobachtet, die eben- 


2 


| 


Abb. 5. ee von Di- Er by Abb. 6. Spektrum von Di- 
methylcadmium bei 25 MHz 7 methylselenid bei 25 MHz 


falls symmetrisch zur Mittellinie gelegen sind und durch die indirekte Spin- 
Spin-Kopplung der Protonen mit dem Kohlenstoffisotop C!3 (Spin 1/,, natür- 
liche Häufigkeit 1,1%) entsprechend einem A,Z-System zustande kommen. 
Die weiteren schwachen Linien im Spektrum von Dimethyleadmium rühren 
von unbekannten Verunreinigungen her. 


Die aus den Spektren erhaltenen Kopplungskonstanten zwischen den 
Protonen und den Metallisotopen bzw. den Protonen und dem C!*-Isotop wur- 
den ebenfalls in der Tab. 5 zusammengestellt. 


Die Werte der bestimmten Kopplungskonstanten und chemischen Ver- 
schiebungen der untersuchten Äthyl- und Methylverbindungen des Cadmiums, 
Zinns und Selens wurden in der Tab. 5 zusammengestellt. Ein Vergleich der 
Kopplungskonstanten JMeH(CH,)-Me und JMeH(CH,)-Me der Äthylverbin- 
dungen zeigt, daß die Kopplung der Protonen der CH,-Gruppe mit den Metall- 
isotopen größer ist als die entsprechende Kopplung der Protonen der CH,- 
Gruppe. Dieses Ergebnis läßt sich nicht mit den Vorstellungen über den Mecha- 
nismus der weitreichenden Kopplung über mehrere Bindungen verstehen. Es. 
muß ein zusätzlicher Mechanismus wirksam sein, der wahrscheinlich in einer 
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unmittelbaren Kopplung zwischen dem Metallatom und den Protonen der 
H 


Den Beitrag, den dieser Mechanismus zur Kopplungskonstanten liefert, 
ist um so größer, je größer die Ordnungszahl des Me-Atoms im periodischen 
System oder, mit anderen Worten, je größer die Elektronenzahl des Me- 
Atoms ist. 


Aus den bestimmten inneren chemischen Verschiebungen 6 berechneten wir 
die Elektronegativitäten der Metalle nach der von Dailey und Shoolery empi- 
risch gefundenen Gleichung !?) 


e, = 0,01929 6 + 1,71.) 


Für Atome bzw. Molekülgruppen mit kleiner Elektronegativität muß die 
Konstante von 1,71 in der ursprünglichen Gleichung nach Shoolery!#) durch 
2,1 ersetzt werden. Sowohl die mit einer Konstanten von 1,71 (e,) als auch die 
mit einer Konstanten von 2,1 berechneten Elektronegativitäten (e;) wurden 
in der Tab. 6 zusammengestellt. Wenn wir annehmen, daß Elektronegativi- 
täten kleiner als 2 klein sind, dann erhält man für Zinn und Selen eine gute 
Übereinstimmung der nach der Dailey-Shoolery-Gleichung berechneten 
Elektronegativitäten mit den von Pritchard und Skinner als ,,beste“ 
Werte für die Elektronegativitäten angegebenen Daten") (e), während für Cad- 
mium ein zu großer berechneter Wert erhalten wird. 


as Tabelle 6 
Elektronegativitat von Cadmium, Zinn und Selen 
Element r(C—X) Juc13 | E, e 
ae 2,18 126,2 1,3 1,7 1,45 1,85 1,4 
- Sn 2,18 127,0 1,5 1,9 1,49 1,89 1,8—1,9 
u é Se 1,977 140,2 2,4 2,8 2,43 2,51 2,4 


Außer der Bestimmung der Elektronegativitäten aus den inneren chemischen 
Verschiebungen der Äthylverbindungen führten wir die Berechnung der Elek- 
tronegativitäten aus den Kopplungskonstanten zwischen den Protonen und 
dem C3.Isotop der Methylverbindungen erstens nach der Formel von Muller 
und Pritchard’) 


Juc» = 22,6 E, + 40,1r(C — X) +5, 


12) B.P. Dailey u. J. N. Shoolery, J. Amer. chem. Soc. 77, 3977 (1955). 

13) In der in®) angegebenen Gleichung von Dailey und Shoolery beträgt der Faktor 
vor Ö 0,02315 entsprechend einer Arbeitsfrequenz von 30 MHz. 
14) J.N.Shoolery, Lecture Notes, Varian Associates Conference (Nov. 1957). 
35) H.O. Pritchard u. H. A. Skinner, Chem. Rev. 55, 745 (1955). 
15) N. Muller u. D. E. Pritchard, J. chem. Physics 31. 1471 (1959). 


Th 


dur 
Kol 
tro! 
recl 
„In 
miu 
Kol 
suc 
Ele 
wat 
die 
gen 
ent 
bzv 
auf 
stä 
höl 
wa 
deı 
baı 
Ve: 
ent 
| Da 
| un 
die 
Tä 
erf 
ke 
tol 
ge; 
an 
sel 
Cy 
At 
40 
= 
| 
¢ 


der 


aktor 


und zweitens nach der von Ranft angegebenen Formel!?) 
rs 


durch, wo r (C— X) die Bindungslänge zwischen dem Met FRE und dem 
Kohlenstoffatom ist. Die nach diesen Formeln erhaltenen Werte für die Elek- 
tronegativitäten wurden ebenfalls in die Tab. 6 aufgenommen. Die zur Be- 
rechnung benutzten Bindungslängen für Zinn und Selen wurden dem Buch 
„Interatomie Distances‘ 1?) entnommen. Die Bindungslänge zwischen Cad- 
mium und Kohlenstoff wurde aus den covalenten Radien von Cadmium und 
Kohlenstoff nach 1?) berechnet. 

Das Verhalten der gemessenen chemischen Verschiebungen der unter- 
suchten Alkylverbindungen läßt sich qualitativ mit den Vorstellungen über die 
Elektronegativitäten der an den Bindungen beteiligten Atome erklären. Mit 
wachsender Elektronegativität in der Reihenfolge Cadmium, Zinn, Selen nimmt 
die Elektronendichte in der Umgebung der CH,-Gruppe der Methylverbindun- 
gen und in der Umgebung der CH,-Gruppe der Äthylverbindungen ab; und 
entsprechend tritt eine Verschiebung der CH,-Linien der Methylverbindungen 
bzw. der CH,-Linien der Äthylverbindungen in Richtung schwächeren Feldes 
auf. Dabei liegen die CH,-Linien der Methylverbindungen bei höherer Feld- 
stärke als die CH,-Linien der entsprechenden Äthylverbindungen infolge der 
höheren Elektronendichte in der Umgebung der CH,-Protonen. Wie |zu er- 
warten ist, unterscheiden. sich die chemischen Verschiebungen der CH,-Protonen 
der Athylverbindungen 6, wenig voneinander, während die 6, durch die Nach- 
barschaft der Metallatome stark beeinflußt werden. Dabei sind die chemischen 

Verschiebungen 6, von Diäthylcadmium und Tetraäthylzinn kleiner als die 
entsprechenden ö,-Werte, während beim Diäthylselenid 6, größer als 6, ist. 
Das ist darauf zurückzuführen, daß die Elektronegativitäten von Cadmium 
und Zinn kleiner als die Elektronegativität des Kohlenstoffs (2,55) sind und 
die Elektronegativität von Selen ungefähr gleich der des Kohlenstoffs ist. 


C-X 


6. Experimentelles 


Die Spektren der untersuchten Substanzen wurden mit einem Triib-, 
Täuber-Spektrometer bei 25 MHz aufgenommen. Die Ausmessung der Spektren 
erfolgte mit Hilfe des bekannten Seitenbandv erfahrens. Die Frequenzgenauig- 
keit des zur Erzeugung der Seitenbänder verwendeten Niederfrequenzgenera- 
tors betrug 1%, woraus die relativ hohe Meßungenauigkeit der in Tab. 5 an- 
gegebenen Werte resultiert. 

Die Kopplungskonstanten und die innere chemische Verschiebung wurden 
an unverdünnter Substanz gemessen. Zur Bestimmung der chemischen Ver- 
schiebung der CH,-Gruppen in bezug auf Cyclohexan wurde ungefähr 5°, 
Cyelohexan der reinen Substanz als innere Referenzsubstanz zugegeben mit 
Ausnahme von Dimethylcadmium, dessen chemische Verschiebung ö, in einer 
40proz. Lösung von Dimethyleadmium in Cyclohexan bestimmt wurde. 


) J. Ranft, Vortrag zur Arbeitstagung Hochfrequenzspektroskopie, Leipzig 1960. 
18) Tables of interatomic distances and configuration in molecules and ions, London, 
The chem. Soc., Burlington House, W. 1, 1958. 
19) V,Schomaker u. D. P. Srevenson, J. Amer. chem. Soc. 63, 37 org 
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Zur Aufnahme der Spektren wurden Proben mit einem Volumen von 
2—4 mm? verwendet. Die Proben wurden nicht besonders entgast oder abge- 
schmolzen mit Ausnahme der Dimethyl- und Diäthylprobe des Cadmiums, die 
infolge der Zersetzlichkeit der Äthyl- und Methylverbindung des Cadmiums 
an der Luft unter einer Stickstoffatmosphäre gefüllt und abgeschmolzen 
wurden. 

Die verwendeten metallorganischen Verbindungen wurden am Organisch- 
Chemischen Institut der Universität Leipzig hergestellt. Herrn Dr. Weißen- 
fels, unter dessen Anleitung die Substanzen angefertigt wurden, danke ich 
für sein freundliches Entgegenkommen in bezug auf alle unsere Wünsche 
sowie für seine Mühe bei der Beschaffung der Ausgangssubstanzen. 


Herrn Prof. Dr. Lösche bin ich für sein stetes förderndes Interesse an 
dieser Arbeit und Herrn Dipl.-Phys. Ranft für wertvolle Hinweise zu Dank 
verpflichtet. 

Leipzig, Physikalisches Institut der Universität... 00... 


ur Bei der Redaktion eingegangen am 8. Dezember 1960. J oe 
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Neu erscheint 


Die physikalischen Grundlagen der Elektroakustik 


Von Prof. Dr.-Ing. WALTER REICHARDT 
Etwa 190 Seiten mit etwa 110 Abbildungen. 8°. Geb. etwa 12,— DM 
(Mathematisch-Naturwissenschaftliche Bibliothek Band 31) 
Inhalt: Die Probleme der Elektroakustik. Die elektromechanischen Analogien. Das 
Schallfeld. Ankopplung des Schallfeldes an einen mechanischen Schwinger. Das BaBreflex- 
gehäuse als Beispiel für die Kombination zweier Analogien und der Kopplung mit dem 
Schallfeld über zwei strahlende Flächen. Elektromechanische Wandler. Schallwandler. 
Das Buch soll diejenigen Grundlagen zusammenfassen, die gebraucht werden, um elektro- 
akustische Probleme quantitativ zu erfassen, d. h. die physikalischen Voraussetzungen 
und ihre mathematische Einkleidung klarzustellen. Den elektromechanischen Analogien, 
der Vierpoldarstellung von mechanischen Schwingern und der Zweipoldarstellung eines 


angekoppelten Schallfeldes wurde dabei besondere Aufmerksamkeit gewidmet. 


2 a Bestellungen nimmt jede Buchhandlung entgegen. 


B. G. Teubn 


Bruno H. Bürgel zum Gedenken 


Wir suchen ständig 


Bürgel-Arbeitskreises im Deutschen Kul- 
turbund, Kreis Potsdam, von L. Wörner, 
P. Nell t¢, J. Deutsch und W. Tripmacker. 


Herausgegeben im Auftrag des Bruno-H.- . 
4 


3. Auflage. 
92 Seiten mit 1 Bildnis, 8 Tafeln und 
2 Textzeichnungen. kl. 8°, Leinen DM 4,— 


Wer war Bruno H. Bürgel, was bedeutet 
er uns? Eine Antwort gibt dieses Buch 
mit seinen Beiträgen und Werkproben. 
Freunde und Verehrer rufen die Erinne- 
rung wach und errichten uns dies Bild 
eines vorbildlichen, aufrechten, der Zu- 
kunft zugewandten Humanisten. 

Vom Hilfsarbeiter steigt Bürgel auf zum 
Schriftsteller, der unzähligen Tausenden 
den Weg zum Wissen um die Naturkräfte 
und zu einer lichten Zukunft weist. Wie 
erunbeirrt und zäh, dem gesellschaftlichen 
Fortschritt zugewandt, der Verwirklichung 
seiner Lebensaufgabe zustrebt, das ver- 
dient unsere Hochachtung; wegen seiner 
Güte, seines tiefen Verständnisses für 
menschliche Nöte verehren wir ihn; und 
seine sprachlich überzeugende. eindring- 
liche Darstellungsart, ganz auf die Bildung 
des Lesers, des einfachen Menschen ge- 
richtet, ist und bleibt beispielhaft für 
populärwissenschaftliche Gestaltung. 


JOHANN AMBROSIUS BARTH /VERLAG / LEIPZIG 


Bücher 
und 
Zeitschriften 


aus allen Gebieten der Natur- und 
Geisteswissenschaften sowie 
alte und neue Grafik 
und erbitten Angebote von Einzel- 
werken und Bibliotheken. Gleichzeitig 
bemühen wir uns um Beschaffung 
wissenschaftlicher Literatur sowie 
grafischer Blätter und erwarten Ihre 
Anfragen. 


Zentral- Antiquariat der Deutschen 
Demokratischen Republik 
Leipzig C1, Talstr. 29 
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 Quantentheorie des Atoms 
4: 
Von Prof. Dr. A. RUBINOWICZ, Warschau 
Vonenng nach der 2. polnischen Auflage des Werkes ,,Kwantowa Teoria 
 Atomu“. VII, 486 Seiten mit 38 Figuren im Text. 1959. gr. 8°. 


Me Broschiert DM 29.80, Leinen DM 31.80 


Dieses Lehrbuch ist aus den langjährigen Vorlesungen des Verfassers hervorgegangen. Im 
> essten Teil des Buches wird die Entwicklung der älteren Quantentheorie dargestellt. Es 
ots” ir a, sdli damit dem Leser vor Augen geführt werden, wie zwangsläufig und wie reich an Momenten 
‘og ® » "voll ‘dramatise her Spannung die Entwicklungsgeschichte der Quantentheorie gewesen ist. 


. Im zweiten Teil werden die Grundlagen und Anwendungen der neuen Quantentheorie 
re gebracht. Dem Kapitel II ist ein mathematischer Anhang beigegeben, der die Sommer- 
Sdi@iche Polynommethode zur Lösung der Eigenwertprobleme der Quantentheorie in 
Auer allgemeinen Fassung darstellt. Die Aufgabensammlung am Schluß des Buches hat 
den Zweck, dem Leser die Möglichkeit zu geben, seine Kenntnisse zu überprüfen und 
sich mit ihren Anwendungen vertraut zu machen. Außerdem-enthält sie einige spezielle 
im Buchtext nicht berücksichtigte Tatsachen und Überlegungen. 


Einführung in die Theorie der Elektronenoptik 


er Von Dr. phil. JOHANNES PICHT 


0. Prof. für theoretische Physik und Optik, Pädagogische Hochschule 
ae Potsdam-Sanssouei (früher Techn. Hochschule Berlin-Charlottenburg) 
2., erweiterte Auflage 
VII, 274 Seiten mit 70 Abbildungen im Text. 1957. gr. 8°. 

: Broschiert DM 29,20, Leinen DM 31.— 


Die neue Auflage ist auf den neuesten Stand der Wissenschaft gebracht worden, wodurch 
sich eine Erweiterung des Umfanges ergab. Der Verfasser war bemüht, den Charakter einer 
Einführung in das Gebiet zu wahren, insbesondere die mathematischen und gedanklichen 
Überlegungen soweit durchzuführen wie dies zum ausreichenden Verständnis der theoreti- 
schen Zusammenhänge erforderlich schien. 


JOHANN AMBROSIUS BARTH / VERLAG / LEIPZIG 
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